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The  purpose  of  this  research  is  to  develop  algorithms  for  mapping  terrain 
characteristics  from  a  100  meter  grid  representation  to  arcs  and  nodes  of  a 
transportation  network.  The  algorithms  capture  elevation  and  trafficability  parameters 
as  they  relate  to  both  the  arc  itself  and  to  the  ofF-arc  characteristics.  These  network 
parameters  are  directly  usable  in  appropriate  optimization  algorithms  to  determine 
minimum  travel  time  path  and  minimum  distance  path  through  the  network  for  a 
variety  of  maneuver  unit  types  and  missions. 
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I.  INTRODUCTION 


A.  PURPOSE  AND  GOALS 

The  purpose  of  this  research  is  to  develop  algorithms  for  mapping  terrain 
characteristics  from  a  100  meter  grid  representation  to  arcs  and  nodes  of  a 
transportation  network.  The  algorithm  must  capture  elevation  and  trtfficability 
parameters  as  they  relate  to  both  the  arc  itself  and  to  the  off*arc  characteristics.  These 
network  parameters  must  be  directly  usable  m  appropriate  optimization  algorithms  to 
determine  minimum  travel  time  path  and  minimum  distance  path  through  the  network 
for  a  variety  of  maneuver  unit  types  and  missi  ons. 

After  reviewing  the  literature  and  discussing  the  problem  with  several  terrain 
model  expats,  it  was  determined  that  no  automated  procedure  for  accomplishing  this 
mapping  exists.  Therefore  the  procedures  and  associated  algorithms  presented  in  this 

thesis  represent  original  research  and  development.  j 

i 

R.  AIR  LAND  RESEARCH  MODEL  METHODOLOGIES 

The  purpose  of  the  Air  Land  Advanced  Research  Model  (  ALARM  )  is  to  j 

develop  a  systemic  corps  level  combat  model,  that  is,  one  that  operates  without  human 
intervention.  This  idea  is  not  new,  but  ALARM  takes  a  different  approach  than 
previous  efforts  to  simulate  the  military  decision  process.  j 

In  most  combat  models,  the  structures  for  planning  and  for  execution  are 
included  in  the  same  program.  In  A  LARK  the  planning  function  and  the  execution 
function  will  be  separate  programs.  The  Phoning  model  will  be  a  carefully  designed 
collection  of  programs  that  provide  orders  to  the  Execution  model  at  the  appropriate 
times. 

A  second  major  difference  between  the  a  LARM  design  and  existing  models  is  in 
the  method  used  to  make  decisions.  All  of  he  models  reviewed  when  developing  the 
ALARM  strategy  depended  either  on  threshold  parameters  or  on  a  series  of  “IF 
THEN*  decision  rules  to  make  decisions. 

Combat  models  using  a  threshold  parameter  strategy  usually  require  a  large 
number  of  threshold  parameters  or  groups  .  f  parameters.  In  such  models  a  threshold 
parameter  (  or  group )  is  needed  for  every  possible  decision  that  any  unit  represented  in 
the  combat  simulation  might  have  to  make.  This  leads  to  numerous  parameters  and 
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thresholds.  Most  parameters  of  this  type  have  no  reasonable  real  world  counterparts. 
In  other  words*  parameters  used  fbr  thresholding  are  frequently  contrived  and  based  on 
unexplainable  and  undefendable  measures.  This  makes  it  difficult  to  create  data  bases 
for  such  models,  and  difficult  to  place  any  confidence  in  the  analytical  results  from  the 
models. 

Models  that  use  a  series  of  decision  rules  are  lacking  in  that  they  require  a  very 
limited  view  of  the  decision  process.  A  decision  can  be  made  only  if  a  series  of  decision 
rules  to  support  the  decision  exist  within  the  software  that  is  developed  for  the  model, 
and  if  those  specific  rules  are  satisfied.  When  a  situation  is  encountered  for  which  there 
is  no  specific  decision  rule,  a  default  rule  must  be  implemented.  These  default  rules  are 
frequently  inappropriate.  The  more  complicated  the  decision  rules  become,  the  better 
the  model  usually  is,  but  the  more  difficult  it  becomes  to  maintain  the  model  and  the 
harder  it  is  to  analyze  results  from  the  model.  The  cause  and  etTect  relationships 
established  when  running  the  model  are  usually  closely  tied  to  the  decision  rules.  Thus 
the  output  or  results  from  such  a  simulation  are  built  into  the  model.  The  foregoing 
discussion  emphasizes  problems  with  the  methods  currently  accepted  and  used  to 
model  the  decision  process.  These  methods  are  recognized  because  they  have  their 
uses,  applications  and  foundations  in  the  human  decision  making  process.  Human 
beings  do,  on  some  level,  use  the  threshold  initiated  decision  rule  process  in  their 
decision  making  procedures.  These  concepts  are  usually  used  to  start  or  constrain  the 
decision  process  and  not  to  make  the  decision  itself. 

The  premise  of  ALARM  is  to  use  both  threshold  strategy  and  decision  rules,  but 
not  within  the  same  framework  as  other  models  to  date.  Thresholds  will  be  used  to 
determine  when  the  planning  or  decision  making  procedures  should  be  executed. 
Det^sion  rules  will  be  used  to  limit  alternatives,  "Hiey  will  not  be  used  explicitly  to 
make  combat  decisions.  ALARM  will  use  network  based  methodologies  to  itemize 
alternatives  so  that  a  decision  can  be  made.  It  is  believed  that  this  methodology  is 
closer  to  what  actually  happens  in  the  human  decision  making  process. 

;  ALARM  decision  processes,  called  Decision  Tasks,  will  use  three  unique 
methodologies  that  are  the  mainstay  of  the  ALARM  concept  and  make  it  different 
than  other  current  research.  These  methodologies  are  time  domain  networks,  cartesian 
space  networks,  and  the  Generalized  Value  System.  Each  methodology  is  briefly 
described  below. 
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1.  Tint  Domain  Networks 

Tune  domain  networks  are  designed  to  handle  the  planning  function  or 
activity  within  the  ALARM  system.  A  time  domain  network  consists  of  nodes  and  arcs 
(  or  Hnks )  connecting  the  nodes.  The  key  to  time  domain  networks  is  that  arcs  do  not 
represent  distance,  but  represent  the  passage  of  time  or  the  completion  of  a  sub¬ 
activity  which  leads  to  the  completion  of  the  entire  activity  represented  by  the  network, 
i  The  time  domain  network  is  used  to  develop  high  level  mission  requirements 

i 

for  ajl  subordinate  units.  It  is  similar  to  a  PERT  or  CPM  network  in  that  a  sub-activity 
cannjn  be  started  until  all  merging  activities  are  completed.  The  input  to  the  activity 
network  is  a  desired  activity  completion  time  and  an.  acceptable  friendly  force  attrition 
level,  and  may  include  a  required  level  of  attrition  to  enemy  forces. 

2.  Cartesian  Space  Networks 

The  Cartesian  space  network  is  a  series  of  networks  each  representing  a 
different  aspect  of  the  battlefield.  Each  network  represents  physical  connections 
betwjwn  points  on  the  battlefield.  ALARM  will  have  three  or  more  Cartesian  space 
networks.  The  three  networks  that  have  been  identified  are  terrain  and  transportation 
.  networks,  communication  networks,  and  logistics  resupply  networks.  This  thesis  deals 
explicitly  with  the  terrain  and  transportation  networks. 

3.  Generalized  Value  System 

The  generalized  value  System  (CVS)  is  a  procedure'  for  quantifying  the 
capabilities  and  importance  of  entities  on  the  battlefield  at  some  fhture  time  (t+x).  It 
does)  this  by  developing  algorithms  to  predict  fhture  entity  or  situation  states  at  time 

t 

(t+x)  based  on  the  situation  at  time  (t).  Thus  GVS  provides  the  framework  for 
foretjasting  future  states  of  entities  in  continuous  time. 

j  Various  combinations  of  exponential  functions  are  used  to  represent  these 
forecasts.  For  example,  the  power  and  value  of  non-combat  entities  such  as  bridges 
and  |oad  interdiction  points;  combat  support  units;  service  support  units;  and  combat 
unitsi  can  all  be  determined  from  a  continuous  function  with  respect  to  time  using  the 
various  exponential  algorithms.  These  values  can  then  be  used  to  determine  the  specific 
plac^  or  target  against  which  to  plan  an  air  strike  or  a  ground  attack  at  time  (t+x). 
Thusj  different  targets  can  be  selected  by  the  Decision  Task  depending  on  the  time  at 
which  strike  assets  are  available.  The  decision  will  be  based  not  on  the  current  value  of 
the  target,  but  on  its  predicted  future  value  in  relation  to  the  rest  of  the  Air  Land 
battle.  To  do  this,  all  targets  and  units  must  be  assigned  consistent  values  so  that 
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comparisons  can  be  made  across  the  diverse  set  of  targets  available  on  the  Air  Land 
battlefield.  The  common  bond  among  these  varying  targets  is  a  value  system  that 
measures  the  target's  value  in  terms  of  its  usefulness  to  combatants.  Thus  the  value  of 
a  bridge  is  not  measured  simply  by  its  width  or  length,  but  is  tied  to  the  number  of 
enemy  or  friendly  units  that  could  cross  the  bridge  to  maneuver  into  favorable  combat 
positions.  As  time  passes,  the  bridge's  value  could  increase  based  on  heavy  combat 
areas  and  then  could  decrease  when  all  relevant  units  have  successfully  crossed  the 
bridge. 

C.  SCOPE  OF  THE  THESIS 

In  attempting  to  achieve  the  goals  of  this  thesis,  various  terrain  modelling 
techniques  will  be  described  in  Chapter  II  as  a  motivation  for  the  network  algorithms 
developed.  The  algorithms  for  a  single  arc  will  form  the  *  isis  of  the  total  model  and 
will  be  presented  in  Chapter  III.  The  fUll  network  implementation  will  be  presented  in 
Chapter  IV  to  determine  the  minimum  time  and  minimum  distance  paths  on  an 
undirected  graph.  Chapter  V  will  discuss  the  results  of  several  runs  of  the  model  to 
illustrate  its  various  capabilities.  Finally,  Chapter  VI  will  provide  recommendations  for 
continued  enhancements  of  the  model. 
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II.  BACKGROUND 


A.  TERRAIN  MODELLING  TECHNIQUES 

Representing  combat  functions  as  a  syo*em  of  networks  is  a  departure  from  the 
current  practice  of  explicit  terrain  representation.  Current  modelling  techniques  require 
large  data  bases  to  represent  numerous  aspects  of  terrain  such  as  slope,  elevation, 
forestation,  a:id  cities,  fhe  network  approach  is  an  attempt  to  represent  the  terrain, 
primarily  through  use  of  the  transportation  system,  as  an  abstraction  of  its  actual  state. 
The  topology  characteristics  required  by  explicit  terrain  models  will  be  reflected  in  the 
arc  and  node  attributes. 

The  research  model  will  gain  an  economy  of  representation  of  the  topology 
because  it  has  the  flexibility  to  represent  only  those  terrain  features  applicable  to  a 
stated  level  of  resolution.  The  corps  rear  is  sparsely  populated  with  combat  and 
support  units  relative  to  the  overall  size  of  the  corps  sector.  Furthermore,  large 
portions  of  the  corps  sector  will  never  support  unit  movement  or  facilities  and  have  no 
reason  to  be  modelled.  A  network  representation  allows  one  to  select  only  those 
features  which,  because  of  their  nature,  give  the  controlling  force  a  marked  advantage 
over  the  other  force  (  e.g.  key  terrain  ).  Additionally,  those  roads  or  cross-country  arcs 
which  can  either  expedite  or  delay  movement  can  be  identified  and  modelled. 

In  a  corps  level  model,  direct  fire  combat  requires  observation,  detection,  and 
lines  of  sight  along  a  narrow  strip  of  terrain  relative  to  the  overall  size  of  the  sectors  of 
the  two  forces.  Existing  methods  of  terrain  representation  may  still  be  required  to 
model  the  details  of  such  combat.  The  majority  of  units  moving  in  a  corps  sector, 
especially  the  CS/CSS  assets,  do  so  without  receiving  direct  fire.  The  movement  of  the 
units  can  be  modelled  via  a  network  abstraction  of  the  transportation  system. 
Numerous  algorithms  exist  which  allow  one  to  exploit  the  features  and  structure  of 
such  a  mathematical  model. 

Many  of  the  terrain  modelling  techniques  in  current  use  restrict  the  model  to  one 
level  of  resolution  and  force  the  terrain  data  base  to  store  attributes  for  large  sections 
of  the  battlefield  which  are  not  used.  To  put  current  terrain  modelling  procedures  in 
perspective,  the  following  is  a  short  discussion  of  three  terrain  modelling  methods:  hex 
terrain  (  used  in  CORDIVEM  ),  digital  terrain  (DYNTACS)  and  network  model. 
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B.  HEX  TERRAIN 

The  Corps-Division  Evaluation  Model  (CORDIVEM)  is  a  corps  level  model  in 
which  the  terrain  is  represented  as  a  series  of  hexagons,  (hexes),  each  approximately  3.5 
kilometers  in  diameter.  The  terrain  underlying  each  hex  is  characterized  by  a  three  digit 
array  which  represents  urbanization,  forestation,  and  roughness.  Each  individual 
attribute  is  an  average  of  the  Defence  Mapping  Agency's  data  base  of  points  12.5 
meters  apart  which  underlie  each  CORDIVEM  hex.  The  attribute  based  on  the 
resulting  average  is  then  assigned  a  code  in  CORDIVEM  and  assumed  homogeneous 
throughout  the  hex  [Ref.  1:  p.  16].  Figure  2:1  illustrates  the  codes  for  each  hex 
attribute. 


Figure  2.1  CORDIVEM  Terrain  Features. 


The  geometry  of  the  hex  is  used  to  represent  unit  movement,  roads,  rivers  and 
obstacles.  Unit  movement  and  a  road  system,  if  one  exists,  are  modelled  from  the 
center  of  one  hex  to  the  center  of  an  adjacent  hex.  Movement  on  a  cross-country  route 
is  penalized  by  reducing  a  unit's  allowable  speed.  Rivers  are  modelled  along  the  edge  of 
a  hex,  and  coded  as  indicated  in  Figure  2.1.  An  implicit  bridge  is  assumed  at  the 
intersection  of  all  roads  and  rivers.  The  bridge  can  be  explicitly  modelled  if  it  may 
affect  movement  factors  or  possibly  be  destroyed.  An  obstacle  occurs  on  the  edge  of 


13 


the  hex  and  is  explicitly  inserted  into  or  removed  from  the  model.  All  obstacles  are 
created  before  the  execution  of  the  model;  the  dynamic  construction  and  reduction  of 
obstacles  are  not  modelled. 

C.  DIGITIZED  TERRAIN 

The  Dynamic  Tactical  Simulation  (DYNTACS)  model  is  a  two-sided,  dynamic 
Monte  Carlo  simulation  capable  of  modelling  the  combat  process  from  the  individual 
crew  to  battalion  engagements  which  utilizes  a  concept  referred  to  as  digitized  terrain 
[Ref,  2:  p.  9].  Macro  terrain  in  DYNTACS  is  usually  represented  as  100  meter  squares 
with  attributes  for  the  comers  being  supplied  by  data  available  from  the  Waterways 
Experiment  Station.  These  comer  attributes  provide  information  about  such  features  as 
elevation,  forestation,  and  location.  Each  square  is  then  divided  diagonally,  resulting  in 
a  series  of  equal  sized  triangles,  varying  in  slope  and  orientation. 

This  approach  assumes  that  the  terrain  modelled  by  each  triangle  is 
homogeneous  and  that  elevation  changes  spacings.  Sudden  changes  in  the  terrain  or  its 
surface  are  not  detected;  instead,  they  are  represented  with  geometric  overlays. 

D.  NETWORK  MODEL  FOR  MOBILITY  IN  ALARM 

The  Cartesian  space  network  compacts  the  terrain  database  still  further  by 
limiting  battlefield  movement  to  a  subset  of  battlefield  locations.  The  allowable 
locations  are  represented  as  a  network  overiayed  on  the  battlefield  map.  Nodes  of  the 
network  represent  physical  locations  on  the  map.  Arcs  of  the  network  represent 
possible  movement  paths  between  nodes  such  as  roads,  trails,  or  likely  cross  country 
routes. 

Each  node  and  arc  of  the  network  has  mobility  attributes  which  determine  the 
limiting  speed  of  combatants  moving  through  that  node  or  along  that  arc.  The 
attributes  which  are  stored  in  the  network  representation  can  be  either  terrain  features 
(  which  will  support  an  explicit  mobility  model  )  or  mobility  multipliers  (  for  a 
simplified  implicit  model ). 

Added  flexibility  can  be  gained  by  storing  several  sets  of  mobility  attributes  for 
some  arcs.  The  main  attribute  set  might  represent  movement  along  a  road,  while  a 
second  attribute  set  could  represent  cross  country  mobility  in  the  near  vicinity  of  the 
road  if  the  road  itself  cannot  be  used. 

The  main  problem  with  network  models  for  battlefield  movement  is  that 
i'ovement  is  restricted  to  the  network,  and  this  may  limit  tactical  flexibility. 
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Preparation  of  the  network  must  be  carefbily  coordinated  with  the  tactical  scenario  to 
ensure  sufficient  movement  opportunities.  This  limitation  would  seem  to  be  most 
severe  in  the  direct  fire  battle  zone  where  cross  country  maneuver  is  common.  Areas 
away  from  the  front  lines  can  more  readily  be  represented  by  a  network  model.  A 
major  advantage  of  the  network  model  is  that  it  allows  efficient  network  optimization 
algorithms  to  be  applied  to  the  problem  of  route  selection  for  moving  combatants. 

To  date,  data  for  network  attributes  have  been  derived  manually  using  map 
analysis.  This  technique  is  both  extremely  manpower  intensive  and  inaccurate.  The 
requirement  for  an  automated  process  to  generate  network  attributes  directly  from 
digitized  terrain  data  bases  serves  as  the  motivation  for  this  research.  The  algorithms  to 
determine  attributes  for  a  single  arc  are  described  in  the  next  chapter. 


III.  SINGLE  ARC  METHODOLOGY 


A.  OBJECTIVE  AND  PROBLEM  DEFINITION 

The  objective  of  this  research  is  to  develop  algorithms  for  mapping  terrain 
characteristics  from  a  100  meter  grid  representation  to  arcs  and  nodes  of  a 
transportation  network.  The  algorithm  must  capture  elevation  and  trafficability 
parameters  as  they  relate  to  both  the  arc  itself  and  to  the  off-arc  characteristics.  These 
network  parameters  must  be  directly  usable  in  appropriate  optimization  algorithms  to 
determine  minimum  travel  time  path  and  minimum  distance  path  through  the  network 
for  a  variety  of  maneuver  unit  types  and  missions. 

Methodology  for  the  single  arc  attributes  which  is  used  for  the  Cartesian  space 
network  is  presented  in  this  chapter  and  is  related  to  the  computer  program  for  the 
single  arc  attributes  in  Appendix  C. 

B.  MODEL  DATA  FILE 

The  program  used  four  disk  data  files,  one  for  100  meter  grid  square  data,  one 
for  the  node  characteristics,  one  for  the  arc  characteristics  which  represents  the 
networks  overlayed  on  the  gridded  terrain,  and  one  for  the  speed  data  related  to 
maneuvering  units. 

1.  100  Meter  Grid  Square  Data 

The  grid  square  data  represents  the  altitude  and  characteristics  of  each  point 
spaced  100  meters  apart  on  the  terrain.  The  source  of  the  data  is  the  Combined  Arms 
Center  at  Fort  Leavenworth,  which  is  the  same  data  base  originally  used  for 
CORDIVEM.  The  gridded  terrain  is  used  to  determine  the  attributes  of  networks 
overlayed  on  the  terrain.  A  typical  network  overlayed  on  the  digitized  terrain  is  shown 
in  Figure  3.1. 

Two  characteristics  are  used  to  describe  a  point  of  gridded  terrain:  three  digit 
altitude  and  an  integer  surface  feature  type.  The  altitude  is  the  height  from  sea  level  or 
simple  map  altitude.  The  characteristics  of  the  surface  feature  type  in  the  network  are 
used  for  determining  the  arc  width  for  maneuver  along  the  arc.  Table  1  shows  how 
data  for  the  100  meter  resolution  surface  feature  types  were  encoded. 
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Figure  3.1  Transportation  Network. 
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TABLE  1 

CHARACTERISTICS  OF  SURFACE  FEATURE  TYPES 


Integer  Value 

Surface  Feature 

0 

No  Slevation/No  feature  data 

1 

Forest 

2 

Urban 

3 

Harsh 

4 

Null (Elevation/No  feature  data) 

5 

Hater 

6 

Heath(Waste  land  with  shrubs) 

7 

Open 

2.  Node  Characteristics 

In  the  network,  nodes  are  used  to  represent  actual  terrain  features  on  a  map. 
Four  characteristics  are  used  to  describe  a  node:  an  integer  node  identification  number, 
three  digit  latitudinal  coordinate,  three  digit  longitudinal  coordinate,  and  node  type. 
The  coordinates  are  simple  map  grid  coordinates.  The  node  type  is  entered  as  a  single 
integer  value.  Eventually,  each  node  has  an  altitude  linearly  interpolated  from  100 
meter  grid  square  data.  Table  2  shows  the  integer  values  and  the  terrain  feature  it 
represents. 

The  primary  function  of  the  node  in  the  transportation  network  is  to  relate  the 
junction  of  arcs  in  the  data  representation  to  physical  locations  on  the  map.  In 
addition,  they  represent  possible  objectives  in  a  combat  mission. 

3.  Arc  Characteristics 

The  arcs  in  the  network  represent  feasible  routes  of  maneuver  from  one 
location  (node)  to  another.  For  a  route  to  be  considered  feasible,  it  must  provide  for 
the  maneuver  of  at  least  one  column  of  tracked  vehicles,  wheeled  vehicles,  or 
dismounted  troops.  Three  characteristics  are  used  to  represent  an  arc:  head  node 
identification,  tail  node  identification,  and  arc  type.  These  characteristics  provide  data 
for  the  algorithms  to  compute  appropriate  arc  traversal  times,  distances  and  flow  rates 
from  node  A  to  node  B. 
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TABLE  2 
TYPE  OF  NODE 


Integer  Value 

Terrain  Feature 

1 

City 

2 

Village 

3 

Road  Junction 

4 

Hill  Top 

5 

Other 

The  head  node  identification  is  the  identification  number  of  the  node  where 
the  arc  originates.  The  tail  node  identification  is  the  identification  number  of  the  node 
where  the  arc  terminates.  The  characteristics  of  the  arc  is  an  integer  cede  for  one  of 
seven  possible  arc  types  as  shown  in  Table  3. 


TABLE  3 
TYPE  OF  ARC 


Integer  Value 

Terrain  Feature 

1 

Autobahn 

2 

Concrete  or  Asphalt  Road 

3 

Dirt  Road 

4 

Railroad 

5 

Forest 

6 

Open  Country 

7 

Bridge  or  Tunnel 
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4.  Spied  Dau 

The  speed  for  unit  maneuver  is  a  (Unction  of  arc  type,  unit  type  and  unit 
(Urination.  The  possible  speed  represents  the  speed  on  the  specific  arc  (Ur  the  specific 
ur;t  for  each  unit  formation.  It  is  assumed  that  the  formation  width  of  a  single  column 
is  less  than  the  width  of  the  road  itself.  If  the  width  of  the  unit  formation  is  less  than 
the  width  of  the  road  itself,  then  arc  speed  is  on-road  speed.  Otherwise,  arc  speed  is 
off-road  speed,  because  on-road  speed  is  greater  than  off-road  speed  with  same  arc 
type  and  unit  type. 

C.  STEPS  IN  PROCESS 

To  translate  100  meter  grid  square  data  to  attributes  of  arcs,  the  network 
structure  is  used.  The  network  structure  can  exploit  the  mathematical  nature  of  the 
variables  of  the  decision  process  by  representing  these  variables  as  characteristics  of  the 
arcs  and  nodes  in  the  network  overlayed  on  the  gridded  terrain.  It  is  an  objective  in 
developing  this  network  to  adequately  describe  the  terrain  and  combat  effects  on 
maneuverability  necessary  for  the  implementation.  The  process  of  translation  and 
implementation  from  100  meter  grid  square  data  to  attributes  of  arcs  is  described  for 
single  arc  geometry  and  later  for  multiple  arc  geometry  and  flow  rates  in  the  network. 

1.  Single  Arc  Attributes 
a.  Gnmttry 

Consider  an  arc  which  links  node  A  to  node  B.  To  represent  the 
coordinates  of  a  node  in  three  dimensions,  the  Cartesian  coordinate  (x,  y,  z)  system  is 
used.  Each  node  has  an  altitude  interpolated  from  the  100  meter  grid  square  data. 
From  the  coordinates  of  node  A  and  node  B,  the  total  distance  between  node  A  and 
node  B,  and  the  angle  of  inclination  of  the  arc  is  computed.  The  angle  of  inclination  is 
determined  by  the  slope  between  the  x-axis  and  y-axis.  The  definition  of  angle  of 
inclination  of  a  line  that  crosses  the  x-axis  is  the  smallest  positive  angle  that  the  line 
makes  with  the  positively  directed  x-axis  [Ref.  3:  p.  50]. 

If  the  angle  of  inclination  is  between  45  degrees  and  135  degrees,  then  the 
x-axis  is  used  as  reference  axis  for  calculation  of  distance  between  each  point  along  the 
arc.  If  the  angle  of  inclination  is  greater  than  135  degrees  or  less  than  45  degrees,  then 
the  y-axis  is  used  as  the  reference  axis.  Using  this  angle,  coordinate?  of  the  cross 
points  between  a  point  on  the  arc  and  a  point  on  the  each  100  meter  grid  line  for  the 
reference  axis  is  computed.  After  determining  the  coordinates  of  the  points  along  the 
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arc,  the  distance  between  each  pair  of  points  is  calculated  step  by  step  until  the  cross 
point  reaches  the  tail  node  based  on  the  distance  formula  for  the  plane.  Figure  3.2 
shows  an  example  of  crossing  points  between  the  arc  (  which  is  Node  1  to  Node  2  ) 
and  the  x*axis  100  meter  grid  line  along  the  arc. 
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Figure  3.2  Crossing  Points. 
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(1)  Along  The  Arc.  To  determinfi  le  gradient  pattern  along  the  arc,  the 
slope  of  the  terrain  between  each  pair  of  points  along  the  arc  is  calculated.  Coordinates 
and  altitudes  of  each  point  along  the  arc  and  the  difference  of  altitudes  between  each 
sequential  pair  of  points  are  used  to  calculate  the  slopes  between  each  pair  of  points 
along  the  arc.  The  coordinates  of  each  point  along  the  arc  are  computed  by  using  the 
the  point-slope  equation  of  the  line.  As  a  measures  of  terrain  steepness  along  the  arc, 
three  factors  are  computed: 

•  Slope  Change  *  sum  of  the  absolute  value  of  slopes  between  each  pair  of 
adjacent  sub-arcs. 

•  Total  Slope  Change  -  sum  of  slope  change  along  the  arc. 

•  Average  Slope  Change  -  total  slope  change  divided  by  the  number  of  pair  of 
sub-arcs  along  the  arc. 

The  slope  of  the  surface  along  the  arc  is  related  to  the  feasible  routes 
of  maneuver  for  each  type  of  unit  from  one  node  to  another.  Some  arcs  provide  for 
the  maneuver  of  dismounted  troops,  but  tracked  and  wheeled  vehicle  movement  may 
be  infeasible.  Some  arcs  provide  for  the  movement  of  tracked  vehicles  to  climb  steep 
slopes  but  wheeled  vehicle  movement  is  not  possible.  The  maneuverability  on  a  given 
slope  depends  on  unit  type  (  tracked  vehicles,  wheeled  vehicles  or  dismounted  troops  ) 
and  soil  strength  which  is  measured  by  rating  cone  index  points:  the  higher  the  index 
numbers,  the  greater  the  soil  strength.  The  maximum  negotiable  slope  can  be  used  after 
obtaining  soil  strength  data  for  arcs  as  shown  in  Appendix  G.  [Ref.  4:  p.  5-19] 

Table  4  shows  an  example  of  general  information  and  slopes  along  the 
arc  from  Node  1  to  Node  2  as  shown  in  Figure  3.2.  XM  represents  the  slope  between 
starting  node  and  terminal  node  from  the  coordinates  of  the  x,  y  plane  on  the  arc. 
ANGLE  represents  the  angle  of  inclination  of  the  arc.  ACT.DIST  is  actual  distance 
between  two  nodes  along  the  surface,  but  DISTANCE  is  the  distance  on  the  plane  or 
on  the  map.  It  is  assumed  that  the  distance  of  the  arc  is  DISTANCE.  AVERAGE 
SLOPE  represents  the  difference  between  altitude  of  Node  l  and  Node  2  divided  by  the 
distance  of  the  arc. 

(2)  Perpendicular  to  Arc.  After  determining  the  slope  of  the  terrain  along 
the  arc,  the  slopes  of  the  terrain  along  lines  perpendicular  to  the  arc  are  computed.  The 
starting  point  is  the  cross  point  between  the  x-axis  grid  line  and  the  arc,  if  the  angle  of 
inclination  is  greater  than  45  degrees  and  less  than  135  degrees.  Otherwise,  the  starting 
point  is  the  cross  point  between  the  y-axis  grid  line  and  the  arc. 
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TABLE  4 

SLOPE  AND  DISTANCE 


NODS  2 

X  Y 

9000  9000 

NODE 

X 

8200 

1  XN  ANGLE 

Y 

9800  -1.00  138.00 

DISTANCE 

1131.37 

ACT.  DIST 

1131.95 

POINTS  ALONG 

ARC 

POSITION 

COORDINATE (  X  .  Y  ) 

ALTITUDE 

POINT 

1 

9000.0  9000.0 

367.0 

POINT 

2 

8900.0  9100.0 

367.0 

POINT 

3 

8800.0  9200.0 

367.0 

POINT 

4 

8700.0  9300.0 

372.0 

POINT 

5 

' 

8600.0  9400.0 

371.0 

POINT 

6 

8SC0.0  9S00. 0 

367.0 

POINT 

7 

8400.0  9600.0 

361.0 

POINT 

8 

3300.0  9700.0 

354.0 

POINT 

9 

8200.0  9800.0 

348.0 

POSITION 

DELTA  H  DISTANCE 

SLOPE 

ACTUAL  DIST 

POINT 

1  - 

2 

0.0  141.4 

0.000 

141.4 

POINT 

2  - 

3 

0.0  141.4 

0.000 

141.4 

POINT 

3  - 

4 

5.0  141.4 

0.035 

141.5 

POINT 

4  - 

5 

-1.0  141.4 

-0.007 

141.4 

POINT 

5  - 

6 

-4.0  141.4 

•0.028 

141.5 

POINT 

6  - 

7 

-6.0  141.4 

-0.042 

141.5 

POINT 

7  - 

8 

-7.0  141.4 

-0.049 

141.6 

POINT 

8  - 

9 

-6.0  141.4 

-0.042 

141.5 

POSITION 

SLOPE  CHANGE 

POINT 

1 

m 

2 

0.00 

POINT 

2 

m 

3 

0.04 

POINT 

3 

4 

0.04 

POINT 

4 

- 

S 

0.02 

POINT 

S 

m 

6 

0.01 

POINT 

8 

- 

7 

0.01 

POINT 

7 

8 

0.91 

TOTAL 

SLOPE 

CHANGE  t  0.1697 

AVERAGE  SLOPE  CHANGE  i  0.0212 

NODE 

2 1 

ALTITUDE  1  NODE 

li 

ALTITUDE  2 

367.0 

348 

.0 

AVERAGE  SLOPE 

<  -0.0168 

The  altitude  of  the  starting  point  is  calculated  by  linear  interpolation 
from  the  nearest  left  intersection  between  the  x-axis  and  y-axis  to  the  nearest  right 
intersection  between  the  x-axis  and  y-axis  line.  In  the  following  discussion,  it  is 
assumed  that  the  angle  is  135  degrees.  Table  5  and  Figure  3.3  are  used  to  illustrate  the 
following  discussion. 
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The  sequence  of  calculations  it  to  determine  the  slope  of  left  ofT-rotd 
and  right  off-road  along  the  perpendiculars.  On  the  left  side  the  first  distance  is  from 
the  starting  point  to  the  first  cross  point  between  the  perpendicular  line  and  y-axis  grid 
line.  The  second  distance  is  from  the  same  starting  point  as  the  first  one  to  the  second 
cross  point,  etc.  Toe  altitude  of  the  cross  point  between  the  perpendicular  line  to  the 
arc  and  the  y-axis  grid  line  is  computed  by  the  same  interpolation  procedure  as  for  the 
starting  point.  Using  the  first  distance  and  altitude  of  the  starting  point  and  the  first 
cross  point  between  the  perpendicular  line  and  y-axis  grid  line,  the  first  slope  is 
calculated.  The  process  is  continued  until  terminated  by  the  stopping  rule  described 
later.  After  computing  each  slope  on  the  perpendicular  line  for  the  initial  point  (  which 
is  the  head  node  ),  slopes  along  the  perpendicular  for  subsequent  points  along  the  arc 
nre  computed,  terminating  with  the  tail  node.  After  determining  the  slope  on  the  left 
side,  the  right  side  slopes  are  computed  using  the  same  procedures.  Before  calculating 
the  slope  on  the  perpendicular  line  to  the  arc,  coordinates  of  each  cross  point  between 
the  perpendicular  line  to  the  arc  and  the  y-axis  grid  line  are  determined.  The 
coordinates  of  each  point  on  the  perpendicular  line  to  the  arc  are  computed  using  the 
same  point-slope  equation  for  the  coordinates  along  the  a*c.  Also,  the  difference  in 
altitude  of  each  pair  of  points  is  computed. 

The  slopes  on  the  perpendicular  lines  are  related  to  the  feasible  arc 
width  for  movement  from  left  boundary  to  right  boundary.  The  arc  width  measured 
along  the  perpendicular  line  to  the  arc  is  considered  for  the  traflicable  path  for  each 
specific  type  of  unit.  Table  5  shows  an  example  of  slopes  used  to  determine  the 
tnfftcable  width  of  an  arc  for  the  movement  of  vehicle  units.  For  example,  POINT  4  - 
1  represents  the  fiist  (left)  cross  point  between  the  y-axis  100  meter  grid  line  and  the 
line  perpendicular  to  the  arc  at  the  fourth  point  from  the  starting  node.  POINT  4  -  5 
is  the  boundary  point  to  left  of  the  arc,  POINT  4  •  6  is  the  first  point  to  the  right,  and 
POINT  4  -  11  is  the  boundary  point  to  right  of  the  arc  as  shown  in  Figure  3.3. 
SLOPE  represents  DELTA  H  divided  by  DISTANCE. 

(3)  Arc  Width.  Two  types  of  stopping  rules  are  used  for  determining  the 
arc  width:  a  code  change  rule  and  a  slope  change  rule.  Briefly  stated,  an  arc  width 
boundary  is  established  along  a  line  perpendicular  to  an  arc  when  the  terrain  code 
changes  from  type  i  to  type  j  as  described  later.  Also,  a  boundary  may  be  established 
because  the  slope  change  along  the  perpendicular  exceeds  an  input  threshold  value. 
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POINT 

7-11 

POINT 

7-12 

00 
00 
00 
oo 
00 
00 
00 
00 
0 

8300 

8200 

8800 

8900 

9000 

9100 

9200 

9300 

8500 

8400 

8300 

200 

100 

000 

900 

800 

700 

800 

900 

000 

100 

8400 

8300 

8200 

8100 

8000 

7900 

7800 

7700 

8600 

8700 

8800 

8900 

8300 

8200 

3100 

8000 

7900 

7800 

7700 

7600 

7500 

8500 

8600 

8700 


8900 
9100 
9000 
9200 
9100 
9300 
9200 
9100 
9000 
8900 
8800 
9400 
9500 
9600 
9700 
9800 
9900 
9300 
9200 
9100 
9000 
8900 
8800 
8700 
8600 
9500 
9600 
9700 
9800 
9900 
9400 
9300 
9200 
9100 
9000 
8900 
8800 
8700 
9600 
9700 
9000 
9900 
95C0 
94C  0 
9300 
9200 
9100 
9000 
8900 
8800 
8700 
9700 
9800 
9900 


351.00 


348.00 


‘■I-I'K'I 


381..  00 


482.00 


'  VVR'I 


499.00 


351.00 


469,00 


381.00 

384.00 

400.00 

404.00 

407.00 


440.00 


351.00 


348.00 


707.11 


1131.37 


VJ 
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TABLE  5 

SLOPE  ON  PERPENDICULAR  LINE  TO  ARC  (CONT'D.) 


POINT 

8-  1 

8200 

9600 

361.00 

7.0 

141.42 

0.05 

POINT 

3-  2 

8100 

9500 

367.00 

13.0 

282.84 

0.05 

POINT 

3-  3 

8000 

9400 

374.00 

20.0 

424.26 

0.05 

POINT 

8-  4 

7900 

9300 

381.00 

27,0 

565.69 

0.05 

POINT 

8-  5 

7800 

9200 

388.00 

34.0 

707.11 

0.05 

POINT 

8-  6 

7700 

9100 

387.00 

33,0 

848.53 

0.04 

POINT 

8-  7 

7600 

9000 

390.00 

36.0 

989.95 

0.04 

POINT 

8-  3 

7500 

3900 

404.00 

50.0 

1131.37 

0.04 

POINT 

8-  9 

8400 

9800 

348.00 

-6 . 0 

141.42 

-0.04 

POINT 

8-10 

3500 

9900 

348.00 

-6.0 

282.84 

-0.02 

POINT 

9-  1 

8100 

9700 

354.00 

6.0 

141.42 

0.04 

POINT 

9-  2 

8000 

9600 

361.00 

13.0 

282.84 

0.05 

POINT 

9-  3 

7900 

9500 

367.00 

19.0 

424.26 

0.04 

POINT 

9-  4 

7800 

9400 

374.00 

26.0 

565.69 

0.05 

POINT 

9-  5 

7700 

9300 

381.00 

33.0 

707.11 

0.05 

POINT 

9-  6 

7600 

9200 

382.00 

34.0 

348.53 

0.04 

!  POINT 

9-  7 

(  8300 

9900 

348.00 

0.0 

141.42 

0.00 

I 
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i  The  stopping  rules  for  the  arc  width  are  a  function  of  unit  mission  ( 

unit. formation  )  and  unit  type.  When  the  unit  formation  is  a  single  column,  only  the 
on-road  portion  of  the  arc  is  used.  When  the  unit  formation  for  movement  is  multiple 
column,  both  on-road  and  off-road  attributes  may  be  used.  The  arc  width  for  a  given 
perpendicular  is  tne  distance  from  the  left  side  boundary  to  nght  side  boundary,  each 
established  by  one  of  the  two  rules. 

For  example,  in  case  of  a  vehicle  unit,  an  open  area  (  code  7  )  and  no 
feature  data  area  (  code  4  )  rr.ay  be  considered  feasible  for  routes  of  maneuver  for  off¬ 
road  movement,  witn  forest,  urban,  marsh,  water,  and  heath  areas  considered 
infeasible.  In  the  case  of  dismounted  troops,  however,  forest,  urban,  heath  and  open 
area  may  be  considered  feasible  for  movement  routes  off-road,  but  water,  and  marsh 
areas  arc  infeasible. 

Consider  the  matrices  shown  in  Tables  6  and  7.  At  the  first  cross  point 
between  the  perpendicular  line  and  y-axis  grid  line,  it  determines  the  code  from  the 
code  of  the  nearest  point  on  the  y-axis  grid  line.  If  the  value  contained  in  the  two 
dimensional  array  of  the  previous  code  and  current  code  equals  'no  boundary',  the 
process  continues.  If  the  value  of  previous  code  and  current  code  is  'boundary',  it  will 
stop.  For  example,  a  tracked  vehicle  can  go  from  open  area  to  open  area  but  can  not 
go  from  open  area  to  water  or  marsh  area.  The  dismounted  unit  can  go  from  forest  to 
forest,  open  or  heath  area,  but  will  not  go  from  forest  to  water  or  marsh  area.  That  is, 
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Figure  3.3  Points  on  the  Perpendicular. 

if  CODE(  i,  j  )  equals  1,  then  the  point  is  not  a  boundary,  otherwise,  the  point  is  a 
boundary.  Subscript  i  represents  the  previous  code  on  the  perpendicular  line  to  the  arc, 
and  subscript  j  represents  the  current  code  on  the  same  line  as  subscript  i. 

The  slope,  as  well  as  terrain  code  change,  could  be  a  barrier  to  certain 
types  of  units.  The  slope  of  the  perpendicular  line  to  the  arc  affects  the  movement  of 
units.  The  arc  boundary  due  to  slope  occurs  when  the  slope  is  less  than  a  specified 
negative  threshold  or  greater  than  a  positive  threshold.  The  threshold  of  slope  is  the 


maximum  negotiable  slope  for  each  unit  type.  The  threshold  for  vehicle  units  may  be 
less  than  that  for  dismounted  troops.  Eventually,  arc  soil  strength  data  will  be  used  to 
determine  these  thresholds  as  shown  in  Appendix  G. 
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Each  arc  width  on  the  perpendicular  line  is  determined  by 
consideration  of  the  code  change  rule  and  slope  change  rule  simultaneously.  The  arc 
width  is  the  distance  from  the  stopping  point  on  the  left  side  to  the  point  on  the  right 
side  of  the  arc.  The  line  connecting  each  stopping  point  on  the  left  (  right )  side  is  the 
arc's  left  (  right )  boundary  line  as  shown  in  Figure  3.4. 

2.  Determine  An:  Attributes 

Each  arc  width  on  the  perpendicular  line  to  the  arc  is  determined  by  the  code 
change  and  slope  change  rules  previously  discussed.  The  arc  widths  may  be  different 
for  each  perpendicular  line  along  the  arc.  However,  one  unique  arc  width  is  required 
from  a  head  node  to  a  tail  node,  since  it  is  difficult  to  change  the  unit  formation  during 
movement  on  an  arc.  Insertion  of  new  nodes  on  the  arc  can  be  considered  for  long 
arcs  when  half  of  the  arc  is  narrow  but  the  other  half  of  the  arc  is  wide. 

The  program  determines  the  arc  width  by  adding  the  left  minimum  distance, 
the  right  minimum  distance,  and  road  width  itself.  The  minimum  distance  on  the  left 
side  is  determined  by  the  shortest  perpendicular  distance  from  the  arc  to  the  left,  and 
similarly  for  the  right  minimum  distance.  When  the  left  and  right  minimum  distances 
°qual  zero,  the  arc  width  is  considered  to  be  the  road  width  itself  since  no  off-road 
movement  is  possible  for  that  unit  type.  The  minimum  distance  determination  of  arc 
width  is  more  realistic  than  consideration  of  some  average  of  the  distances.  The 
minimum  distance  has  an  direct  effect  on  the  flow  rate  of  an  arc.  Table  8  shows,  an 
example  of  the  minimum  distance  off  arc,  arc  width  and  those  coordinates  for  the  plot 
in  Figure  3.3. 

3.  Flow  Rate 

Using  the  arc  width  and  possible  speed  for  a  specified  unit  type  and  mission, 
the  rate  of  flow  is  computed.  Possible  speed  is  a  function  of  arc  type,  unit  type  and 
unit  mission  type.  The  arc  type  represents  several  conditions  of  the  maneuver  paths 
between  nodes  and  is  an  important  element  for  the  speed  of  maneuver.  There  are 
several  categories  of  arc  types  as  shown  in  Table  3  .  The  unit  types  currently 
considered  are  tracked  vehicles,  wheeled  vehicles  and  dismounted  troops  as  shown  in 
Table  18,  Appendix  F.  The  types  of  unit  formations  (  missions  )  currently  considered 
are  multiple  column  and  single  column  as  shown  in  Table  19,  Appendix  F. 

The  following  equation  is  used  to  determine  the  flow  rate  on  an  arc  and  is 
measured  in  terms  of  battalions  per  hour  for  each  unit  type  [Ref.  5:  p.  20]. 

FR  |  -  (SPbn  X  wp  /  (DWbn  X  DDbn)  (eqn  3.1) 
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ARC  WIDTH 


a 


TABLE  E 

MINIMUM  DISTANCE  OFF  ARC 


LEFT  MIN  DISTANCE 
RIGHT  MIN  DISTANCE 
ROAD  WIDTH 
ARC  WIDTH 


INITIAL  POINT  t 
END  POINT  < 


MINIMUM  DIJ 

s  141.4 
i  141.4 

i  6.0 

i  288.8 

POINT  1 

X  LEFT  Y  LEFT 


DISTANCE 


POINT  2 

X  RIGHT  Y  RIGHT 


8900.0 

8100.0 


8900.0 

9700.0 


9099.0 

8299.0 


9100.0 

9900.0 


The  equation  is  used  for  ail  arcs,  j,  to  compute  the  rate  of  flow,  FRj,  based  on 
the  speed  of  a  battalion  sized  unit,  SPbn,  the  width  of  the  arc,  W.,  the  doctrinal  width 
of  the  battalion,  DWbn,  and  its  doctrinal  depth,  DDbn.  The  doctrinal  width  and  depth 
of  the  battalion  depend  on  the  unit  type.  Since  the  doctrinal  width  and  depth  of  a 
battalion  of  dismounted  troops  is  different,  from  those  of  a  tank  battalion,  different 
values  were  assigned  for  these  parameters  as  a  function  of  unit  type.  Table  9  shows  an 
example  of  flow  rate  and  maximum  number  of  elements  on  the  arc  width  for  Arc  1 
(from  Node  1  to  Node  2)  as  shown  in  Figure  3.2. 


The  maximum  number  of  elements  which  can  move  in  multiple  column 
formation  on  the  arc  width  is  computed.  The  number  is  determined  from  the  arc  width 
and  the  distance  between  elements  ( i.e.,  the  maximum  number  of  elements  equal  the 
arc  width  divide  by  the  specified  distance  between  elements  for  each  formation  type  ). 

4.  Arc  Traversal  Time 

For  purposes  of  this  planning  model,  it  is  not  necessary  to  develop  extremely 
accurate  representations  of  attacking  force  movement.  General  estimates  of  expected 
movement  rates  based  on  unit  type  and  the  characteristics  of  the  terrain  being  crossed 
are  sufficient.  For  simplicity,  arc  traversal  time  is  the  arc  length  divided  by  arc  speed 
for  the  attacking  force  unit. 

Assume  the  formation  width  of  a  single  column  is  less  than  the  width  of  road 
itself.  If  the  width  of  the  unit  formation  is  less  than  the  width  of  the  road  itself,  then 
arc  speed  is  on-road  speed.  Otherwise,  arc  speed  is  off-road  speed.  On-road  speed  is 
assumed  to  be  larger  than  off-road  speed.  Arc  speeds,  in  terms  of  on-road  and  off-road 
speeds,  are  provided  in  Appendix  A. 


IV.  NETWORK  IMPLEMENTATION 


A.  INTRODUCTION 

Methodology  for  determination  of  minimum  time  and  minimum  distance  paths 
using  a  network  representation  of  terrain  is  presented  in  this  chapter.  This  chapter  is 
related  to  the  computer  program  for  the  Cartesian  space  network  in  Appendix  D.  The 
Cartesian  space  network  is  used  to  find  a  shortest  path  between  two  given  nodes  on  an 
undirected  graph  represented  by  a  list  of  arcs.  Input  data  for  this  network  program  is 

the  output  of  the  computer  program  of  single  arc  attributes  given  in  Appendix  C. 

/ 

B.  DATA  FILE 

The  network  representation  of  the  maneuver  area  is  placed  on  two  disk  files 
which  are  the  output  files  from  the  computer  program  for  single  arc  attributes:  one  for 
vehicle  units  and  one  for  dismounted  troops.  Eight  characteristics  are  used  to  represent 
an  arc:  arc  identification,  tail  node  identification,  head  node  identification,  traversal 
time,  flow  rate,  length,  width,  and  arc  speed.  These  characteristics  provide  sufficient 
data  for  the  algorithms  to  determine  the  minimum  time  or  distance  path  from  starting 
node  to  terminal  node.  The  time  for  the  arc  is  an  integer  representing  the  value  of 
traversal  time  for  the  movement  of  each  unit.  The  flow  rate  is  an  integer  value 
representing  flow  rate  for  arc,  j,  times  one  hundred  (  FR.  *  100  ).  The  integer  value  is 
used  for  computational  reasons  in  the  shortest  path  algorithm.  The  length  is  the 
dista^c  >etween  head  node  and  tail  node  in  kilometers.  The  width  of  the  arc  is  the 
width  v  oad  itself  plus  the  left  and  right  off-road  widths  in  kilometers  for  which  the 
terrain  *  support  movement.  The  arc  speed  represents  the  average  speed  for  each 
unit  along  the  arc. 

Two  k  ds  of  width  determination  rules  are  used  for  the  input  data  files  from  the 
single  arc  ^tributes  program.  The  width  determination  rule  is  used  in  the  following 
situation.  Recall  from  Chapter  III  that  either  a  terrain  code  or  slope  change  may 
terminate  the  search  for  the  arc  width  on  each  side  of  the  arc.  If  the  first  point  on  a 
perpendicular  to  an  arc  terminates  the  search,  the  width  determination  rule  is  invoked. 
For  example,  suppose  the  distance  to  the  first  point  is  100  meters.  If  rule  I  is  used,  the 
off-road  width  is  assigned  to  be  100  meters  for  that  perpendicular.  If  it  is  desired  to 
consider  off-road  width  to  be  zero,  then  Rule  2  is  used. 
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For  example,  if  the  unit  was  a  convoy  of  wheeled  vehicles,  then  Rule  2  would 
likely  be  selected.  For  tracked  vehicles  or  dismounted  troops,  it  may  be  appropriate  to 
select  Rule  1.  Table  10  shows  an  example  of  the  data  for  the  network  representation 
of  vehicle  unit  from  width  determination  rule  1. 


TABLE  10 

DATA  FOR  VEHICLE  UNIT  -  RULE  1 


ARC 

TAIL  HEAD 

TIME 

Unit* 

EUR 

~  ■  ~  -  -i 

DIST 

Km 

r"  -  •  — 

WIDTH 
Km 

SPEED 

Km/Hr 

1 

1 

2 

5 

240 

1.1 

0.289 

25 

2 

1 

8 

6 

94 

1.4 

0.113 

25 

3 

2 

3 

3 

106 

2.0 

0.128 

25 

4 

2 

9 

58 

1 

2.6 

0.010 

5 

5 

2 

10 

76 

20 

3.5 

0.120 

5 

6 

3 

4 

9 

97 

2.2 

0.117 

25 

7 

3 

10 

56 

20 

2.5 

0.123 

5 

3 

4 

5 

5 

144 

1.3 

0.174 

25 

9 

4 

10 

8 

102 

1.9 

0.123 

25 

10 

4 

11 

9 

76 

1.8 

0.115 

20 

11 

5 

6 

4 

184 

1.1 

0.222 

25 

12 

5 

11 

7 

89 

1.3 

0.134 

20 

13 

5 

12 

5 

102 

1.3 

0.123 

25 

14 

6 

12 

a 

98 

1.9 

0.113 

25 

15 

6 

7 

8 

103 

1.6 

•0.155 

20 

16 

7 

12 

12 

84 

2.2 

0,127 

20 

17 

7 

13 

10 

70 

1.8 

0.106 

20 

18 

8 

9 

9 

71 

1.7 

0.107 

20 

19 

8 

18 

40 

17 

1.8 

0.108 

5 

20 

9 

10 

59 

28 

2.7 

0.170 

5 

21 

9 

14 

35 

19 

1.6 

0.115 

5 

22 

9 

18 

48 

24 

2.2 

0.148 

5 

23 

9 

19 

52 

19 

2.4 

0.115 

5 

24 

10 

11 

5 

171 

1.3 

0.206 

25 

25 

10 

14 

3 

80 

1.6 

0.121 

20 

26 

10 

IS 

37 

20 

1.7 

0.123 

5 

27 

11 

12 

6 

70 

1.2 

0.105 

20 

28 

11 

15 

7 

93 

1.6 

0.U2 

25 

29 

11 

16 

a 

115 

1.9 

0.138 

25 

30 

12 

13 

8 

96 

1.6 

0.145 

20 

31 

12 

16 

6 

94 

1.5 

0.114 

25 

32 

13 

16 

43 

20 

2.0 

0.125 

5 

33 

13 

17 

40 

20 

1.8 

0.124 

5 

34 

14 

15 

7 

119 

1.6 

0.144 

25 

35 

14 

19 

47 

23 

2.1 

0.140 

5 

Unit*  :  Hour  x  100 

Unit**  :  (  Battalions  /  Hour  )  x  J00 


TABLE  10 

DATA  FOR  VEHICLE  UNIT  •  RULE  l  (CONT’D.) 


36 

15 

16 

8 

70 

1.5 

0.105 

20 

37 

15 

21 

5 

76 

1.0 

0.115 

20 

38 

16 

17 

59 

17 

2.7 

0.105 

5 

39 

16 

23 

9 

171 

2.1 

0.206 

25 

40 

16 

24 

42 

21 

1.9 

0.130 

5 

41 

17 

24 

5 

108 

1.3 

0.130 

25 

42 

18 

19 

57 

18 

2.6 

0.110 

5 

43 

18 

26 

58 

24 

2.6 

0.147 

5 

44 

18 

31 

95 

17 

4.3 

0.107 

5 

45 

19 

20 

5 

240 

1.1 

0.289 

25 

46 

19 

26 

27 

48 

1.3 

0.292 

5 

47 

20 

21 

9 

119 

2.2 

0.143 

25 

48 

20 

22 

8 

136 

2.0 

0.164 

25 

49 

20 

27 

7 

105 

1.7 

0.126 

25 

50 

20 

28 

45 

24 

2.1 

0.147 

5 

51 

21 

22 

21 

20 

1.0 

0.124 

5 

52 

22 

23 

6 

89 

1.2 

0.134 

20 

53 

22 

29 

7 

101 

1.3 

0.153 

20 

54 

23 

24 

7 

121 

1.7 

0.146 

25 

55 

23 

25 

8 

111 

1.8 

0.133 

25 

56 

23 

29 

3 

121 

0.8 

0.146 

25 

57 

24 

25 

8 

78 

1.6 

0.118 

20 

58 

25 

30 

6 

114 

1.1 

0.172 

20 

59 

26 

27 

33 

21 

1.5 

0.129 

5 

60 

26 

31 

65 

21 

3.0 

0.131 

5 

61 

27 

28 

42 

33 

1.9 

0.202 

5 

62 

27 

31 

14 

115 

2.7 

0.173 

20  • 

63 

27 

32 

7 

93 

1.7 

0.113 

25 

64 

27 

33 

58 

24 

2.6 

0.147 

5 

65 

28 

29 

43 

47 

2.0 

0.285 

5 

66 

28 

33 

39 

33 

1.8 

0.202 

5 

67 

28 

34 

47 

18 

2.2 

0.110 

5 

68 

29 

34 

10 

91 

2.3 

0.110 

25 

69 

30 

34 

9 

113 

2.2 

0.137 

25 

70 

31 

32 

11 

99 

2.6 

0.120 

25 

71 

32 

33 

9 

85 

1.7 

0.128 

20 

72 

33 

34 

49 

28 

2.2 

0.170 

5 

C.  DETERMINING  THE  PARAMETERS 

The  program  requires  several  parameters  for  determining  the  shortest  distance  or 
time  path.  Changeable  parameters  are  unit  type,  unit  formation,  starting  point  and 
terminal  point,  arc  travel  time  with  obstacles,  and  whether  minimum  time  or  distance  is 
to  be  calculated.  During  the  program  runs,  the  parameters  are  provided  by  user 
interaction  with  the  program. 


1.  (Jolt  tad  Path 

The  path,  can  b«  selected  in  two  ways:  minimum  time  or  minimum  distance. 
Minimum  time  determines  the  least  time  route  through  a  network  from  a  user  selected 
starting  node,  ISTART,  to  a  terminal  node,  LAST.  Minimum  distance  determines  the 
shortest  distance  path  between  two  selected  nodes. 

The  unit  type  is  vehicle  unit  or  dismounted  troops  based  on  a  battalion  sized 
unit.  Vehicle  unit  is  assumed  to  be  tracked  vehicles  for  the  examples  described  later. 
When  the  unit  is  determined,  the  program  reads  the  data  file  appropriate  for  the 
selected  unit. 

2.  Unit  Formation 

The  unit  formation  is  described  as  either  a  single  or  multiple  column 
formation.  For  dismounted  troops,  a  :v\  >  column  formation  is  defined  as  two 
columns,  with  multiple  column  formation  being  more  than  two  columns. 

The  width  of  the  unit  formation  is  the  number  of  columns  times  the  input 
spacing  between  columns.  The  depth  of  the  unit  formation  is  the  number  of  rows  times 
the  input  spacing  between  rows.  If  the  width  of  the  formation  is  larger  than  the  arc 
width,  then  the  program  assigns  a  Targe*  number  to  the  travel  time,  which  is  the  arc 
cost  in  the  shortest  path  algorithm.  Otherwise,  the  program  uses  the  original  value ‘for 
time  previously  determined  as  the  arc  cost. 

3.  Obstacles 

Obstacles  (  e.g.,  minefields  )  can  change  the  travel  time  and  the  shortest  path 
in  the  program.  The  program  requires  information  about  obstacles  on  the  arcs.  If 
there  is  an  obstacle  on  an  arc,  then  the  travel  time  is  appropriately  increased.  For  the 
purpose  of  examples  in  this  thesis,  the  original  travel  time  is  multiplied  by  four  when 
obstacles  are  present. 

Dv  SHORTEST  PATH  ALGORITHM 

The  determination  of  the  minimum  time  and  distance  path  through  the  sub¬ 
network  is  one  of  a  class  of  many  such  problems  often  referred  to  as  "shortest  path" 
network  problem.  The  general  shortest  path  problem  is  to  determine  the  least  cost 
route  through  a  network  starting  at  node  ISTART  and  ending  at  node  LAST.  The 
cost  of  a  route  is  some  function  of  the  characteristics  of  the  arcs  and  nodes  that  make 
up  the  route  from  node  ISTART  to  node  LAST.  In  the  case  of  this  minimum  time 
path  problem,  the  cost  of  traversing  an  arc  is  the  amount  of  time  it  takes  a  unit  to 
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traverse  the  arc.  The  cost  of  the  route  is  the  run  of  the  com  of  all  the  arcs  in  the  path 
from  the  supply  node  to  the  terminal  node. 

Before  selecting  a  method  of  finding  the  shortest  path,  the  network  must  be 
examined  to  insure  that  there  is,  in  fret,  a  solution  in  all  cases.  Because  no  arc  in  the 
network  will  have  a  negative  cost  associated  with  it.  the  addition  of  another  arc  to  any 
path  will  never  reduce  the  total  time  traveled.  Therefore,  negative  cycles  cannot  exist  in 
the  network.  The  only  arcs  which  have  a  zero  cost  are  those  leaving  the  supply  node 
and  entering  the  terminal  node.  Furthermore,  boundaries  are  selected  to  run  parallel  to 
avenues  of  approach  through  a  sector,  so  it  is  valid  to  assume  that  there  will  be  at  least 
one  set  of  connected  nodes  which  extends  laterally  through  the  sector.  Thus,  the' 
method  of  construction  of  the  subnetworks  insures  there  will  always  exist  at  least  one 
path  between  supply  and  terminal  nodes.  Therefore,  there  must  be  a  solution  to  the 
shortest  path  problem,  and  there  will  not  be  any  cycles  in  a  minimum  path(  no  node 
will  be  visited  twice).  If  the  assumption  that  boundaries  are  drawn  which  include  a 
connected  path  from  the  start  to  the  terminal  node  is  invalid,  a  solution  may  not  exist 
to  the  problem. 

Because  of  the  many  applications  for  this  class  of  problems,  there  are  many 
algorithms  available  for  its  solution.  Most  of  these  algorithms  consist  of  two 
procedures:  a  label  correcting  procedure  and  a  search  procedure. 

In  the  label  correcting  procedure  each  node  is  initially  assigned  an  infinite  cost 
with  the  exception  of  the  starting  node,  s  which  is  given  a  cost  of  zero.  Letting  c(u)  be 
the  cost  assigned  to  node  u,  d(u,\  /  be  the  cost  to  traverse  the  arc  from  node  u  to  node 
v,  and  pred(u)  be  the  node  previous  to  u  in  a  path,  then  the  following  rule  is  applied  to 

i 

change  the  costs  associated  with  node  v: 

> 

'  If  c(u)  +  d(u,v)  <  c(v)  (eqn  4.1) 

then  c(v)  ■  c(u)  +  d(u,v)  and  pred(v)  -  u. 

Though  inefficient  if  applied  indiscriminately,  if  this  rule  is  continuously  applied 
until  no  cases  of  Equation  4.1  being  true  can  be  found,  the  chain  cf  pred(v)  to  pred(u) 
"  s  will  describe  the  minimum  path  from  every  node  in  the  network  to  the  starting 
node,  s.  [Ref.  6:  p.  61} 

The  search  procedure  is  used  to  decide  in  which  order  the  nodes  are  to  be 
scanned  to  apply  the  labeling  rule.  This  step  is  where  most  algorithms  differ  and  also 
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where  the  efficiency  of  the  algorithm  it  achieved.  In  most  caaes  one  of  three  search 
techniques  are  used:  Dijkstra’s  algorithm,  depth  first  search,  and  breadth  first  search, 
but  many  hybrid  techniques  are  also  available.  Dijkstra’s  algorithm  gives  priority  of 
search  to  the  adjacent  node  which  is  the  shortest  distance  away.  Depth  first  search 
gives  priority  to  the  most  recent  node  searched.  Breadth  first  search  gives  priority  to 
the  oldest  node  searched.  The  Dijkstra's  algorithm  was  selected  for  use  in  the  program 
because  all  arc  lengths  are  at  least  greater  than  zero  and  are  not  the  same. 

As  the  name  implies,  the  shortest  path  algorithm  solves  arc  cost  in  terms  of  the 
minimum  travel  time  and  distance  for  an  attacker  across  an  arc.  For  simplicity,  the 
attacker  force  is  treated  as  a  rectangle  for  all  computations  and  travel  time  is  the  arc 
length  divided  by  average  arc  speed  of  the  attacking  force.  The  former  is  an  arc 
attribute,  while  the  latter  is  determined  from  unit  movement  capabilities.  The  algorithm 
is  described  below. 

Input:  Cartesian  Space  network  G(V,E),  nodes  V,  arcs  E  with  their 

characteristics  as  described  in  Section  B. 
output:  Shortest  path  distances  or  times(  d(v) )  from  starting  node,  s, 

to  terminal  node,  t. 

Step  1.  Initialization: 

a)  set  all  labels  to  a  very  large  value. 

d(v)  -  infinity 

b)  set  cost  of  starting  node  to  zero. 

d(s)  -  0 

c)  place  starting  node  in  the  Set  S. 

Set  S  -  {s} 

b)  remove  the  Set  S  from  Set  V,  the  set  of  nodes  of-G. 

Set  T  -  Set  V  -  {s} 

Step  2.  For  each  node  v  adjacent  to  the  node  u  (s.t.  u  6  S,  v  6  T). 

{  If  d(u)  +  l(u,v)  <  d(v) 
then  d(v)  *  d(u)  +  l(u,v)  } 

Step  3.  Determination  of  Set  S  and  Set  T. 

a)  determine  V  . 

min 

Vmin  "  ar8min  {d(v)} 

b)  save  the  value  of  cost  (d(v)}  in  the  Set  S. 

Set  S  -  Set  S  +  {V^} 

c)  remove  the  Set  S  from  the  set  of  nodes  of  G  (  Set  T ). 
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Set  T  -  Set  T  •  {«} 

d)  If  all  arcs  in  the  graph  have  not  been  analyzed,  go  to  Step  2. 

Step  5.  Return  Set  S  to  the  minimum  time  path  solution  routine.  Stop. 

l(u,v)  is  the  length  of  the  arc  between  nodu  u  and  adjacent  node  v. 

Set  S  is  the  set  of  nodes  such  that  d(v)  gives  the  true  shortest 
path  length  from  starting  node  s  to  some  node  v. 

Set  T  is  Set  V  •  Set  S 

Set  V  is  the  set  of  nodes  of  undirected  graph  C.  [Ref.  7:  p.  129] 

When  Set  T  is  empty,  each  node  in  the  network  will  have  a  cost  assigned  which 
represents  the  minimum  cost  for  ail  possible  paths  from  the  node  to  the  starting  node,s. 
The  minimum  path  from  any  node  to  the  start  node  can  then  be  found  by  tracing  the 
chain  of  predecessor  values  back  to  the  start  node.  When  the  algorithm  is  completed, 
the  minimum  time  or  distance  path  is  defined  by  the  chain  of  predecessor  values  from 
the  supply  r.ode.s,  to  the  terminal  node.  After  determining  the  node  number  along  the 
minimum  path,  the  arc  number  is  found  for  computing  total  travel  time  or  total 
distance  from  the  start  node  to  the  terminal  node.  The  total  travel  time  for  the  unit  is 
computed  as  follows. 

» 

T^-(IT«,)  +  (UL/SP)  (eqn  4.2) 

The  equation  is  used  for  total  travel  time  for  a  unit,  T11r1i.  based  on  the  total 
travel  time  for  an  element(  summation  of  j)  and  unit  length,  UL,  divided  by 
average  speed  along  the  path,  SP.  The  second  term  of  Equation  4.2  represents 
additional  travel  time  for  the  entire  unit  to  complete  the  path. 

The  following  equation  is  used  to  determine  average  speed  along  the  path. 

SP  -  £(d 1,  x  Vj)/D  (eqn 4.3) 

The  average  speed  SP  is  computed  based  on  the  summation  of  length  of  arC|,  d, 
multiplied  by  the  speed  on  arcj,  V{,  and  that  quantity  divided  by  the  total  path  length, 
D.  This  algorithm  finds  the  minimum  time  or  distance  path  through  the  sector  by 
using  the  unit  formation.  Any  type  of  arc  can  be  on  this  path  and  speeds  used  to 
determine  travel  time  arc  adjusted  for  the  type  of  arcs  along  the  minimum  time  path. 
Table  11  shows  an  example  of  results  for  the  vehicle  unit  from  the  network  program 
given  in  Appendix  D. 
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TABLE  11 

OUTPUT  FOR  VEHICLE  UNIT  -  RULE  1 


MIN  TIME  PATH 


VEHICLE  UNIT 

FORMATION  i  50  ROWS  1  COLUMNS 

ROW  SPACE  50  METERS 

NO  MINEFIELD  ARC 


SUM  OF  TRAVERSAL  TIME  s  0  HOUR  36  MINUTE 
TOTAL  TRAVERSAL  TIME  :  0  HOUR  42  MINUTE 

NODE  NUMBER  ALONG  MIN  TIME  PATH 
N  NODE 
1  1 

2  2 

3  3 

4  .4 

5  5 

6  12 


9  29 

10  34 

ARC  NUMBER  ALONG  MIN  TIME  PATH 
N  ARC 
1  1 

2  3 

3  6 

4  S 

5  13 

6  31 

7  39 

8  56 

9  68 

TOTAL  PATH  LENGTH  «  14.60  KM 


V.  ANALYSIS 


A.  INTRODUCTION 

This  chapter  describes  the  results  of  varying  parameter  values  of  the  single  arc 
attributes  program  and  Cartesian  Space  Network  program  used  in  determining 
minimum  distance  or  minimum  time  path.  The  results  in  this  chapter  are  final  outputs 
from  the  single  arc  attributes  program  through  the  Cartesian  Space  Network  program. 
A  10  X  10  Km  sector  of  European  terrain  is  used  for  the  analysis  to  demonstrate  the 
algorithms.  There  are  34  nodes  and  72  arcs  in  the  transportation  network.  The  start 
node  is  Node  1  and  the  terminal  node  is  Node  34.  The  unit  size  is  battalion  level 
consisting  of  50  tanks  in  a  tracked  vehicle  unit  and  1000  men  in  a  dismounted  troop 
unit.  Minimum  distance  path  is  first  discussed,  followed  by  a  comparative  analysis  of 
the  minimum  time  path  generated  under  various  conditions. 

B.  MINIMUM  DISTANCE  PATH 

This  algorithm  determines  the  path  that  provides  the  shortest  distance  through 
the  network.  The  arc  cost  is  defined  by  the  distance  between  head  node  and  tail  node. 
The  tactical  significance  of  this  algorithm  would  be  the  determination  of  the  path 
through  the  network  that  provides  the  best  route  to  reach  the  given  point  in  terms  of 
distance.  A  path  of  this  nature  might  be  used  by  raiding  or  reconnaissance  elements  of 
a  dismounted  troop  unit. 

The  minimum  distance  path  is  Path  D  through  the  network  shown  in  Figure  5.1 
The  distance  of  the  path  is  10.9  Km  from  the  start  node  (  Node  1  )  to  the  terminal 
node  (  Node  34  ).  Figure  5.1  shows  the  minimum  distance  path  and  several  minimum 
time  paths  described  in  the  next  section. 

C.  MINIMUM  TIME  PATH 

The  purpose  of  this  section  is  to  demonstrate  the  full  capabilities  of  the 
algorithms  developed  in  terms  of  minimum  time  paths.  The  following  parameters  are 
considered  in  the  analysis: 

•  Unit  type  -  vehicles  and  dismounted  troops. 

•  Formation  size  -  specified  by  number  of  rows  and  columns  in  the  formation. 

•  Width  determination  rule  -  specifies  the  method  used  to  determine  total  arc 
width  in  the  special  case  described  in  Chapter  IV. 

•  Obstacle  (  minefield  )  representation. 
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•  Bridge  representation. 

The  results  of  several  model  runs  for  various  parameter  values  are  given  in  Tables 
12  and  13.  The  total  travel  time  for  the  lead  element  and  the  unit,  the  path  length,  and 
the  minimum  time  path  from  Figure  5.1  are  presented  for  each  case.  The  impact  of  the 
width  determination  mle  selected  is  shown  in  Table  12  for  vehicle  units  and  in  Table  13 
for  dismounted  troop  units. 


TABLE  12 

MIN  TIME  PATH  -  VEHICLE  UNIT 


Case 

— 

Formation 
Row  x  Col 

A.  1 

50  X  1 

A. 2 

25  X  2 

A. 3 

10  X  5 

A. 4 

5  X  10 

B.l 

50  X  1 

B.2 

25  X  2 

B.3 

10  X  5 

B.4 

5  X  10 

C.l 

50  X  1 

C.2 

25  X  2 

C.3 

10  X  5 

C.4 

5  X  10 

D.l 

50  X  1 

D.2 

25  X  2 

D.3 

10  X  5 

D.4 

5  X  10 

E.l 

50  X  1 

E.2 

25  X  2 

E.3 

10  X  5 

E.4 

5  X  10 

F.l 

50  X  1 

F.2 

25  X  2 

F  .3 

10  X  5 

F.4 

5  X  10 

Bridge  Travel  Time  Length  Path 
y/n  element/unit  Km 


36  /  42 
36  /  39 
36  /  36 
no  feasible 


36  /  43 
36  /  39 
36  /  36 
no  feasible 


36  /  43 
36  /  40 
36  /  36 
no  feasible 


36  /  42 
no  feasible 
no  feasible 
no  feasible 


36  /  43 
no  feasible 
no  feasible 
no  feasible 


36  ./  43 
no  feasible 
no  feasible 
no  feasible 
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1.  Vehicle  Unit 

The  formation  of  vehicle  units  is  either  a  single  column  or  multiple  column 
formation.  The  width  of  the  formation  for  a  single  column  is  five  meters,  while  that  for 
the  multiple  column  formation  is  25  meters  for  each  column.  For  example,  the  width  of 
a  50  X  1  formation  is  five  meters,  of  a  25  X  2  formation  is  25  meters,  and  of  a  10  X  5 
formation  is  100  meters.  Note  from  Table  12  that  for  vehicle  units,  Path  A  and  Path  B 
are  selected  for  the  minimum  time  path  under  various  situations,  because  these  paths 
do  not  have  arcs  which  are  in  a  forest  or  marsh  area. 

a.  No  Minefield 

For  the  no  minefield  and  no  bridge  situation,  the  minimum  time  path  is 
Path  A  in  the  two  cases  which  use  width  determination  rules  1  and  2  (  see  Cases  A.l, 
A.2,  A.3,  and  D.l  ).  The  path  from  the  results  of  Rule  1  yields  an  optimal  route  for 
movement  until  the  formation  size  is  10  X  5  (  i.e.,  formation  width  is  100  meters  in 
which  case  no  feasible  route  exists(  see  Case  A.4  in  Table  12  )). 

For  Rule  2,  a  path  exists  only  for  the  single  column  formation  (  Case  D.l  ). 
For  wider  formations,  there  is  no  feasible  route  for  movement  (  see  Cases  D.2,  D.3, 
D.4  in  Table  12  ).  The  distance  of  the  minimum  time  path  is  14.6  Km  and  travel  time 
for  an  element  is  36  minutes,  as  shown  in  Table  12  .  Note  that  the  unit  travel  times 
vary  as  a  function  of  formation  size. 

b.  Minefield 

The  purpose  of  runs  labeled  Cases  B  and  E  in  . Table  12  is  to  demonstrate 
the  effect  of  insertion  of  obstacles  on  the  network.  For  these  runs,  a  minefield  is 
positioned  on  Arc  8  which  is  between  node  5  and  node  12  in  Figure  5.1.  Path  B 
bypasses  the  arc  which  has  a  minefield  for  both  Rule  1  and  2.  The  reason  is  that,  for 
these  runs,  arc  cost  is  four  times  larger  when  a  minefield  is  on  the  arc  as  compared  to 
the  no  minefield  case  ( i.e.,  the  travel  time  is  four  times  larger  for  a  minefield  arc ). 

Other  than  a  change  from  Path.  A  to  Path  B,  the  effects  of  the  width 
determination  rule  and  formation  size  are  similar  to  the  no  minefield  case. 

c.  Bridges 

Bridges  tend  to  constrain  movement  of  units  to  only  the  on-road  portion  of 
the  arc.  For  these  runs  (  shown  as  Cases  C  and  F  in  Table  12  ),  a  bridge  was  inserted 
on  arc  8  instead  of  a  minefield.  Even  though  Rule  1  (  Case  C  runs )  are  included,  only 
Rule  2  would  be  applied  in  the  real  situation.  The  minimum  time  path  bypasses  the 
arc  which  has  the  bridge  for  both  rules.  The  width  of  the  bridge  is  assumed  to  be  four 
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meters.  Because  the  width  of  the  bridge  is  less  than  the  width  of  the  formation  in  all 
cases,  the  minimum  time  path  algorithm  does  not  include  the  bridge  arc  on  the  optimal 
path.  The  results  for  this  case  are  the  same  as  for  the  minefield  case. 

2.  Dismounted  Troop  Units 

The  formations  for  dismounted  troop  units  are  considered  as  either  one  double 
column  (  six  meters  wide  )  or  multiple  column  (  five  meters  for  each  column  ).  For 
example,  a  500  X  2  formation  is  six  meters  wide,  while  a  50  X  20  formation  is  100 
meters  wide.  For  a  comparison  of  vehicle  versus  troop  units  with  the  same  formation 
widths,  consider  Cases  A.2,  G.2  and  D.2,  F.2  in  Tables  12  and  13.  Note  that  under 
Rule  1,  Path  C  is.  optimal  for  troop  units  (  Case  G.2  )  while  Path  B  is  optimal  for 
vehicle  units  (  Case  A.2  ).  For  the  wider  formations  (  Cases  G.4  through  G.6  )  the 
optimal  path  shifts  to  Path  B  under  Rule  1.  For  the  tighter  arc  width  restrictions  of 
Rule  2,  Path  B  is  optimal  in  ail  cases  for  dismounted  troop  units. 

Finally,  the  effect  of  long  units  for  travel  time  is  evident  from  Case  G.l.  Note 
that  an  additional  44  minutes  is  required  for  the  entire  unit  to  complete  Path  C. 

This  chapter  has  demonstrated  various  aspects  of  the  algorithms  developed  in 
this  thesis.  Except  for  those  areas  described  in  Chapter  VI,  the  model  is  ready  for  use 
in  a  full  production  mode  for  any  terrain  areas  which  have  digitized  terrain  data  and 
transportation  networks  in  the  form  described  in  Chapter  III. 


TABLE  13 

MIN  TIME  PATH  -  DISMOUNTED  TROOP  UNITS 


Case 

Formation 
Row  x  Col 

Mine 

y/n 

Bridge 

y/n 

Travel  Time 
element/unit 

Length 

Km 

Path 

G.l 

500 

X 

2 

1 

no 

no 

4:22 

/ 

5:06 

11.9 

C 

G.2 

200 

X 

5 

1 

no 

no 

4:22 

/ 

4:39 

11.9 

C 

G.3 

50 

X 

20 

1 

no 

no 

4:22 

/ 

4:27 

11.9 

C 

G.4 

20 

X 

so 

1 

no 

no 

4:40 

/ 

4:43 

14.2 

B 

G.5 

10 

X 

100 

1 

no 

no 

4:40 

/ 

4:42 

14.2 

B 

G.6 

7 

X 

140 

1 

no 

no 

4:40 

/ 

4:41 

14.2 

B 

F.l 

500 

X 

2 

2 

no 

no 

4:40 

/ 

5:31 

14.2 

B 

F.2 

200 

X 

5 

2 

no 

no 

4:40 

/ 

5:00 

14.2 

B 

F.3 

50 

X 

20 

2 

no 

no 

4:40 

/ 

4:45 

14.2 

B 

F.4 

20 

X 

50 

2 

no 

no 

4:40 

/ 

4:43 

14.2 

B 

F.5 

10 

X 

100 

2 

no 

no 

4:40 

/ 

4:42 

14.2 

B 

F.6 

7 

X 

140 

2 

no 

no 

4:40 

/ 

4:41 

14.2 

B 
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VI.  SUMMARY  AND  FUTURE  RESEARCH  AREAS 


It  is  essential  that  any  combat  simulations  of  the  Airland  battle  provide  a 
realistic  representation  of  terrain  and  environment.  Previous  approaches  have  typically 
utilized  either  hex  or  grid  squares  within  the  model  to  describe  terrain  conditions.  In 
some  cases  network  overlays  were  used  for  convoy  movements,  but  no  models  exist 
which  capture  essential  terrain  characteristics  without  the  internal  use  of  some  grid 
pattern. 

In  order  to  take  advantage  of  the  various  algorithms  available  using  networks 
without  the  need  for  gridded  terrain  representation,  it  is  necessary  to  describe  essential 
terrain  parameters  as  arc  and  node  attributes.  A  family  of  algorithms  were  developed 
and  demonstrated  in  this  thesis  which  accomplishes  this  objective.  These  algorithms 
can  be  implemented  for  any  terrain  area  for  which  appropriate  data  exists.  The 
resulting  network  is  then  usable  in  a  model  such  as  ALARM  for  both  the  planning  and 
execution  phases  of  combat  representation. 

Areas  of  future  research  include  the  development  of  a  maximum  flow  path 
algorithm.  This  area  requires  a  significant  effort  in  determining  appropriate  arc  costs 
and  development  of  maximum  flow  methodology  for  an  undirected  network.  In 
addition,  several  refinements  are  possible  to  represent  other  aspects  of  terrain 
characteristics  not  explicitly  contained  in  the  current  model.  Finally,  additional  test 
runs  of  the  model  with  different  terrain  areas,  unit  types,  combat  formations, and 
boundary  rules  need  to  be  made. 
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Arc 


Node 
H  T 


Road  Off-Road  Off-Road  Total  Speed 
Width  Left  Right  Arc  Width  On  off 


1 

2 

6 

.  141 

141 

288 

40 

25 

1 

8 

6 

100 

7 

113 

4C 

25 

2 

3 

6 

109 

12 

127 

40 

25 

2 

9 

2 

4 

4 

10 

5 

5 

2 

10 

2 

108 

10 

120 

5 

5 

3 

4 

6 

111 

0 

117 

40 

25 

3 

10 

2 

1 

120 

123 

5 

5 

4 

S 

6 

112 

56 

174 

40 

25 

4 

10 

6 

111 

7 

123 

40 

25 

4 

11 

5 

6 

104 

115 

30 

20 

5 

6 

6 

133 

83 

222 

40 

25 

5 

11 

4 

119 

11 

134 

30 

20 

5 

12 

6 

9 

108 

123 

40 

25 

6 

12 

6 

101 

11 

118 

40 

25 

7 

6 

4 

129 

22 

155 

30 

20 

7 

12 

4 

1-3 

0 

127 

30 

20 

7 

13 

4 

101 

0 

106 

30 

20 

8 

9 

5 

102 

0 

107 

30 

20 

8 

18 

2 

6 

100 

108 

5 

5 

9 

10 

2 

112 

56 

170 

5 

5 

9 

14 

2 

101 

13 

115  ' 

5 

5 

9 

18 

2 

137 

9 

148 

5 

5 

9 

19 

2 

5 

108 

115 

5 

5 

10 

11 

6 

100 

100 

206 

40 

25 

10 

14 

5 

115 

1 

121 

30 

20 

10 

15 

2 

8 

113 

123 

5 

5 

11 

12 

5 

100 

0 

105 

30 

20 

11 

15 

6 

105 

0 

112 

40 

25 

11 

16 

6 

15 

117 

138 

40 

25 

12 

13 

4 

141 

0 

145 

30 

20 
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TABLE  14 

VEHICLE  UNIT  -  RULE  1  (CONTD.) 


31 

12 

16 

6 

101 

7 

114 

40 

25 

32 

13 

16 

2 

116 

7 

125 

5 

5 

33 

13 

17 

2 

7 

114 

124 

5 

5 

34 

14 

15 

1  * 

100 

38 

144 

40 

25 

35 

14 

19 

2 

122 

16 

140 

5 

5 

36 

15 

16 

5 

100 

0 

105 

30 

20 

37 

15 

21 

4 

100 

10 

115 

30 

20 

38 

16 

17 

2 

103 

0 

105 

5 

5 

39 

16 

23 

6 

100 

100 

206 

40 

25 

40 

16 

24 

2 

128 

0 

130 

5 

5 

41 

17 

24 

6 

114 

10 

130 

40 

25 

42 

18 

19 

•> 

2 

108 

0 

110 

5 

5 

43 

18 

26 

138 

7 

147 

5 

5 

44 

18 

31 

2 

e 

100 

107 

5 

5 

45 

19 

20 

6 

141 

141 

289 

40 

25 

46 

19 

26 

1  2 

228 

62 

292 

5 

5 

47' 

20 

21 

6 

130 

8 

143 

40 

25 

48 

20 

22 

6  - 

109 

49 

164 

40 

25 

49 

20 

27 

6 

107 

13 

126 

40 

25 

50 

20 

28 

2 

16 

128 

147 

5 

5 

51 

21 

22 

2 

12 

109 

124 

5 

5 

52 

22 

23 

5 

103 

26 

134 

30 

20 

53 

22 

29 

5 

14 

134 

153 

30 

20 

54 

23 

24 

6 

117 

23 

146 

40 

25 

55 

23 

25 

6 

115 

13 

133 

40 

25 

56 

23 

29 

6 

109 

31 

146 

40 

25 

57 

24 

25 

5 

13 

101 

118 

30 

20 

58 

25 

30 

4 

56 

112 

172 

30 

20 

59 

26 

27 

2 

9 

11? 

129 

5 

5 

60 

26 

31 

2 

119 

10 

131 

5 

5 

61 

27 

28 

2 

100 

100 

202 

5 

5 

62 

27 

31 

5 

112 

56 

173 

30 

20 

63 

27 

32 

6 

6 

101 

113 

40 

25 

64 

27 

33 

2 

138 

7 

147 

5 

5 

65 

28 

29 

2 

141 

141 

285 

5 

5 

66 

28 

33 

2 

100 

100  1 

202 

5 

5 

67 

28 

34 

2 

108 

0 

110 

5 

5 

68 

29 

34 

6 

0 

104 

110 

40 

25 

69 

30 

34 

6 

118 

12 

137 

40 

25 

70 

31 

32 

6 

114 

0 

120 

40 

25 

71 

32 

33 

4 

100 

24 

128 

30 

20 

72 

33 

34 

2 

112 

56 

170 

j 

5 

5 

TABLE  15 

VEH^LE  UNIT  •  RULE  2 


VM 


M 


m 

r 


m 


Arc 

Nod* 

H  T 

Ro*d 

Width 

Off-Ro*d 

L*ft 

Off**  Rood 
Right 

Total 

Arc  Width 

Sp**d 
On  Off 

1 

1 

2 

6 

0 

0 

6 

40 

25 

2 

1 

8 

6 

0 

0 

6 

40 

25 

3 

2 

3 

6 

0 

0 

6 

40 

25 

4 

2 

9 

2 

0 

0 

2 

5 

5 

5 

2 

10 

2 

0 

0 

2 

5 

5 

6 

3 

4 

6 

0 

0 

6 

40 

25 

7 

3 

10 

2 

0 

0 

2 

5 

5 

8 

4 

5 

6 

0 

0 

6 

40 

25 

9 

4 

10 

6 

0 

0 

6 

40 

25 

10 

4 

11 

5 

0 

0 

5 

30 

20 

11 

5 

6 

6 

0 

0 

6 

40 

25 

12 

5 

11 

4 

0 

0 

4 

30 

20 

13 

5 

12 

6 

0 

0 

6 

40 

25 

14 

6 

12 

6 

0 

0 

6 

40 

25 

15 

7 

6 

4 

0 

0 

4 

30 

20 

16 

7 

12 

4 

0 

0 

4 

30 

20 

17 

7 

13 

4 

0 

0 

4 

30 

20 

18 

8 

9 

5 

0 

0 

5 

30 

20 

19 

8 

18 

2 

0 

0 

2 

5 

5 

20 

9 

10 

2 

0 

0 

2 

5 

5 

21 

9 

14 

2 

0 

0 

2 

5 

5 

22 

9 

18 

•j 

0 

0 

2 

5 

5 

23 

9 

19 

2 

0 

0 

2 

5 

5 

24 

10 

11 

6 

0 

0 

S 

40 

25 

25 

10 

14 

5 

0 

0 

5 

30 

20 

26 

10 

15 

2 

0 

0 

2 

5 

5 

27 

11 

12 

5 

0 

0 

5 

30 

20 

28 

11 

15 

6 

0 

0 

6 

40 

25 

29 

11 

16 

6 

0 

0 

6 

40 

25 

30 

12 

13 

4 

0 

0 

4 

30 

20 

31 

12 

16 

6 

0 

0 

6 

40 

25 

32 

13 

16 

2 

0 

0 

2 

5 

5 

33 

13 

17 

2 

0 

0 

2 

5 

5 

34 

14 

15 

6 

0 

0 

6 

40 

25 

35 

14 

19 

2 

0 

0 

2 

5 

5 

36 

15 

16 

5 

0 

0 

5 

30 

20 

37 

15 

21 

4 

0 

0 

4 

30 

20 

38 

16 

17 

2 

0 

0 

2 

5 

5 

39 

16 

23 

6 

0 

0 

6 

40 

25 

40 

16 

24 

2 

0 

0 

2 

5 

5 

41 

17 

24 

6 

0 

0 

6 

40 

25 

42 

13 

19 

2 

0 

0 

2 

5 

5 

43 

13 

26 

2 

0 

0 

2 

5 

■  5 

44 

13 

31 

2 

0 

0 

2 

5 

5 

45 

19 

20 

6 

0 

0 

6 

40 

25 

46 

19 

26 

2 

0 

0 

2 

5 

5 

47 

20 

21 

6 

0 

0 

6 

40 

25 

48 

20 

22 

6 

0 

0 

6 

40 

25 

49 

20 

27 

6 

0 

0 

6 

40 

25 

50 

20 

28 

2 

0 

0 

2 

5 

5 

TABLE  16 

DISMOUNTED  TROOP  UNIT  *  RULE  1 


Nod* 

H  T 


2  10 
3  4 

3  10 

4  5 

4  10 

4  11 

5  6 

5  11 

5  12 

6  12 

7  6 

7  12 

7  13 

8  9 
8  18 

9  10 

9  14 

9  18 

9  19 

10  11 
10  14 

10  15 

11  12 

11  15 
11  16 

12  13 
12  16 

13  16 

13  17 

14  15 

14  19 

15  16 

15  21 

16  17 

16  23 

16  24 

17  24 

18  19 
18  26 

18  31 

19  20 

19  26 

20  21 

20  22 
20  27 

20  28 


Road 

Width 

Off-Ro*d 

Lift 

6 

990 

141 

6 

100 

7 

6 

1204 

218 

2 

4 

105 

2 

1187 

975 

6 

1118 

203 

2 

1202 

160 

6 

1230 

670 

6 

1216 

59 

5 

105 

1042 

6 

1462 

664 

4 

1304 

1078 

6 

1092 

1083 

6 

1114 

916 

4 

1185 

538 

4 

1357 

1110 

4 

1115 

912 

5 

113 

0 

2 

6 

77 

2 

1230 

1062 

2 

1109 

920 

2 

137 

34 

2 

5 

108 

6 

1100 

1000 

5 

1267 

1037 

2 

1141 

1133 

5 

1104 

903 

6 

1160 

949 

6 

1194 

1179 

4 

j.  J.60 

949 

6 

1110 

915 

2 

130 

75 

2 

36 

1147 

6 

1102 

908 

2 

1165 

959 

5 

1102 

902 

4 

1105 

915 

2 

1129 

47 

6 

1100 

1000 

2 

1409 

1153 

6 

1253 

1036 

2 

108  • 

0 

2 

181 

73 

2 

21 

100 

6 

141 

1414 

2 

1253 

1036 

6 

1425 

1174 

6 

1204 

997 

6 

1175 

975 

2 

1297 

1288 

sp««d 

On  Off 


1137 

113 

1428 

111 

2164 

1327 

1364 

1906 

1280 

1152 

2132 

2386 

2181 

2036 

1727 

2471 

2032 

118 

85 

2294 
2030 

173 

115 

2106 

2309 

2276 

2012 

2114 

2379 

2112 

2030 

207 

1185 

2016 

2126 

2009 

2024 

1178 

2106 

2563 

2295 
110 
256 
123 

1562 

2291 

2605 

2207 

2155 

2587 


TABLE  16 

DISMOUNTED  TROOP  UNIT  -  RULE  I  (CONTD.) 


SI 

21 

22 

2 

1107 

1094 

2203 

52 

22 

23 

5 

1134 

953 

2092 

53 

22 

29 

5 

1359 

1345 

2709 

S4 

23 

24 

6 

1263 

1050 

2338 

55 

23 

25 

6 

1190 

960 

2175 

56 

23 

29 

6 

1197 

995 

2198 

S7 

24 

25 

5 

1020 

1008 

2033 

56 

25 

30 

4 

1174 

1118 

2296 

59 

26 

27 

2 

1183 

1174 

2359 

60 

26 

31 

2 

1306 

474 

1782 

61 

27 

28 

2 

1100 

1000 

2102 

62 

27 

31 

5 

1230 

112 

1347 

63 

27 

32 

6 

6 

101 

113 

64 

27 

33 

2 

138 

1269 

1408 

65 

28 

29 

2 

1556 

1414 

2972 

66 

28 

33 

2 

1100 

1000 

2102 

67 

28 

34 

2 

1185 

0 

1187 

68 

29 

34 

6 

94 

1036 

1136 

69 

30 

34 

6 

1301 

1071 

2378 

70 

31 

32 

6 

204 

1504 

1714 

71 

32 

33 

4 

100 

1903 

2007 

.72 

33 

34 

2 

112 

1062 

1176 
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TABLE  17 

DISMOUNTED  TROOP  UNIT  -  RULE  2 


Node 

H  T 

Road 

Width 

Off-Road 

Left 

Off-Road 

Right 

Total 

Arc  Width 

Speed 

On  Off 

1 

2 

6 

990 

141 

1137 

4 

3 

1 

8 

6 

0 

0 

6 

4 

3 

2 

3 

6 

1204 

218 

1428 

4. 

3 

2 

9 

2 

0 

0 

2 

4- 

2 

2 

10 

2 

1187 

975 

2164' 

4 

2 

3 

4 

6 

1118 

203 

1327' 

4 

3 

3 

10 

2 

1202 

160 

1364 

4 

2 

4 

5 

6 

1230 

670 

x906 

4 

3 

4 

10 

6 

1216 

59 

1280 

4 

3 

4 

11 

5 

0 

1042 

1047 

4 

3 

5 

6 

6 

1462 

664 

2132 

4 

* 

w 

5 

11 

4 

1304 

1078 

2386 

4- 

3 

5 

12 

6 

1092 

1083 

2181 

4 

3 

6 

12 

6 

1114 

916 

2036 

4 

3 

7 

6 

4 

1185 

538 

1727 

4 

3 

7 

12 

4 

1357 

1110 

2471 

4 

3 

7 

13 

4 

1115 

912 

2032 

4 

3 

8 

9 

5 

0 

0 

5 

4 

3 

8 

18 

2 

0 

0 

2 

4 

2 

9 

10 

2 

1230 

1062 

2294 

4 

2 

9 

14 

2 

1109 

920 

2030 

4 

2 

9 

18 

2 

0 

0 

2 

4 

2 

9 

19 

2 

0 

0 

2 

4 

2 

10 

11 

6 

1100 

1000 

2106 

4 

3 

10 

14 

5 

1267 

1037 

2309 

4 

3 

10 

15 

2 

1141 

1133 

2276 

4 

2 

11 

12 

5 

1104 

903 

2012 

4 

3 

11 

♦  m 

x.J 

6 

1160 

949 

2114 

4 

3 

11 

16 

6 

1194 

1179 

2379 

4 

3 

12 

13 

4 

1160 

949 

2112 

4 

3 

12 

16 

6 

1110 

915 

2030 

4 

3 

13 

16 

2 

0 

75 

77 

4 

2 

13 

17 

2 

0 

1147 

1149  • 

4 

2 

14 

15 

6 

1102 

908 

2016 

4 

3 

14 

19 

2  ■ 

1165 

959 

2126 

4 

2 

15 

16 

5 

1102 

902 

2009 

4 

3 

15 

21 

4 

1105 

915 

2024 

4 

3 

16 

17 

2 

1129 

47 

1178 

4 

2 

16 

23 

6 

1100 

1000 

2106 

4 

3 

16 

24 

2 

1409 

1153 

2563 

4 

2 

17 

24 

6 

1253 

1036 

2295 

4 

3 

18 

19 

2 

0 

0 

2 

4 

2 

18 

26 

■2 

0 

0 

2 

4 

2 

18 

31 

2 

0 

0 

2 

4 

2 

19 

20 

6 

0 

1414 

1420 

4- 

3 

19 

26 

2 

1253 

1036 

2291 

4 

2 

20 

21 

6 

1425 

1174 

2605 

4 

3 

20 

22 

6 

1204 

997 

2207 

4 

3 

20 

27 

6 

1175 

975 

2155 

4 

3 

20 

28 

2 

1297 

1288 

2587 

4 

2 

TABLE  17 

DISMOUNTED  TROOP  UNIT  -  RULE  2  (CONT'D.) 


51 

21 

22 

2 

1107 

1094 

2203 

4 

2 

52 

22 

23 

5 

1134 

953 

2092 

4 

3 

53 

22 

29 

5 

1359 

1345 

2709 

4 

3 

54 

23 

24 

6 

1283 

1050 

2338 

4 

3 

55 

23 

25 

6 

1190 

980 

2175 

4 

3 

56 

23 

29 

6 

1197 

995 

2198 

4 

3 

57 

24 

25 

5 

1020 

1008 

2033 

4 

3 

58 

25 

30 

4 

1174 

1118 

2296 

4 

3 

59 

26 

27 

2 

1183 

1174 

2359 

4 

2 

60 

26 

31 

2 

1306 

475 

1783 

4 

2 

61 

27 

28 

2 

1100 

1000 

2102 

4 

2 

62 

27 

31 

5 

1230 

112 

1347 

4 

3 

63 

27 

32 

6 

0 

0 

6 

4 

3 

64 

27 

33 

2 

0 

1269 

1271 

4 

2 

65 

28 

29 

2 

1556 

1414 

2972 

4 

2 

66 

28 

33 

2 

1100 

1000 

2102 

4 

2 

67 

28 

34 

2 

1185 

0 

1187 

4 

2 

68 

29 

34 

6 

0 

1036 

1042 

4 

3 

69 

30 

34 

6 

1301 

1071 

2378 

4 

3 

70 

31 

32 

6 

204 

1504 

1714 

4 

3 

71 

32 

33 

4 

0 

1903 

1907 

4 

3 

72 

33 

34 

2 

0 

1062 

1064 

4 

2 

55 


APPENDIX  B 

DOCUMENTATION  FOR  COMPUTER  PROGRAM 
SINGLE  ARC  ATTRIBUTES 

**************  DOCUMENTATION  ******************* 

This  FORTRAN  program  determines  the  attributes  for  individual  arc  from  the 

100  meter  grid  square  data,  node  characteristics  and  arc  characteristics  of  the  network 

overlayed  on  the  gridded  terrain. 

Assumption  : 

•  If  the  angle  of  inclination  is  between  45  degrees  and  135  degrees,  then  use  X- 
axis. 

•  If  the  angle  of  inclination  is  greater  than  135  degrees  or  less  than  45  degrees, 
then  use  the  Y-axis. 

FUNCTION  XCORD  :  Compute  X  grid  coordinates. 

FUNCTION  YCORD  :  Compute  Y  grid  coordinates. 

SUBROUTINE  INIT  :  Set  the  initial  conditions. 

SUBROUTINE  INPUT1:  Read  the  data  from  data  1,  data  2,  and  data  3. 
SUBROUTINE  INPUT2:  Read  the  data  about  the  arc. 

SUBROUTINE  DST  :  Compute  distance  along  arc. 

SUBROUTINE  SLOP1  :  Compute  slopes  between  each  pair  of  point  along  arc. 
SUBROUTINE  SLOP2  :  Compute  slopes  between  each  pair  of  point  off  arc. 
SUBROUTINE  MIND  :  Compute  minimum  boundary  distance  from  the  arc. 
SUBROUTINE  FLOW  :  Compute  flow  rate  for  each  arc. 

SUBROUTINE  TIME  :  Compute  min  traversal  time  for  each  arc. 

SUBROUTINE  PRT  :  Print  the  results  of  calculation. 

*************  VARIABLE  DEFINITION  ************** 

CASE  :  option  of  military  unit  type 
CA  SE  1;  vehicle  unit 
CASE  2;  dismounted  troops 
PRINT  :  option  of  print 

PRINT  1;  print  all  results  of  calculation. 

PRINT  2;  print  the  minimum  distance  off  arc. 

PRINT  3;  print  the  input  data  for  checking. 

PRINT  4;  print  the  coordinates  of  boundary  line. 

PRINT  5;  print  the  data  for  network. 
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UNITF  :  type  of  formation 

UNITF  1;  multiple  column  formation 
UNITF  2;  a  single  column  formation 
THRESHOLD  of  slope  for  boundary  line 

THRES1  :  threshold  for  positive  slope 
THRES2  :  threshold  for  negative  slope 
ARCTP  :  arc  type  from  a  head  node  to  a  tail  node. 

ACTD  :  actual  surface  distance  of  each  pair  of  points. 

AD  :  total  actual  surface  distance  along  the  arc. 

ADIST  :  distance  between  head  node  and  tail  node  (in  kilometers). 
ALTD  :  altitude  of  each  point  (  raw  data  ). 

ALT1  :  altitude  of  the  head  node  on  the  arc. 

ALT2  :  altitude  of  the  tail  node  on  the  arc. 

ANGLE  :  angle  between  arc  and  X-axis  (in  degree). 

AREAL  :  real  variable  by  using  array. 

ASLPCH  :  average  slope  change  along  the  arc. 

ATI  ME  :  minimum  traversal  time  for  each  arc. 

AVGSLP  :  average  slope  along  the  arc. 

DELTAX  :  distance  between  node  1  and  node  2  along  X-axis. 
DELTAY  :  distance  between  node  1  and  node  2  along  Y-axis. 
DIST  :  distance  between  node  l  and  node  2  on  the  map. 

DISTA  :  distance  between  each  pair  of  point  along  arc. 

DISTB  :  distance  between  each  pair  of  point  off  arc. 

DHT  :  difference  between  height  of  each  pair  of  point  along  arc. 
DHT2  :  difference  between  height  of  each  pair  of  point  off  arc. 
DSTN  :  distance  of  each  element  in  line  (in  meters). 

DW  :  doctrinal  width  of  battalion  during  deployment. 

DD  :  doctrinal  depth  of  battalion  during  deployment. 

FORWV  :  minimum  formation  'width  of  vehicle  unit. 

FORWD  :  minimum  formation  width  of  dismounted  troops. 

FR  :  flow  rate  of  the  arc  (battalion  /  hour). 

HT  :  height  of  intercept  between  arc  and  Y-axis  along  arc. 

HT2  :  height  of  intercept  between  arc  and  X-axis  off  arc. 

INTS  :  single  integer  variable. 

INTA  :  integer  variable  by  using  array. 
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IX, IY  :  crossing  points  between  an  arc  and  100  meter  X,  Y  grid  line. 
KNODE  :  number  of  crossing  point  along  the  arc. 

KR  :  coefficient  of  route  availability. 

LMIN  :  minimum  distance  of  the  left  side  off  arc. 

MAXNO  :  max  number  of  elements  that  move  in  line  along  the  arc. 
NARC  :  total  number  of  arc  in  the  data  file. 

NCODE  :  characteristics  of  each  point  (  raw  data  ). 

NCODEX  :  dummy  variable  for  characteristics  of  each  point. 

NODE  :  total  number  of  node  in  the  data  file. 

NODEA  :  identification  number  of  head  node-  on  each  arc. 

NODEB  :  identification  number  of  tail  node  on  each  arc. 

NSPCL  :  variable  of  using  width  determination  rule  2  for  boundary. 
NX  :  dummy  variable  for  altitude  of  each  point  (  raw  data  ). 
RMIN  :  minimum  distance  of  the  right  side  off  arc. 

ROADW ' :  width  of  road  itself. 

SLOPE1  :  slope  between  each  pair  of  point  along  arc. 

SLOPE2  :  slope  between  each  pair  of  point  off  arc. 

SLPCH  :  slope  change  between  each  pair  of  point  along  the  arc. 

.  SREAL  :  single  real  variable. 

SP  :  movement  speed  for  an  element  of  unit. 

TSLPCH  :  total  slope  change  along  the  arc. 

UTMX  :  U.T.M.  grid  coordinates  of  the  point  (  East  -  West ). 

UTMY  :  U.T.M.  grid  coordinates  of  the  point  (  South-North  ). 

WIDTH  :  width  of  arc  (  road  width  +  off-road  width  ). 

XA, YA  :  intercept  point  between  arc  and  Y-axis. 

XB, YB  :  intercept  point  between  reference  line  and  X-axis. 

XL,YL  :  XfY  coordinates  of  minimum  distance  of  left  side  off  arc. 
XR,YR  :  X,Y  coordinates  of  minimum  distance  of  right  side  off  arc. . 
XNODE  :  X  coordinates  of  each  node. 

YNODE  :  Y  coordinates  of  each  node. 

XM  :  slope  between  X-axis  and  an  arc. 

YM  :  slope  between  Y-axis  and  an  arc. 


CARTESIAN  SPACE  NETWORK 
**************  DOCUMENTATION  ******************* 

This  FORTRAN  program  determines  the  minimum  time  path  and  minimum 
distance  path  from  the  undirected  graph  by  using  a  shortest  path  algorithm.  Input 
data  of  this  program  is  output  data  from  the  single  arc  attributes  program. 
SUBROUTINE  OPTION:  Determine  the  initial  options. 

SUBROUTINE  INIT  :  Set  the  initial  conditions. 

SUBROUTINE  FORMTN:  Specify  the  formation  width  and  depth. 

SUBROUTINE  SHORTP:  Determine  the  minimum  time  or  distance  path. 
SUBROUTINE  PRINT  :  Print  the  results  of  calculation. 

*************  VARIABLE  DEFINITION  ************** 

CASE  :  Option  of  Military  Unit  Type 
CASE  1;  vehicle  unit 
CASE  2;  dismounted  troop  unit 
UNITF  :  Type  of  Formation 

UNITF  1;  multiple  column  formation 
UNITF  2;  a  single  column  formation 
PARAM  :  Parameter  for  Shortest  Path 

PARAM  1;  minimum  time  path 
PARAM  2;  minimum  distance  path 
BIG  :  an  arbitrary  large  real  number. 

DIST  :  distance  between  head  node  and  tail  node  on  the  map. 

FLOW  :  flow  rate  of  each  arc. 

HEAD  :  head  node  of  each  arc. 

INODE  :  tail  node  of  each  arc. 

INTS  :  single  integer  variable. 

INTA  :  integer  variable  by  using  array. 

ISPATH  :  integer  vector  of  length  N. 

ISTART  :  start  node  of  the  minimum  time  or  distance  path. 

IWORK2  :  integer  vector  of  length  N. 

IWORK3  :  integer  vector  of  length  M. 

JNODE  :  head  node  of  each  arc. 

LARC  :  arc  number  along  the  minimum  parameter  path. 

LAST  :  terminal  node  of  the  minimum  time  or  distance  path. 

LINK  :  arc  number  in  the  graph  ( link(head  node,  tail  node) ). 
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M  :  number  of  arcs  (  edges  )  in  the  graph. 

MINE  :  arc  number  which  is  on  the  minefield. 

N  :  number  of  nodes  in  the  graph. 

NARC  :  total  number  of  arc  in  the  data  file. 

NODE  :  total  number  of  node  in  the  data  file. 

NO  :  arc  number  in  the  graph  (  one  dimensional  array  ). 

NP  :  the  integer  NP  has  the  value  zeio  if  a  shortest  path  is  found 

between  ISTART  and  LAST,  otherwise  it  has  the  value  one 
NUMP  :  the  number  of  nodes  in  the  shortest  path  found  between 
ISTART  and  LAST. 

PATHL  :  total  distance  between  ISTART  and  LAST. 

REALS  :  single  real  variable. 

REALA  :  real  variable  by  using  array. 

RESM  :  character  variable  RESM  has  the  value  Y  if  there  is  a  mine 
field  in  the  graph,  otherwise  it  has  the  value  N. 

TAIL  :  tail  node  of  each  arc. 

TIME  :  minimum  traversal  time  of  each  arc. 

TPATH  :  total  traversal  time  for  the  entire  unit  between  ISTART 
and  LAST. 

WIDTH  :  width  of  arc  (  road  width  +  ofF-road  width ). 

WK4  :  real  vector  of  length  N. 

XLEN  :  total  value  of  parameter  between  ISTART  and  LAST. 
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*•  * 


APPENDIX  C 

COMPUTER  PROGRAM  FOR  SINGLE  ARC  ATTRIBUTES 


This  appendix  contains  the  computer  program  used  to  support  the  algorithms  in 

Chapter  III. 

PROGRAM  TSIM 
JULY  5,  1987  (20*00) 

*************  VARIABLE  DECLARATION  ************* 

INTEGER  ALTO , ARCTP , CASE , CODE , PRINT , UTMX , UTMY , UNITE 

REAL  LMIN 

CHARACTER  *1  BLANK 

DIMENSION  ALTD(100 .100) ,NC0DE(10O ,100) ,NX(10) ,NC0DEX(10) 
DIMENSION  XN0DE(90  ,YNODE(90) ,C0DE(0»7,0«7) 

DIMENSION  DISTA  90) ,XA(90) ,YA(90) 


DIMENSION  HT(90) ,DHT(9Q) ,SLOPE1(90) 

DIMENSION  DISTBUOO.IOO)  ,XB(100,106) ,YB(100.100) 
DIMENSION  HT2?100 ,100) .DHT2(100,1G0) ,SL0PE2(1Q0 ,: 
DIMENSION  N0(90),SLPCH(90J 
DIMENSION  XL(loj,YL(10),XR(10),YR(10' 


4Md  \  *V  I  <  \  AW  f  t  \  4W  /  *  *i\\AW  f 

DIMENSION  SP(7 .3 , 2) ,FR(5 , 5) ,MAXNO(3 . 3) 

DIMENSION  ACTD(90),ATIME(3 
DIMENSION  ROADH(7) 

COMMON  /  INTS  /  NSTEP  ,1X1, 1 Y1 .  1X2 , 1 Y2 , JA , JB , K ,  NARC ,  NCOUNT , 

NODEA ,  NODES  ,  J2 ,  J3  ,  J4 ,  PRINT .  UTMX ,  UTMY , 

CASE , ARCTP , FR , AD . UNITF , NSPCL 
COMMON  /  INTA  /  NCODE, NO, ALTO, CODE, MAXNO 

COMMON  /  SREAL  /  XM , ANGLE , DIST , TSLPCH, ASLPCH , ALT1 , ALT2 , AVGSLP , 

,  LMIN, RMIN, XL, YL.XR.YR.YM.RY, WIDTH. THRES1 ,THRES2 
COMMON  /  AREAL  /  XA , YA , HT , DHT , DISTA , DISTB , SLPCH , SL6PE 1 , SLOPE2 , 

XB , YB , HT2 , DHT2 , SP , ACTD , ATIME , ROADW 


*  1.  SET  THE  INITIAL  OPTIONS. 

* 


WRITE (6, 601) 

FORMAT ( 1  ',////////) 

PRINT*,  '  WHICH  TYPE  OF  UNIT  DO  YOU  DESIRE?  (  1  OR  2  ) ' 

PRINT*.  ' 


PRINT*,  1 
PRINT*,  ' 
PRINT*,  ' 
PRINT*,  ' 
PRINT*,  ' 
PRINT*,  ' 
PRINT*,  ' 
READ? 5,*) 
PRINT*,  ' 


CASE 
**  NOTE 
CASE 


******************************* i 

*  *  i 

*  1.  VEHICLE  UNIT  *' 

*  *  i 

*  2.  DISMOUNTED  TROOPS  *' 

*  *  i 
******************************* i 

:  FOR  YOUR  REFERENCE,  CURRENT  ANSWER  IS1 


IF  (  CASE  .EQ.  1  .OR.  CASE  .EQ.  2  )  THEN 
GO  TO  17 

ELSE  IF  (  CASE  .NE.  1  .OR.  CASE  .NE.  2  )  THEN 
PRINT*,  '  ' 

PRINT*,  '***  ERROR  «  ENTER  THE  NUMBER  1  OR  2  ***' 

PRINT*,  '  ' 

GO  TO  16 
END  IF 
WRITE(6f602) 

FORMAf(f ,',/////////). 

PRINT* ,  1  1 

PRINT*,  1  WHICH  TYPE  OF  FORMATION  DO  YOU  DESIRE?' 
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*  *  *  *  *  .  *  * 


18 

603 


604 


k 

19 


PRINT* 

PRINT* 

PRINT* 

PRINT* 

PRINT* 

PRINT* 

PRINT* 

PRINT* 

READ?5 

PRINT* 


FOR  YOUR  REFERENCE ,  CURRENT  ANSWER  IS1 


************************************* i 

*  *i 

*  1.  MULTIPLE  COLUMN  FORMATION  *• 

*  *i 

*  2.  SINGLE  COLUMN  FORMATION  *' 

*  *  • 

************************************* • 

I  UNITF 

**  mote 

UNITF 

IF  (  UNITF  .EQ.  1  .OR.  UNITF  .EQ.  2  )  THEN 
00  TO  18 

ELSE  IF  (  UNITF  .NE.  1  .OR.  UNITF  .NE.  2  )  THEN 
PRINT* ,  '  1 

PRINT*,  '***  ERROR  s  ENTER  THE  NUMBER  1  OR  2  ***' 

PRINT*,  '  ' 

GO  TO  17 
END  IF 
WRITE(6.603) 

FORMAT(f  ',///////,/) 

PRINT* ,  '  1 

PRINT*,  1  DO  YOU  WANT  TO  USE  SPECIAL  RULE  FOR  BOUNDARY?' 

PRINT*,  '  YES  (  1  )  OR  NO  (  2  )' 

READ? 5,*)  NSPCL 
PRINT*,  '**  NOTE 
NSPCL 
WRITE (6 . 604) 

FORMAT '  ,/////////////) 

PRINT*,  1  1 

WHICH  OUTPUT  DO  YOU  WISH  TO  PRINT?  (  1  -  5  ) 1 


FOR  YOUR  REFERENCE,  CURRENT  ANSWER  IS1 


PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 


******************************************** ■ 

*  1.  GENERAL  INFORMATION  FOR  EACH  ARC  *' 

*  2.  FLOW  RATE  &  TRAVERSAL  TIME  *' 

*  3.  INPUT  DATA  FOR  CHECK  *' 

*  4.  COORDINATES  OF  BOUNDARY  LINE  *' 

*  5.  NETWORK  DATA  FOR  NETWORK  PROGRAM  *' 

*************************************** 


*********** • 
>  FOR  YOUR  REFERENCE,  CURRENT  ANSWER  IS1 , 


.EQ.  1  .OR.  PRINT  .EQ.  2  .OR.  PRINT 
"  .EQ.  4  .OR.  PRINT  .EQ.  5  )  THEN 


EQ.  3 


READ? 5,*)  PRINT 
PRINT*,  ’**  NOTE 

*  PRINT 
IF  (  PRINT 

*  .OR.  PRINT 
GO  TO  19 

ELSE  IF  (  PRINT  .NE.  1  .OR.  PRINT  .NE.  2  .OR.  PRINT 

*  .OR.  PRINT  .NE.  4  .OR.  PRINT  .NE.  5  )  THEN 
PRINT*,  1 

PRINT*,  '***  ERROR  ;  ENTER  THE  NUMBER  1  THRU  5  ***' 
PRINT* ,  ' 

GO  TO  18 
END  IF 

THRESHOLD  FOR  BARRIER  DUE  TO  SLOPE  OFF  ROAD 
THRES1  =  0.3 
THRES2  =  -0.3 


.NE.  3 


2.  SET  THE  INITIAL  CONDITIONS. 


TOTAL  NUMBER  OF  ARCS  IN  NETWORK. 

NARC  =  72 

TOTAL  NUMBER  OF  NODES  IN  NETWORK. 

NODE  »  34 

t  CALL  INIT (NX , NCODEX , XL , YL , XR , YR , XA , YA , SLPCH , HT , NO , DHT , CASE , 

*  DISTA , SLOPE 1 , XB , YB , HT2 , DHT2 , DISTB , SL0PE2 , CODE ) 

3.  READ  THE  DATA  FROM  THE  DATA  FILE. (DATA  1,  DATA  2,  DATA  3,  DATA  4) 
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\ 


•yjjj-i  \  i^i  >n  >rr  rm  mi\  Iftw  art  wutoJtfuwtfifiis 


\JW 


««««  *  *  **  *  **  *  ** 


DATA  1  ;  ALTITUDE  AND  CODS 
DATA  2  ;  NODS  CHARACTERISTICS 
DATA  3  ?  ARC  CHARACTERISTICS 
DATA  4  ;  SPEED  FOR  MOVEMENT 
CALL  INFUT1 (NARC .NODE , ALTD .NCODB, XNODE , YNODE, 

UNITF , CAS4 , THRBS1 , THRES2 , SP , PRINT) 

NCOUNT  ■  1 

IY  -  100 

DO  10  IB  ■  1.  NARC 

CALL  INPUT2 (XNOD* , YNODE , IY , 1X1 , IY1 , 1X2 , IY2 , NODEA , NODEB , 
ARCTP , UTMX , UTMY ) 


4.  DETERMINE  THE  ANGLE  BETWEEN  ARC  AND  X-AXIS. 


DBLTAX  *  1X2  -  1X1 

DELTAY  *  IY2  -  IY1 

IF  (  DELTAX  .NE.  0.)  THEN 

XM  ■  DELTAY  /  DELTAX 

ELSE  IF  (  DELTAX  .EO.  0.)  THEN 

XM  »  9999999. 

END  IF 

RADIAN  *  ATAN(  XM  ) 

DIST  *  SQRT(  DELTAX  **  2  +  DELTAY  **  2  ) 
ANGLE1  *  (RADIAN  *  180.)  /  3.141592 
IF  (  ANGLE 1  ,GE.  0.)  THEN 
ANGLE  »  ANGLE 1 
ELSE  IF  (  ANGLE 1  .LT.  0.)  THEN 

ANGLE  *  180.  +  ANGLE 1 
END  IF 


5.  COMPUTE  THE  DISTANCE  BETWEEN  EACH  PAIR  OF  POINTS  ALONG  ARC. 


CALL  DST ( IY1 , 1X1 , IY2 , 1X2 , ANGLE ,DIST,XM,XA,YA,DISTA,K) 


6.  COMPUTE  AVERAGE  SLOPE  ALONG  THE  ARC. 


•LX1  *  1X1  /  100 

LY1  ■  IY1  /  100 

LX2  »  1X2  /  100 

LY2  ■  IY2  /  100 

ALT1  »  ALTD(  LX1,  LY1  ) 

ALT2  »  ALTD(  LX2,  LY2  ) 

AVGSLP  ■  (  ALT2  -  ALT1  )  /  DIST 


7.  COMPUTE  SLOPES  BETWEEN  EACH  PAIR  OF  POINTS  ALONG  THE  ARC. 


CALL  SLOP1 ( ANGLE . ALT1 , ALT2 , 1X1 , IY1 , IY2 , 1X2 , DISTA , XM , XA , YA , ALTD , 
*  HT, DHT, SLOPE 1,NA,ACTD, AD) 


*  8.  AS  A  MEASURE  OF  TERRAIN  "  STEEPNESS  "  ALONG  THE  ARC,  COMPUTE 

*  TOTAL  AND  AVERAGE  SLOPE  CHANGE. 


ASLPCH  «  0. 

TSLPCH  *  0. 

NB  a  NA  -  1 
DO  11  ID  »  1  NB 

SLPCH(ID)  =  ABS(  SL0PE1(ID+1)  -  SLOPEI(ID)) 
TSLPCH  «  TSLPCH  +  SLPCH(ID) 

11  CONTINUE 


! 


] 

j 


KIN  DISTANCE 


*  9.  COMPUTE  SLOPES  FROM  ARC  TO  POINT  OFF  ARC,  AND 


* 

CALL  SLOP2 

* 

A 

10.  DETERMINE  COORDINATES  FOR  THE  KIN  DISTANCE . (  LEFT  AND  RIGHT  ) 

A 

CALL  MIND(K,XM,YM,XA,YA,LMIN,RMIN,XL,YL,XR,YR) 

A 

A 

11.  COMPUTE  THE  RATE  OF  FLOW  ALONG  THE  ARC. 

A 

CALL  FLOW(SP, WIDTH, ARCTP, CASE, FR,MAXNO) 

* 

A 

12.  COMPUTE  THE  MIN  TRAVERSAL  TIME  ALONG  THE  ARC. 

A 

CALL  TIME ( ARCTP , SP , UNITF , CASE , DIST , WIDTH , ATIME ) 

* 

* 

13.  PRINT  THE  RESULT  OP  CALCULATION. 

CALL  PRT 
NCOUNT  »  NCOUNT  +  1 
10  CONTINUE 
STOP 
END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAitAAAAAAAAAAAAAAAAAAAAAAA 

*  A.  SUBROUTINE  FOR  SETING  THE  INITIAL  CONDITIONS. 


SUBROUTINE  INIT (NX .NCODEX ,  XL ,YL,XR, YR.XA , YA , SLPCH , HT , NO .DHT , CASE 
*  DISTA , SLOPE 1 ,XB, YB,HT2, DHT2 , DISTB , SLOPE 2 , CODE ) 


INTEGER 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


YR(10) ,NX(10) .NCODEX (10) 


CASE, CODE 
XL(  ‘ 

DIJ 
HT( 

DisTBUoo'i66);3a(ioo;i66),YB(ioo.ioo) 

HT2? 100 , 100 ) . DHT 2 (100,100), SLOPE2 ( 100 , 100 ) 
NO(90),SLPCH(90) 

CODE(6«7,0«7) 


*  INITIALIZATION 
DO  10  L 
XL 
YL 
XR 
YR 
NX 

ncOd: 

10  CONTINUE 

DO  11  LA  > 
XA(LA 
YA(LA 
HT(LA 
NO(LA. 
DHT(LA 
SLPCH 
DISTA 


FOR  VARIABLES. 
-  1,  10 
=  0. 
a  0. 
a  0. 
a  0, 

*  0 

X(L)  =  0 


\ 


i 


I 


0. 


SLOPEl(LA) 

11  GONfXNUI 
DO  12  LB 
DO  12  LC 

XB 
YB 
HTL 
DHT2 

D1S  , _ _ 

SLOPE2(LB,LC> 

12  CONTINUE 

*  SET  CODS  OF  BOUNDARY  FOR  VEHICLE  UNIT. 

IF  7  CASE  .SQ.  1)  THEN 
DO  13  I  *  0,7 

DO  13  J  *  0,7 

CODE(I, J)  *  0 

13  CONTINUE 

CODE (4, 4)  -  1 
CODE (4, 7)  *  1 
CODE (7, 4  *  1 
0008(7,7)  *  1 

*  SET  CODE  OF  BOUNDARY  FOR  DISMOUNTED  TROOPS. 

ELSE  IF  (  CASE  .EQ.  2)  THEN 
DO  14  I  *  0*7 

DO  14  J  *  0,7 

CODB(I,J)  *  1 

14  CONTINUE 

DO  IS  IA  *  0,7 

L  *  0 


CODE 

CODE 

CODE 

IS  CONTINUE 

END  IF 
RETURN 
END 


IA,L 

IA,3 

IA,S 


0 

0 

0 


*  B.  SUBROUTINE  FOR  READING  THE  DATA  FILES (  DATA  1,  DATA  2,  DATA  4  ) 

A********************************************************************** 


100 


22 

21 

101 

20 

* 


SUBROUTINE  INPUT1 (NARC , NODE , ALTD , NCODE . XNODE , YNODE , 

*  UNITF, CASE, THRES1,THRES2,SP, PRINT) 

INTEGER  ALTD, CASE, PRINT, UNITF 

rHiRArrvn  *1  rt inv 

DIMENSION  ALTD (100, 100). NCODE (100 .100) 

DIMENS ION  NX ( 10 ) , NCODEX ( 10 ) , XNODE ( 90 ) , YNODE ( 90 ) 
DIMENSION  SP(7 ,3,2) 

READ  THE  ALTITUDE  AND  CODE  BY  TWO  DIMENSIONAL  ARRAY. 

J  *  0 

DO  20  IN  *  1,99 
1*1 
J  *  J  +  1 
DO  21  IT  *  1.10 

READ(1 .100)  (  NX(I1).  NCODEX(Il) ,  11*1,10  ) 

FORMAT  14,12,9(15,12) J 
DO  22  IS  *  1,10 

ALTD(I.J)  *  NX(IS) 

NCODE ( I , J )  *  NCODEX (IS) 

1*1  +  1 
CONTINUE 
CONTINUE 

READ(l.lOl)  BLANK 
FORMAT  (Al) 

CONTINUE 

READ  THE  NODE  CHARACTERISTICS  FOR  SINGLE  ARC  ATTRIBUTES. 
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110 

23 


130 

26 


716 

769 

710 

711 
24 

712 

713 

714 


DO  23  IA  •  1(  NODS  , 

Mlite5<u>-YH0D,<“) 

CONTINUE 

DO  26  IN  *  1,7 

DO  26  JM  *  1,3 

DO  26  KM  ■  1,2 

READ  THE  SPEED  FOR  MOVEMENT  BY  THREE  DIMENSIONAL  ARRAY. 

RE  AD (4, 130)  SP(IM, JM,KM) 

FORMAT  F5. 6) 

CONTINUE 

IF  (  PRINT  .EQ.  1  )  THEN 

WRITE (6  ,716) CASS , UNITF , THRBS1 . THRES2 
FORMAT ( '  * ,4X, 'UNIT  TYPE  i f , 14,7 ,SX, 'UNIT  FORMATION' , 

*  '  «» /I4  ,/,5X'?OS.  THRESHOLD  t',F7.3, 

*  /,5X,'NEG.  THRESHOLD  *',F7.3  ) 


ELSE  IF 


!•* 

.  * 


(  PRINT  .EO.  2  )  THEN 
WRITE (20,769) CASE , UNITF , TKRES1 . THRES2 
FORMAT ( '  1 ,4X, 'UNIT  TYPE  « f , 14, / ,5X, 'UNIT  FORMATION', 

'  i ' ,14  , / , 5X , ' POS .  THRESHOLD  »',  F7.3, 

/,5X, 'NEC.  THRESHOLD  t',F7.3  ) 

ELSE  IF  (PRINT  .SQ.  3  )  THEN 
WRITE (6 .710) 

FORMAT? '  ' ,3X,71('-'),//,4X,'  INPUT  DATA  ', 

'(  ALTITUDE  AND  CHARACTERISTICS  )',/) 


NX  *  1 
DO  24  NA 
DO  24  NB 


1.2 

1,100 


WRITE?6.71l)  NX.NB.NA,  ALTD(NB.NA) .NCODE(NB.NA) 
FORMAT^  'IS,'  ({,  i3,I3.'  J  ' ,i5,I5) 

NK  *  NX  +  1 
CONTINUE 

WRITE (6  712) 

FORMAT ( |  ' ,3X,// ,4X, 1  NODES  ALONG  ARC  ', 

ELSE  IF  (  PRINT  .EQ.  5  )  THEN 
WRITE?6 . 714) 

FORMAT ( '  '  ,oX, 'HEAD  TAIL',3X," 


.  .  .  .'TIME  FLON  RATE  DIST  WIDTH  SPEED') 

END  IF 
RETURN 
END 

*********************************************************************** 


*  C.  SUBROUTINE  FOR  READING  THE  DATA  FILE  (  DATA  3  ) 
*********************************************************************** 

SUBROUTINE  INPUT2 ( XNODE , YNODE . IY , 1X1 , IY1 , 1X2 , IY2 , NODEA , NODEB , 

*  ARCTP.UTMX.UTKY) 

INTEGER  ARCTP , UTMX , UTHY 

DIMENSION  XNODE(90) ,YN0DE(90) 

*  READ  THE  ARC  CHARACTERISTICS  FOR  SINGLE  ARC  ATTRIBUTES. 

READ(3.120)  NODEA, NODEB, ARCTP 
120  FORMAT (13, 14, 14) 

NX1  =  XNODE (  NODEA  ) 

NY1  *  YNODE (  NODEA  ) 

UTMX*  NX1  /  100 
UTMY*  NY1  /  100 
UTMX*  UTMX  *  10 
UTMY*  UTMY  *  10 

MX1  *  MOD(  NXl.IY  )  *  100 
MY1  *  MOD(  NY1,IY  )  *  100 
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NX2  -  XNODE(  NODEB  ) 

NT2  *  YNODE(JIODEB  ) 

HX2  «  HOO(  NX2.IY  )  *  100 
H! n  -  H00rHT2.1T  >  *  100 

IF  (  HT2  .OT.  HTI)  THIN 
1X1  *  HX1 
IT1  ■  HTI 
1X2  «  MX2 
1Y2  ■  MT2 

ZHANGE  THE  ORDER  OF  NODE  FOR  USING  ENGLES  OF  INCLINATION. 


CHANGE  T1 
ELSE  IF 


THE  ORDER  OF 
F7HT2  .LT. 
ixi  «  NX2 


THEN 


m  ■ 

1X2  * 
IT2  « 
N00E1 
NODE  2 
NODEA 


ELSE  IF 
IF 


HT2 

MX1 

HTI 

»  NODEA 

■  NODEB 
*  N0DE2 

■  N0DE1 

!  .10.  HTI) 
!  .GE.  HX1) 


THEN 

THEN 


IT1  ■ 
1X2  ■ 
IY2  * 
ELSE  IF 
1X1  * 

:yi  - 

1X2  « 
IY2  » 
N0DE1 
N0DE2 
NODEA 
NODEB 
END  IF 
IF 


HTI 

HX2 

^HX2  .LT.  HX1)  THEN 

HT2 

HX1 

HTI 

■  NODEA 
-  NODEB 
•  N0DE2 
«  N0DE1 


END  IF 
RETURN 

*  D.  FUNCTION  FOR  CALCULATING  GRID  COORDINATES (  X  ). 


FUNCTION  XCORD(XK.IA.IY.IX) 

IF  (  XH  .NE.  0.)  THEN 

XCORD  »  ((  1,/  XH  )  *  ( 
ELSE  IF  (  XH  .EQ.  0.)  THEN 


IA  -  IT  ))«•  IX 


XCORD  «  IX 
END  IF 
RETURN 
END 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa*aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
*  E.  FUNCTION  FOR  CALCULATING  GRID  COORDINATES (  Y  ). 


FUNCTION  YCORD(YH.XXB.YC.XC) 

YCORD  »  (  YH  *  (  XXB  -  XC  ))+  YC 

RETURN 

END 


*  F.  SUBROUTINE  FOR  CONFUTING  THE  DISTANCE  BETWEEN  EACH  PAIR 

*  OF  POINTS  ALONG  ARC. 


30 

•k 


31 

k 


32 


SUBROUTINE  DST ( IY1 , 1X1 , IY2 , 1X2 , ANGLE ,DXST,XM,XA,YA, 
*  DIS*A,x) 


DIMENSION  DISTA(90) ,XA(90) ,YA(90) 


.AND. 


- \ 

NA  «  ( 

(IY2  - 

IY1  +  49 

ELSE 

IF  (  ANGL) 

i  .GT. 

135.  .OR 

NA45  =  7 

(1X2  - 

1X1  +  49 

NA135  «  { 

1X1  - 

1X2  +  49 

IF(  ANGLE 

.LT. 

45.)  NA 

IF(  ANGLE 

.GT. 

135.)  NA 

END 

IF 

K 

a  1 

YA(K 

)  a  IY1 

XA(K 

=  1X1 

LE  ,LE.  135.  ) 

)  /  100)  +  1 
.  ANGLE  .LT.  45. 


! 


/  100) 
/  100) 
NA45 
NA135 


THEN 
)  THEN 


ANGLE  OF  INCLINATION  IS  BETWEEN  45  AND  135  DEGREES. 
IF  (  ANGLE  .GE.  45.  .AND.  ANGLE  .LE.  135.  )  THEN 
DO  30  IA  *  0,  IY2-1Q0 ,  100 
A  =  IA  -  IY1 
IF  (A  .GT.  0.)  THEN 
K  =  K  +  1 
IY  *  IY1 
IX  =  1X1 

YA(K)  =  REAL(  IA  ) 

XA(k)  =  XCORD (  XM,IA,IY,IX  ) 

END  IF 
CONTINUE 

ANGLE  OF  INCLINATION  IS  GREATER  THAN  135  DEGREES. 
ELSE  IF  (  ANGLE  .GT.  135.  )  THEN 
'  1X1  +  100  )  /  100 
1X2  /  100  )  +  1 
il  IA  =  11,12,-1 
K  *  K  +  1 
XA(K)  -  IA  *  100. 

YC  »  IY1 
XC  =  1X1 
YM  a  XM 
XXB  =  XA(K) 

YA(K)  ■  YCORD(  YM,  XXB,  YC,  XC  ) 

CONTINUE 

ANCLE  OF  INCLINATION  IS  LESS  THAN  45  DEGREES. 

ELSE  IF  (  ANGLE  .LT.  45.  )  THEN 

11  =  (  1X1  +  100  )  /  100 

12  =  (  '  *“ 


1X2  /  100  )  -  1 

DO  32  IA  =  II,  12,  1 


END 


CONTINUE 
IF 


II, 

K  =  K  +  1 
XA(K)  a  IA  *  100. 

YC  =  IY1 
XC  =  1X1 
YM  =  XM 
XXB  =  XA(K) 

YA(K)  =  Y CORD (  YM,  XXB,  YC,  XC  ) 


COMPUTE  THE  DISTANCE  BETWEEN  EACH  PAIR  OF  POINTS. 
K  =  K  +  1 
IY  »  IY1 
IX  =  1X1 
IA  =  IY2 

IF(  XM  .NE.  0.  )  THEN 

YA(“* 

XA< 

ELSE  IF ( 

YA( 

XA( 

END  IF 
DIST1  =  0. 


•  W  I 

(K)  a  IY2 

K)_  -  XCORD (  XM,  IA,  IY,  IX  ) 


XM 


...I  .EQ. 
K)  =  IY2 
K)  =  1X2 


0.  )  THEN 
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XA(IB+1)  -  XA(IBj j  **2 


^YA^IB+1)  -  YA(IB 
DIST'-  DISH 


**2  ) 


DO  33  IB  ■  1,  HA 

IF(  IB  .LT.  NA  )  THEN 
DISTA(IB)  ■  SQRT(( 

ELSE  IF(  IB  .EQ.  NA  ) 

DISTA(IB)  ■  r — 

END  IF 

DIST1  *  DIST1  +  DISTA(IB) 

33  CONTINUE 
RETURN 
END 

*********************************************************************** 


*  G.  SUBROUTINE  FOR  COMPUTING 

*  ALONG  THE  ARC. 


SLOPES  BETWEEN  EACH  PAIR  OF  POINTS 


*********************************************************************** 


SUBROUTINE  SLOP1 (ANGLE , ALT! , ALT2 , 1X1 , IY1 , IY2 , 1X2 , DISTA , XM , XA , YA , 

*  ALTD  ,HT , DHT , SLOPE1 ,NA, ACTD , AD) 

INTEGER  ALTD 

DIMENSION  ALTD(100 , 100) ,DISTA(90) ,XA(90) ,YA(90) 

DIMENSION  HT(90) ,DHT(9Q) , SLOPEl (90) 

DIMENSION  ACTD (90) 

AD  =0. 

RY  =»  100. 

HT(1)  *  ALT1 

IF  (  ANGLE  .GE.  45.  .AND.  ANGLE  .LE.  135.  )  THEN 
NA  »  ( (IY2  -  IY1  ••  49  )  /  100)  +  1 
DO  40  IC  *  2,  NA 

JX1  *  INT(  XA(IC)  /  100  ) 

JY1  =  INT(  YA(IC)  I  100  ) 

R  *  MOD(XA(IC),RY) 

HTA  =  (R/100. )  *  (  ALTD( JX1 , JY1)  -  ALTD(JX1+1, JY1)) 

IF(  IC  .LT.  NA  )  THEN 

HT(IC  )  »  ALTD(JXl.JYl)  -  HTA 
ELSE  IF(  IC  .EQ.  NA  )  THEN 
HT(IC  )  *  ALT2 
END  IF 

DHT(IC-l)  »  HT(IC  )-  HT(IC-l) 

SLOPEl (IC-1)  =  DHT(IC-1)  7  DISTA(IC-l) 

ACTD(IC-l)  =  DISTA(IC-l)  /  COS(SLOPEl(IC-l) ) 

AD  8  AD  +  ACTD (IC-1) 

40  CONTINUE 

ELSE  IF  (  ANGLE  .GT.  135.  .OR.  ANGLE  .LT.  45.  )  THEN 
NA45  a  ((1X2  -  1X1+  49  )  /  100)  +1 

NA135  =  ((1X1  -  1X2+  49  )  /  100)  +  3 

IF (  ANGLE  .LT.  45.)  NA  =  NA45 
IF(  ANGLE  .GT.  135.)  NA  =  NA135 
DO  41  IC  =  2,  NA 

JX2  =  XA(IC)  /  100 

JY2  »  yA(IC)  /  100 

R  =  MQD(YA(IC) ,RY) 

HTA  a  (R/100.)  *  (  ALTD ( JX2 , JY2 )  -  ALTD(JX2, JY2+1) ) 

IF(  IC  .LT.  NA  )  THEN 

HT(IC)  a  ALTD(JX2,0Y2)  -  HTA 
ELSE  IF (  IC  .EQ.  NA  )  THEN 
HT(IC)  a  SLT2 
END  IF 

DHT (IC-1)  =  HT(IC)-  HT(IC-l) 

SLOPEl(IC-l)  =  DHT(IC-l)  /  DISTA(IC-l) 

ACTD(IC-l)  =  DISTA(IC-l)  /  COS (SLOPEl ( IC-1 ) ) 

AD  =  AD  +  ACTD(IC-l) 

41  CONTINUE 
END  IF 
RETURN 
END 

*********************************************************************** 
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*  H.  SUBROUTINE  FOR  COMPUTING  SLOPES  FROM  ARC  TO  POINTS  OFF  ARC. 


SUBROUTINE  SLOP2 


INTEGER  ALTD , ARCTP , CASE , CODE , PRINT , UTMX , UTMY 

REAL  LMIN 

DIMENS ION  ALTD (100 , 100 ) , NCODE ( 100 , 100 ) , NX ( 10 ) ,NCODEX( 10 ) 
DIMENSION  XNODE(90 ) , INODE (90) , C0DE(Q  j7 ,0 :7) 

DIS1A(90) ,XA(90),YA(90) 

HT[90) ,DHT(90) , SLOPE1 (90) 

STB(100,100  ,XB(100,100) 
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DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


DISTB( 100 , 100) ,XB{100,100),YB(100,100) 

HT2 ( 100 , 100 ) , DHT2 ( 100 , 100 ) , SLOPE2 ( 100 , 100) 
NO(90) ,SLPCH(9Q) 

XL (10) ,YL(10) ,XR(10) ,YR(10) 
SP(7,3,2),FR(J,5), MAXNO (3,5) 

ACID  ( 90 ) ,  ATIME  ( 3 ) 

ROADW ( 7 ) 


1  COMMON  /  INTS  /  NSTEP ,  1X1 , 1 Y1 ;  1X2 ,  IY2 ,  JA ,  JB ,  K  ,NARC ,  NCOUNT , 

*  NODEA.NODEB , J1 , J2 , J3 , J4 , PRINT  rUTMX,  UTMY , 

*  CASE, ARCTP, FR, AD, UNITF.NSPCL 
COMMON  /  INTA  /  NCODE , NO , ALTD , CODE , MAXNO 

COMMON  /  SREAL  /  XM ,  ANGLE ,  DIST ,  TSLPCH ,  ASLPCH ,  ALT! ,  ALT2 ,  AVGSLP , 

*  LMIN, RMIN, XL, YL ,XR,YR,YM,RY, WIDTH, THRES1 ,THRES2 
COMMON  /  ARaAL  /  XA , YA , HT , DHT , DISTA , DISTB , SLPCH , SLOPE 1 , SLOPE2 , 

*  XB,YB,HT2,DHT2,SP,ACTD, ATIME, ROADW 

LMIN  =  10000. 

RMIN  =  10000. 

RY  =  100 . 

IF  (  XM  .NE.  0.)  THEN 
YM  *  -  1.  /  XM 

ELSE  IF  (  XM  .EO.  0.)  THEN 
YM  «  9999999. 

END  IF 

ANGLE  OF  INCLINATION  IS  BETWEEN  45  AND  135  DEGREES. 

IF  (  ANGLE  .GE.  45.  .AND.  ANGLE  .LE.  135.  )  THEN 

COMPUTE  THE  SLOPE  ON  THE  LEFT  SIDE . 

DO  70  IE  =  1,K 
KD  =  1 

IH  *  (INT(XA(IE)  /  100)  *  100)-  100 
IH1  *  IH 
XB(IE,KD)  *  IH 
YC  =  YA(IE) 

XC  »  XA(IE) 

XXB  =  XB(IE,KD) 

YB ( IE , KD )  =  YCORD(  YM,  XXB,  YC,  XC) 

J1  a  XB(IE,KD)  /  100 

J2  =  YB(IE,KD)  /  100 

RA  =  (MOD (  YB(IE ,KD) ,  RY))  /  100. 

IF  (  XM  .GE.  0.  )  THEN 

IF (  RA  .GE.  0.5  )  J4  a  J2  +  1 
IF(  RA  .LT.  0.5  )  J4  =  J2 
ELSE  IF  (  XM  .LT.  0.  )  THEN 
IF(  RA  .GE.  0.5  )  THEN 
J4  =  J2 

ELSE  IF(  RB  .LT.  0.5  )  THEN 
J4  a  J2  -  1 

END  IF 
END  IF 

IF  (  XB ( IE , KD)  .GT.  10000.  .OR.  YB(IE ,KD) 

KD  =  KD  -1 
GO  TO  75 
END  IF 

DISTB(IE,KD)  =  SQRT( (  XB(IE,KD)  -  XA(IE))  **2 


GE.  10000.  )  THEN 


70 


,  +  (  YB(IK.KD)  -  YA(XB)>  **2) 

HTB  ■_  RA.  *  (ALTD j J1 , J||  -  ALTD(J1, J2+1) } 


[E,KD)  -  ALTD(J1,J2)  -  HTB  v 
IS.KD)  ■  KT2(XE,KD)  -  HT(IE) 
ilsfB(IE.KD)  .EO.  0.)  GO  TO  7C 


75 


* 

71 

74 


HT2( 

DHTJ 

IF?  DISTB(IE.KD)  .EO.  , 

SL0PE2CIE ,KD)  ■  DHTZCiMD)  /  DXSTB(IZ,KD) 

XF(  KD  .EO.  1  )  THEN 
HCA  *  7 

ELSE  IF(  KD  .GT.  1  )  THEN 
JA  ■  XB(IE ,KD-l )  /  100 
JB  *  YB(IE.KD-l)  /  100 
NCA  «  NCODE (  JA, JB) 

END  IF 

NCB  ■  NCODE (J1,J4) 

IF  (  CODE (NCA .NCB)  .EQ.  1  .AND.  SL0PE2(IE,KD)  .LE.  THRES1 

*  .AND.  SLOPE2(IE,kD)  .GE.  THRES2  .AND.  KD  .LE.  10)  THEN 
IH  ■  IH  -  100 

KD  *  KD  +  1 
GO  TO  72 

ELSE  IF  (CODS(  NCA, NCB  )  .NE.  1  .OR.  SLOPE2(IE ,KD)  .GT.  THRES1 

*  .OR.  SL0PE2(XB,KD)  .LT.  THRES2  .OR.  KD  .GT.  10  )  THEN 
IH  *  IH1  +  200 

IF(DISTB(IE,KD)  .LT.  LMIN)  LMIN  *  DISTB(IE.KD) 

APPLY  SPECIAL  RULE (WIDTH  DETERMINATION  RULE' 2).  FOR  BOUNDARY. 

IF (  NSPCL  .EO.  1  )  THEN 

IF(  KD  .EQ.  1  )  THEN 
LMIN  *  0. 

END  IF 

END  IF 
GO  TO  71 
END  IF 

COMPUTE  THE  SLOPE  ON  THE  RIGHT  SIDE. 

KD  *  KD  +  1 
XB(IE,KD)  *  IH 
YC  =  YA(IE) 

XC  »  XA(IE) 

XXB  =  XB(IE,KD) 

YB(IE,KD)  *  YCORD(  YM,  XXB,  YC,  XC  ) 

J1  ■  XB(IE,KD)  /  100 

J2  «  YB(IE.KD)  /  100 

RB  a  (MOD  (YB(IE,KD),  RY))  /  100. 

D2  =  ABS (  YA(IE)  -  YB(IE,W))  ) 

IF  (  3M  .GE.  0.  )  THEN 

IF(  RB  .GE.  6.5  )  J4  =  J2  +  1 
IF(  RB  .LT.  0.5  .AND.  D2  .GT.  50.)  J4  a  J2 
ELSE  IF  (  XM  .LT.  0.  )  THEN 

IF(  RB  .GE.  0.5  )  J4  =  J2  +  1 
IF(  RB  .LT.  0.5  ).  J4  =  J2 
END  IF 

IF(XB(IE,KD)  .GT.  10000.  .OR.  YB(IE,KD)  .GE.  10000.  )  THEN 
IF(DISTB(IE,KD-1)  .LT.  RMIN)  THEN 
RMIN  «  DISTB(IE,KD-1) 

END  IF 
GO  TO  77 
END  IF 

DISTB(IE,KD)  =  SQRT ( (  XB(IE,KD)  -  XA(IE))  **2 

*  +  (  YB(IE,KD)  -  YA(IE))  **2) 

HTB  =  RB  *  (ALTD(J1.J2)  -  ALTD(J1,J2+1)) 

HT2(IE,KD)  *  ALTD(J1,J2)  -  HTB 
DHT2(IE,KD)  =  HT2(IE,KD)  -  HT?IE) 

IF (  DISTB(IE.KD)  .EO.  0.)  GO  TO  75 
SLOPE2(IE ,KD)  *  DHT2(IE ,KD)  /  DISTB(IE ,KD) 

IF  (  IH  .EQ.  IHH-260  )  THEN 
NCA  =  7 

ELSE  IF (  KD  .GT.  1  )  THEN 
JA  »  XB(IE ,KD-i)  /  100 
JB  a  YB(IE,KD-1)  /  100 
NCA  ■  NCODE (  JA, JB) 

END  IF 
*  NCODE (J1,J4) 


NCB 


71 


76 


IF  (  CODS (NCA .NCB)  .EQ.  1  .END.  SLOPE2 (IE ,KD)  .LE.  THRES1 
.AND.  SLOPE27lE,KB)  .OE  THRBS2  .AND.  KD  .LE.  20)  THEN 
IH  ■  IH  100 


KD  -  KD  +  1 
GO  TO  74 


ELSE  IF  ( CODE (NCA ,NCB)  .NE.  1  .OR.  SL0PE2(IE ,KD)  .GT. 
.OR.  SL0PE2(IE.KD)  .LT.  THRES2  .OR.  KD  .GT.  20) 

-  .LT.  RMIN) - 

)ISTB(IE,KD) 


i/  • 

.  W*\.  WWW.  U. l.U t  tw  )  .] 

IF  (DISTB^IE , KD^ .LT.  RMIN)  THEN 


THRES1 

THEN 


77 

70 


62 


END  IF 

APPLY  SPECIAL  RULE (WIDTH  DETERMINATION  RULE  2)  FOR  BOUNDARY. 

IF(  NSPCL  .EQ.  1  )  THEN 

IF(  KD  .EQ.  2  )  THEN 
RMIN  *  0. 

END  IF 

END  IF 
KD  *  KD  +  1 
END  IF 

NO 7  IE  5  9  KD  •  1 
CONTINUE 

ANGLE  OF  INCLINATION  IS  GREATER  THAN  135  OR  LESS  THAN  45  DEGREES. 
ELSE  IF  (  ANGLE  .GT.  135.  .OR.  ANGLE  .LT.  45.)  THEN 

COMPUTE  THE  SLOPE  ON  THE  LEFT  SIDE. 

DO  60  IE  =  1,K 
KD  =  1 

IF(  ANGLE  .GT.  135, 

IF(  ANGLE  .LT.  45. 

IH1  =  IH 
YB(IE ,KD)  »  IH 
IY  *  YA(IE) 

IX  >  XA(1E) 

IA  a  YB(IE.KD) 

XB(IE,KD)  a  XCORD(  YM,  IA,  IY,  IX  ) 

IF(  XB(IE,KD).  .LE.  0.)  THEN 
IH  =  IH1  +  200 
GO  TO  64 
END  IF 


)  IH  a 

(INT 

:ya< 

/100) 

*  100) 

-  100 

)  IH  a 

(INT 

[YA( 

!ie) 

/100) 

*  100) 

+  100 

J1 

J2 

RA 


a  XB  ( IE ,  KD 
a  YB(IE  ,KD 


/  100 
/  100 


IF  (  RA 


/  100. 


J3  =  J1  +  1 
ELSE  IF(  RA  .LT.  0.5  )  THEN 
J3  =  J1 
END  IF 

IF  (  XB(IE,KD)  .GT.  10000.  .OR.  YB(IE,KD)  .GE.  10000.  )  THEN 
KD  =  JQ  -1 
GO  TO  65 
IF 


END 


DISTB(IE,KD)  =  SQRT 

HTB  »  RA  *  (ALTD 
HT2(IE,KD)  =  ALTD 
DHT2(IE,KD)  =  HT2( 


:T((  XB(IE,KD)  - 
+  (  YB(IE,KD)  - 
{ J1 ,  J2)  -  ALTD( J 

(J1,J2)  -  - 

IE,KD)  - 


-  XA(IE 
YA(IE 
( Jl+1 , J2 

HTB 
HTflE) 


**2 

**2) 
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IF (  DISTB(IE.KD)  .EQ.  0.)  GO  TO  65 
SL0PE2(IE,KD)  =  DHT2(IE,KD)  /  DISTB(IE,KD) 

IF(  KD  .EQ.  1  )  THEN 
NCA  =  7 

ELSE  IF (  KD  .GT.  1  )  THEN 
JA  *  XB(IE,KD-1)  /  100 
JB  a  YB(IE,KD-1)  /  100 
NCA  ■  NCODE (  JA, JB) 

NCB  =  NCODE (J1,J4) 

IF(  ANGLE  .GT.  135.)  THEN 

IF  (  CODE(NCA,NCB)  .EQ.  1  .AND.  SL0PE2(IE,KD)  .LE.  THRES1 
*  .AND.  SLOPE2(IE,KD)  .GE.  THRES2  .Alb.  KD  .LE.  10)  THEN 

IH  »  IH  -  100 
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ELSE 


SL0PE2(XE,KD)  .GT.  THRES1 
.OR.  KD  .GT.  10)  THEN 


KD  ■  KD  ♦  1 
GO  TO  62 

XF  ( CODE (NCA ,NCB)  .NE.  1  .OR.  SL 
.OR.  SL0PE2(IE,KD)  .LT.  THRES2  . 

XH  *  XH1  +  200 
XF (DXSTB (IE , KD)  .LT.  LHXN)  THEN 
LMIN  ■  DISTB(IE.KD) 

END  IF 

APPLY  SPECIAL  RULE (WIDTH  DETERMINATION  RULE  2)  FOR  BOUNDARY. 

IF (  NSPCL  .EQ.  1  )  THEN 

IF(  KD  .EQ.  1  )  THEN 
LAIN  -  0. 

END  IF 

END  IF 
GO  TO  61 
END  IF 

ELSE  IF (  ANGLE  .LT.  45.)  THEN 

IF  (  CODE (NCA, NCB)  .EQ.  1  .AND.  SLOPE2 (IE , KD )  .LE.  THRES1 

*  .AND.  SL0PE2(IE,KD)  .GE.  THRES2  .Alb.  KD  .LE.  10)  THEN 

IH  »  IH  +  100 
KD  *  KD  +  1 
GO  TO  62 

ELSE  IF (CODE (NCA, NCB)  .NE.  1  .OR.  SL0PE2(IE ,KD)  .GT.  THRES1 

*  .OR.  SL0PE2 ( IE , KD )  .LT.  THRES2  .OR.  KD  .GT.  10)  THEN 

IH  *  IH1  -  200 

IF(DISTB(IE,KD)  .LT.  LMIN)  THEN 
LMIN  =  DISTB(IE.KD) 

END  IF 

APPLY  SPECIAL  RULE (WIDTH  DETERMINATION  RULE  2)  FOR  BOUNDARY. 


IF (  NSPCL  .EQ.  1  )  THEN 

IF(  KD  .EQ.  1  )  THEN 
LMIN  »  0. 

END  IF 


END 

END  IF 
GO  TO  61 
END  IF 
END  IF 

*  COMPUTE  THE  SLOPE  ON  THE  RIGHT  SIDE. 

61  KD  =  KD  +  1 

64  YB(IE,KD j  a  IH 

IY  »  YA(IE) 

IX  »  XA(IE) 

IA  »  YB(IE ,KD) 

XB(IE.KD)  =  XCORD(  YM,  IA,  IY,  IX) 

J1  =  XB(IE ,KD)  /  100 

J2  »  YB(IE,KD)  /  100 

RB  a  (MOD  (XB(IE,KD),  RY))  /  100. 

IF(  RB  .GE.  0.5  )  THEN 
J3  »  J1  +  1 

ELSE  IF(  RB  .LT.  0.5  )  THEN 
J3  a  J1 
END  IF 

IF(XB(IE,KD)  .GT.  10000.  .OR.  YB(IE,KD)  .GE.  10000.  )  THEN 
IF(DISTB(IE,KD-1)  .LT.  RMIN)  THEN 
RMIN  a  DISTB(IE,KD-1) 

END  IF 
GO  TO  67 
END  IF 

DISTB(IE,KD)  «  SQRT ( (  XB(IE,KD)  -  XA(IE)) 

*  +  (  YB(IE.KD)  -  YA(IE ) )  **2) 

HTB  »  RB  *  (ALTD( J1 , J2)  -  ALTD( Jl+1 , J2) ) 

HT2(IE,KD)  *  ALTD(J1,J2)  -  HTB 
DHT2(IE,KD)  »  HT2(IE,KD)  -  HT(IE) 

IF(  DISTB(IE.KD)  .EQ.  0.)  GO  TO  66 
SLOPE2(IE,KD)  a  DHT2(IE,KD)  /  DISTB(IE,KD) 

IF(  ANGLE  .GT.  135  .AND.  IH  .EQ.  XH1*200)  THEN 
NCA  «  7 

ELSE  IF(  ANGLE  .LT.  45  .AND.  IH  .EQ.  XH1-200)  THEN 
NCA  *  7 
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IF  ( 


NCB 


ELSE  XF(  XD 
JA  “ 
JB 
NCA 

END  IF 
NC0DE(J1,J4) 


CODE (NCA, NCB)  .10.  1  .AND.  SL0PE2 ( IE , KD)  .LE.  THRES1 
.AND.  SL0PS2(IE ,KD)  .GE.  THRES2  .AND.  KD  .LS.  20)  THEN 


* 

GT. 

IH  +  100 


.»■*:.  135.)  THEN 


IF(  ANGLE 
IH 

END  IF 

IF(  ANGLE  .LT.  45.)  THEN 
IH  *  IH  -  100 

END  IF 

IF (  IH  .LS.  0  )  GO  TO  60 
KD  *  KD  -  1 
GO  TO  64 

ELSE  IF  ( CODE (NCA.NCB )  .NE.  1  .OR.  SLOPE2(IE,KD)  .GT..  THRES1 
*  .OR.  SLOPE2(IE,KD)  .LT.  THRES2  .OR.  KD  .GT*  20)  THEN 

IF (DISTB (iE , KD )  .LT.  RMIN)  THEN 
RMIN  *  DISTB(IE,KD) 

END  IF 

APPLY  SPECIAL  RULE (WIDTH  DETERMINATION  RULE  2)  FOR  BOUNDARY. 

IF(  NSPCL  .EQ.  1  )  THEN 

IF(  KD  .EQ.  2  )  THEN 
RMIN  *  0. 

END  IF 


END  IF 
KD  *  KD  +  1 
END  IF 

67  NO  ( IE )  *  KD  -  1 

60  CONTINUE 

END  IF 

*  DETERMINE  THE  TOTAL  ARC  WIDTH. 

*  ASSIGN  THE  WIDTH  OF  ROAD  ITSELF. 

ROADW(l)  *  18. 

ROADW? 2 )  «  6. 

ROADW ( 3 )  ■  4. 

R0ADW(4)  *  4. 

ROADW (  5 )  »  2. 

ROADW(6)  «  5. 

ROADW(7)  »  4. 

IF  (ARCTP  .NE.  7)  THEN 

WIDTH  *  LMIN  +  RMIN  +  ROADW (ARCTP) 

ELSE  IF  (ARCTP  .EQ.  7)  THEN 
WIDTH  a  ROADW (ARCTP) 

END  IF 

RETURN 

END 

*********************************************************************** 


*  I.  SUBROUTINE  FOR  COMPUTING  COORDINATES  FOR  THE  MINIMUM  DISTANCE. 
*********************************************************************** 


SUBROUTINE  MIND(K,XM, YM,XA,YA, LMIN, RMIN, XL, YL\XR,YR) 

REAL  LMIN 

DIMENSION  XA(90),YA(90) 

DIMENSION  XL(10),YL(10),XR(10),YR(10) 

MK  a  1 

ANGLE  OF  INCLINATION  IS  BETWEEN  0  AND  89.9999  DEGREES. 
TAN(89.9999)  =  572957.7 

IF  (  KM  .GE.  0.  .AND.  XM  .LT.  572958.0)  THEN 
DO  80  J  a  i  R,  K-l 
INY  a  YA(J) 
in::  *  XA(j) 

YL(MK)  »  yA(J)  +  SQRT(  LMIN**2  /  ((l./YM  **2)+l.)) 
INA  «  YL(MK) 
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80 

A 


Ml  .10.  Q  }  00  TO  80 

■  YAfS?  -  S$RTt'8ttti®2  )  ((l./YH  **2)+l.)) 
_  YR(MK) 

IF  (  INI  .10.  0  )  GO  TO  80 
XR •  llT  YK,XNA,INY,IKX  ) 

CONTINUE 

ANGLE  OF  INCLINATION  IS  BETWEEN  90.0001  AND  179.99  DEGREES. 
ELSE  IF  <  XM  .LT.  0.)  THEN 
DO  81  J  ■  1.  K.  K-l 
INY  -  YA(J) 


I  NX  *  XA(J 
YL(HX) 


Yl(J)  -  SQRT(  LMIN**2  /  ((l./YM  **2)+l.)) 
0  ^00  TO  81 


81 

A 


INA  -  YL(MK) 

IF  (INA  .EO.  .  _ 

XL(HK)  *  XCORD(  for  INA. INY, INX  ) 

YR(MK)  *  YA( J)  ♦  SQRT(  RMIN**2  /  ((l./YH  **2)+l.)) 

INA  -YR(MX) 

IF  (INA  .SO.  0  )  GO  TO  81 
XR(MX)  »  XCORD(  YM, INA, INY, INX  ) 

MX  «  MX  +  1 
CONTINUE 

ANGLE  OF  INCLINATION  IS  90  DEGREES. 

ELSE  IF  (  XM  .GE.  e72958.0)  THEN 
DO  32  J  «  1.  K,  K-l 
INY  -  YA( J) 

INX  *  XA(J 
YL(MK)  *  YA(J) 

INA  *  YL(MK) 

IF  (  INA  .EQ.  0  )  GO  TO  82 
XL (MX)  *  INX  -  LMIN 
YR(MK)  *  YA( J) 

INA  -  YR(MK) 

IF  (_  INA  .EQ.  0  )  GO  TO  82 
XR(MK)  -  INX  +  RMIN 
MK  ■  MX  +  1 
82  CONTINUE 

END  -IF 
RETURN 
END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAkAA 

*  J.  SUBROUTINE  FOR  COMPUTING  THE  RATE  OF  FLOW. 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


SUBROUTINE  FLOW  ( SP ,  WIDTH ,  ARCTP ,  CASE  ,,FR ,  MAXNO  ) 

TNTRGRR  ARCTP  CASE 

DIMENSION  SP(7 . 2) , FR(5. 5) ,DSTN(6) ,MAXN0(3 ,3) 

DIMENSION  DW(3,2).DD(3,2j 

DATA  DSTN  /  25. ,5. ,20. ,5. ,5. ,3.  / 

*  DSTN  IS  DISTANCE  BETWEEN  EACH  ELEMENTS. 

I  *  ARCTP 

*  ASSIGN  DOCTRINAL  WIDTH  AND  DEPTH  OF  BATTALION. 


DW 

1,1 

*  1. 

DD 

1,1 

*  3. 

DW! 

1,2 

■  0.005 

DD 

1.2 

■  4. 

DW 

2,1 

a  1. 

DD 

2,1 

■  3. 

DW 

2,2 

■  0.005 

DD 

2,2 

»  4. 

DW 

3,1 

■  1. 

DD 

3,1 

a  3. 

DW 

3,2 

a  0.006 

DD 

3,2 

«  2.5 

WIDTH1 

a  WIDTH 
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COMPUTE  THE  FLOW  RATE .  (  BATTALION  /  HOUR  ) 
ZF  (  CASE  .EQ.  1  7  THE! 

DO  80  J  ■  1*2 
K  ■  1 


80 


CONTINUE 
DO  81 


D  ■  DW(J,K)  *  DD(J ,K) 

FR(J,K)  »  (  SP(Z;jrX)  *  WZDTHl  )  t  D 


J  •  1,2 
K  ■  2 

IF (WIDTH  .GE.  5.)  THEN 
WZDTHl  -  0.6' 

END  IF 


81 


0.005 

D-  DW(J,R)  *  DD(J,R> 
FR(J,R)  ■  (  SP(I,J,R) 

2  )  THEN 


*  WZDTHl  )  /  D 
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CONTINUE 

ELSE  IF  (  CASE  .EQ. 

j  •  3 
DO  33  K  *  1,2 

XF(K  .EQ.  2  .AND.  WIDTHl  .GE.  0.008)  THEN 
WlDTHl  -  0.006 

ELSE  IF(K  .EQ.  2  .AND.  WlDTHl  .LT.  0.006)  THEN 
WlDTHl  •  WlDTHl 
END  IF 

D  *  DW(J.K)  *  DD(J.K) 


END  IF 


FR 
CONTI 


i  «  DW( J,K) 
»K>  ■ S 


SP(I< 


*  WlDTHl  )  /  D 


*  COMPUTE  THE  MAX  NUMBER  OF  ELEMENT  IN  MULTIPLE  COLUMN  FORMATION. 

IF  (  CASE  .EQ.  1  )  THEN 
KNA  -1 

DO  85  RA  *  1,2 
DO  85  RB  *  1,2 

IF(  RB  .EQ.  1)  THEN 

MAXNO(KA.KB)  *  WIDTH  /  DSTN (KNA) 

ELSE  IF(KB  .EQ.  2)  THEN 

MAXNO(KA,XB)  »  (WlDTHl  *  1100.)  /  DSTN(KNA) 

END  IF 

KNA  *  KNA  +  1 

85  CONTINUE 

ELSE  IF^CASE  .EQ.  2  )  THEN 

KA  *  3 

DO  86  KB  *  1,2 

IF(  KB  .EQ.  1)  THEN 

MAXNO(KA.KB)  *  WIDTH  /  DSTN (KNA) 

ELSE  IF(KB  .EQ.  2)  THEN 

MAXNO(KA.KB)  »  (WlDTHl  *  1000)  /  DSTN (KNA) 

END  IF 

KNA  *  KNA  +  1 

86  CONTINUE 
END  IF 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

*  K.  SUBROUTINE  FOR  COMPUTING  THE  MIN  TRAVERSAL  TIME  ALONG  THE  ARC. 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  TIME (ARCT? , SP , UNITF , CASE , DIST , WIDTH , ATIME ) 

INTEGER  ARCTP, CASE, UNITF 
REAL  KR 

DIMENSION  ATIME(3) ,SP(7 ,3,2) 

DATA  KR  /  0.4/ 

*  KR  IS  COEFFICIENT  OF  ROUTE  AVAILABILITY. 


I 


»  ARCTP 


90 


DXST_  /  1000. 


95 


* 

* 


97 


HTF  .EQ._1  >lI^ 


ADXST  • 

ir  (  uni  ,  .  % 

XF(  CASRKQ.  1  )  THEN 
DO  90  KU  -1,2, 

.  klUaXbi)  -  ADI  ST  /(  RR  *  SP  (  X ,  KU ,  1 )  ) 

CONTINUE 

ELSE  XF(  OSS  .10.  2  )  THEN 

ATXMI(3)  ■  AD  I  ST  /  SP(I,3,1) 

END  XF 

ELSE  XF  (  UNXTF  .EQ.  2  )  THEN 
XF(  CASE  .EQ.  1  )  THEN 
K  *  2 

DO  95  J  »  1,2 

ATIME (J)  ■  ADXST  /  SP(I,J,K) 

CONTINUE 

ELSE  IF(  CASE  .SO.  2  )  THEN 

ATXNE(3)  «  ADXST  /  SP(I,3,2) 

END  XF 
END  XF 

IF  ARC  WIDTH  IS  LESS  THAN  MIN  FORMATION  WIDTH, 

HAKE  THE  TIMS  BIO  NUMBER. 

BIG  *  9.99 
FORWV  »  4. 

FORWD  «  2. 

IF  (  UNXTF  .EQ.  1  )  THEN 

IF(  CASE  .EQ.  I  )  THEN 

IF(  WIDTH  .tT.  FORWV  )  THEN 
DO  97  KN  ■  1,2 
ATIME (KN>  -  BIG 
CONTINUE 
END  IF 

ELSE  IF?  CASE  .EQ.  2  }  THEN 

r(  WIDTH  .LT.  FORWD  )  THEN 


IF(  WIDTH  .LT.  FOR* 

ATXME(3)  *  BIG 
END  IF 
END  XF 
END  IF 
RETURN 

*********************************************************************** 

*  L.  SUBROUTINE  FOR  PRINTING  THE  RESULTS. 
*********************************************************************** 


SUBROUTINE  PRT 


INTEGER 

PITAT 

DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
COMMON  / 


ALTD , ARCTP , CASE , CODE , PRINT , UTMX , UTMY ,  UNITF , SPD 
LMIN 

ALTD (100 , 100 ) , NCODE ( 100 . 100 ) , NX( 10 ) ,NCODEX( 10 ) 
XNODE(90) , YN0DE(90) ,CODE(0 :7 ,0  *7) 

DISTA(90  ,XA(90l,YAi90) 

- J - 90  n  ,  SLOPE  1(9( 


HT ( 90 ) . DHT ( 90 ) , SLOPE1 ( 90 ) 

DISTB ( 100 , 100 ) , XB(100 , 100 ) , YB ( 100 , 100 ) 

HT2 ( 100 . 100 ) , DHT2 (100,100), SLOPE 2 ( 100 , 100 ) 

NO ( 90 ) , SLPCH ( 90 ) 

XL(10) , YL(10) ,XR(10) , YR(1Q) 

SP(7  3,2),FR(5,5),MAXNC(3,3) 

ACTD ( 90 } , ATIME ( 3 ) 

R0ADW(7) 

INTS  /  NSTEP , 1X1 , 171 , 1X2 , 172 , JA , JB . K , NARC , NCOUNT , 

NODEA , NODEB , Ji , J2 , J3 , J4 , PRINT , UTMX , UTMY , 

CASE .ARCTP, FR , AD , UNITF , NSPCL 
COMMON  /  INTA  /  NCODE , NO , ALTD , CODE , MAXNO 

COMMON  /  SREAL  /  XM, ANGLE , DIST , TSLPCH , ASLPCH , ALT1 , ALT2 , AVGSLP , 

LMIN , RMIN , XL , YL , XR , YR.YM , RY , WIDTH , THRES 1 , THRES2 
COMMON  /  AREAL  /  XA , YA ,HT , DRf , DXSTA, DISTB .SLPCH, SLOPE 1 , SLOPE2 , 

XB , YB , HT2 , DHT2 , SP , ACTD , ATIME , RbADW 
OPEN  (  UNIT  ■  10,  FILE  «  'NETWl '  ‘ 

OPEN  (  UNIT  ■  20,  FILE  ■  'NKTW2' 


l!  1 
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OPEN  (  UNIT  -  25,  FILE 


•META'  ) 
•HETB'  ) 


700 


701 

730 


731 

91 


735 

93 

733 

740 

741 

742 
750 


751 
94 

752 


753 

754 
756 


tTXON  FOR  EACH  ARC. 
THEN 

.NODEA, NODES 


OPEN  (  UNIT  *  26;  FILE 

RINT  GENERAL  INFORMA' 

IF  (PRINT  .10.  1  T 

WRITE (6 .700)  NCOUNT,.. _ 

A  FORMAT?*', /.7X,'  kkc  .  *12. 

*  /7,4X,71(,-'1./>,51,  'NODE  '  ,  il .  $X,  'NODE  M2, 

*  SX,'  XM  ANGLE  DISTANCE  ACT.  DISTANci ' , 

■**  f  ,4X, 1  1X1  IY1' , 5X, '  1X2  IY2') 

WRITS* 6 .701)  1X1 , IY1 . 1X2 , IY2.XM, ANGLE , DI5T . AD 
FORMAT^  '^iX,2I6.2X,2I6,lX,F3.2,F8,2,2X,Fil.2,2X,F10. 


2) 


WRITE  (( 

FORMAT* *  ' ,/,4X,71( '-'),//, 20X, 'SOME  POINTS  ALONG  ARC'/, 

*  4X, 'POSITION  ' ,7X, 'COORDINATE*  X  :  ?  ) 1 ,5X, 'ALTITUDE' ,j\ 

DO  91  MA  *  1,X 

WRITE (6 , 731 )  MA.XA(MA) ,YA*MA) ,HT*MA) 

FORMAT*'  \3X, 'POINT  ',13 ,8X, ' ( ' ,2F8.1 ,1X, ' ) • ,  F13.1) 
CONTINUE 
WRITE (6, 732) 

732  FORMAT?1  ',/, 4X, 71* , 

*  fOSITION'Jx  1  DELTA  H*  ,5X,  'DISTANCE*  ,6X,  'SLOPS'  ,5X. 

*  'ACTUAL  DISTANCE') 

DO  92  MB  »  l.K-1 

WRITE  (6,733)  MB , MB* 1 , DHT (MB ) , DISTA(MB ) ,  SLOPEl  *  MB ) ,  ACTD  (MB ) 
FORMAT*'  ',3X,  'POINT  '.13,'  -',13,  2F12.1,  F12.3.H2.1) 
CONTINUE 
WRITE(6,734) 

FORMAT* '  1 ,/,4X. 71* '-'),/. 

*  7X, 'POSITION  ' ,5X, '  SLOPE  CHANGE') 


733 
92 

734 


DO  93  MC 


l.K-2 

WRITE *6, 735)  MC.MC+l .SLPCH(MC) 

FORMAT* »  ',ix,  'POINT  M3,'  -',I3,2X,  FU.2) 
CONTINUE 

WRITE (6, 738)  TSLPCH.ASLPCH 
v  FORMAT* 1  '.  /,4X, 'TOTAL  SLOPE  CHANGE  s  '  F1Q.4, 

*  /,4X  'AVERAGE  SLOPE  CHG.  :  \F10.4,/ ,4X,*1< ) 

WRITE (6. 740)  NODEA.NODEB 

FORMAT* (  ' ./ ,4X, 'NODE' ,13. ' »  ALTITUDE  1  ' , 

*  t  v  4X, 'NODE  , l3 , '  i  ALTITUDE  2  {) 

WRITE (6, 741)  ALT! , ALT2 


WRITE (6 .741)  ALT1.ALT2 
FORMAT?*  '  20X.F5.1.23X.F5.1) 
WRITE  *  6 . 74$ )  AVGSLP 
format 7 1  '  /,4X, 'AVERAGE  SLOPE 


' , F7. 4, /,4X, 71 ('-'),/  ) 


FORMAT*1 
WRITE (6, 750) 

FORMAT*'  ' ,15X, '  SOME  POINTS  OFF  ARC  *  FROM  EACH  POINT  ' 

*  'ON  THE  ARC  )',/ ,7X, ‘COORDINATE  (  X  i  Y  )', 3X ,' ALTITUDE ' , 

*  4X. 'DELTA  H' ,4X, 'DISTANCE' ,5X, 'SLOPE* ,/) 


Ml  »  0 

DO  94  JX  *  l.K 
N1  *  N1+  1 

DO  94  JY  =  1 ,N0(N1 ) 

WRITE (6, 7 51)  JX,JY,XB*JX,JY),YB*JX,JY) ,HT2(JX, JY) , 

*  DHT2 ( JX , JY ) , DISTB ( JX , JY ) , SLOPE2  *  JX , JY ) 

FORMAT*'  ',3X. 'POINT  ',12,'  -',12,'  * ' , 2F7 . 1 , IX. ' ) ' , 4F11 .2 ) 
CONTINUE 

WRITE (6, 7 52)  LMIN.RMIN, WIDTH 

FORMAT*1  ' ,/,4X, 71* '-),/, 25X, 'MINIMUM  DISTANCE  OFF  ARC', 

*  /,4X,'LEFT  MIN  DISTANCE  j'.FS.l, 

*  /,4X, 'RIGHT  MIN  DISTANCE  j',F8.1, 

*  /,4X  'ARC  WIDTH  «',F8.lJ 

WRITE *6 ,753) 

.  FORMAT  * '  ' ,/,27X, 'POINT  1'13X, 'POINT  2', 

*  /,24X ,'X  LEFT  Y  LEFT' ,4X, 'X  RIGHT  Y  RIGHT') 

WRITE(6 .754) (XL*K) . YL?K) ,XR*K) , YR(K) ,  K=l,2) 

^  FORMAT*'  '  ,/,4X, 'INITIAL  POINT  i  >,4F10.1, 

*  ,  / , 4X , ' END  POINT  »  ',4F10.1) 

WRITE *6, 7 56; 

FORMAT*1  ' ,/,4X,71( '*' )) 
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PRINT  FLOW  RATI  AND  7HA1 _ 

I LSI  IF  (  PRINT  .EQ.  2 


no 


Tim. 
THIN 


SSHP;W  HP1K "‘“'7°%^ 

ffi.ft S*?mi8 ! Loli , out 


W.WIDTH.ARCTP 
CNINuN  DISTANCE  OFF 


766 

9S 


768 

96 


7620 


7660 

950 


7680 

763 

764 

765 
767 
* 


*//.4X.'ARC  WIDTH  .{.F8.1, 

J  /.4X.ARC  TYPE  *' Jo, 

*  / j4X , ' UNIT  TYPE* ,4X, 'MISSION  WPi' ,3X, 

4  i  4*  au  mm  — i—  i  — -  -  '  9 


#  w«4*  44«8  »  ai 

FLOW  RATE'  3X, 'MAX  *  ELEMENT') 

IF  ^CASE  .EJ.^IJ  THEN 

DO  95  IK  * 

WRITE (20 , 766 )  I J , IK , FR ( IvJ , IK) , MAXNO( IJ,IK) 
^FORMAT?'  1 , 7X,X2,9X,I2, 10X,F? !2,9X,I4) 
CONTINUE 

ELSE  IF  (  CASS  .EQ.  2)  THEN 
IJ  »  3 

DO  96  IK  *  1.2 

NRITE ( 20 , 768 )  I J , IK , FR( I J , IK) , HAXNO( I J , IK ) 
FORMAT (•  ' ,7X,X2,9X,X2,10X,Fo.l,9X,I4) 
CONTINUE 
END  IF 

WRITE? 20 , 7620 } 

FORMAT ('  '  7,4X, 'UNIT  TYPE  MIN  TRAVERSAL  TINE') 

IF  (  CASE  .EO.  1)  THEN 
DO  950  K  *  1.2 

WRITE (20. 7^60)  K.ATIHE(K) 

FORMAT^  ',7X.i1.9x,f*.2) 

CONTINUE 

ELSE  IF  jj^CASE  .EQ.  2)  THEN 

WRITE (20 .7680)  IK.ATIME(IK) 

_  _  FORMAT ( '  ' ,7X,I2,$X,F9.2) 

END  IF 

WRITE (20, 763) 

FORMAT( J* ./»27X, 'POINT  1*  13X. 'POINT  2', 

A24X,'X  LEFT  Y  LEFT*  ,4X.  'X  RIGHT  Y  RIGHT') 
wftiTE?fc>.764)  XL(1) ,YL(1)!xr(i) ,YR(1)  ' 

FORMAT?'  \ 7, 4X, 'INITIAL  $OINT  «  ‘  .4F10.1) 

WRITE (20, 7WJ  X£(2),YL(2),XR(2),YR(2) 

FORMAT?'  •  3X,'ENd'  POINT  t(,4Fl6.1) 

WRITE (20, 767) 

FORMAT ( '  ' ,/,4X,71( '*' )) 


PRINT  COORDINATES  OF  BOUNDARY  LINE. 
ELSE  IF  (  PRINT  .EQ.  4  )  THEN 
DO  97  IN  *  1,2 


97 

774 

775 
4 


XL(IN) 

YL(IN 

XRiIN) 

YRfIN 


CONTINUE 


UTMX  +(XL 
UTMY  +(YL 
UTMX  +(XR 
UTMY  +  (YR 


IN 

IN 

IN 

IN 


/ 1000. 
/ 1000. 
/ 1000. 
/1000. 


PRINT  NETWORK  DATA  FOR  CARTESIAN  SPACE  NETWORK  PROGRAM. 
ELSE  IF  (  PRINT  .EQ.  5  )  THEN 

IF(  NODEA  .GT.  NODES  )  THEN 
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ARC', 


874 


H0DI1  •  NOOtB 
HODE2  «  HOOKA 

KLSI  IP7  M0DI8  ,GT.  NODKA  )  THIN 
HODK1  •  HOOKA 
NODI2  •  NOOtB 
IND  XP 

DST  ■  OZST  /  1000. 

NO  •  WIDTH  f  1000. 

XP(  CASK  .IQ.  1  )  THIN 

IS  ■  1 

KLSI  XP ("cASI  .IQ.  2  )  THIN 
XS  «  3 
IT  ■  1 
END  XP 

NT INI  -  ATXHt(RS)  *  100 

NPR  •  FR(KS,RT)  *  100 

XARCT  ■  ARCTP 
I UNITE  -  UNXTP 

S«  SP(  XARCT,  KS,  I  UNITE  ) 

CASK  .10.  1  j  fKKN 

,874)  NCOUNT ,  NODEl , N0DE2 . NTXMB , NFR , DST , HD , SPD 


2  )  THEN 

IDI1 .N0DK2 ,NTXMI,NPR,OST,NO,SPD 


nwwm , 

ELSE  IF(  CASE  .E 
WRITE (26. 874)  MCOUNT, 

FORMATS  * , 1X,I3,2IS,2X6,3X,F4.1,3X,F6.3,2X,X4) 

END  IP 

RETURN 

END 
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APPENDIX  D 

COMPUTER  PROGRAM  FOR  CARTESIAN  SPACE  NETWORK 

This  appendix  contains  the  computer  program  used  to  support  the  algorithms  in 
Chapter  IV. 


PROGRAM  HIT 

JULY  13,  1987  (13(00) 

NETWORK  MAIN  PROGRAM 
SHORTEST  PATH  ALGORITHM 

*************  VARIABLE  DECLARATION  ************* 


INTEGER 
INTEGER  CASS , NTYPt , ?A 
CHARACTER  *1  RISK 
DIMENSION  TAIL7i50).HEA 
DIMENSION  FLOW (100) ,DIS 
DIMENSION  LINK? 100,100) 
DIMENSION  NO(lOO) ,SP(10 
DIMENSION  INODE ( ISO ) , JN 
DIMENSION  ISPATH(SO) 
COMMON  /  INTS  /  CASE. 


TAIL, HEAD, TIME, FLOW 

CASE , NTYPE , PARAM . S? , UNITE , ROWS? 

*1  rAsm 

TAIL7iS0) ,HEAD(150) .TIME (150) 
FLOW (100) ,DIST(100) .WIDTH(IOO) 


ARCOST(lSO) 


COMMON  /  INTS  /  CASE. NTYPE, N,M. MARC, ISTART, LAST, PARAM. 

NUMP . UNITF , NVEH , NCOL , NROW , ROWS? .MINE , 
NWD , NDP , NFORWD 

COMMON  /  INTA  /  TAIL, HEAD, TIME. FLOW, LINK, LARC, NO, SP, 
INODE , JNODE . ISPATH 

COMMON  /  REALS/  BIG . TPATH , ELEN , PATHL , FORWD , FORD? , PRT 
COMMON  /  REALA/  WIDTH, DIST.RBSM 


*  (1).  SET  THE  INITIAL  OPTIONS. 

CALL  OPTION  (CASE, UNITF, NCOL, KOWSP, PARAM, FORWD  .NFORWD , 

*  RE SM, MINE) 

*  SET  THE  START  NODE 

ISTART  «  1 

*  .  SET  THE  LAST  NODE 

LAST  «  34 

*  :  TOTAL  NUMBER  OF  ARC  IN  THE  SECTOR 

NARC  »  72 

TOTAL  NUMBER  OF  NODE  IN  THE  SECTOR 

NNODE  »  34 

*  SET  THE  INITIAL  CONDITIONS. 

CALL  XNIT  (ARCOST, NUMP, ISPATH, XLEN, ADJUST, PATHL, SPEED. 

*  LINK, TPATH. NNODE, BIG) 


*  (2).  READ  THE  DATA  FROM  THE  DATA  FILE. (DATA  1.  OR  DATA  2  ) 

*  DATA  1  ;  ARCS  AND  CHARACTERISTICS  FOR  VEHICLE  UNIT 

*  DATA  2  ;  ARCS  AND  CHARACTERISTICS  FOR  DISMOUNTED  TROOPS 

*  DETERMINE  VARIABLE  FOR  READING  INPUT  DATA. (  VEHICLE  OR  DISMOUNTED  ) 

IF(  PARAM  ,EO.  1  )  THEN 

IF(  CASE  .EQ.  1  )  THEN 
KUNIT  ■  1 

ELSE  IF(  CASE  .EQ.  2  )  THEN 
KUNIT  *  2 
END  IF 
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100 

* 

999 

* 


ELSE  IF(  PARAM  .EQ.  2  )  THEN 
KUNIT  *  1 
END  IF 


READ (KUNIT, 100,  END=999) 
FORMAT (2X 


T,100,  END=999)  (NO(I  ,TAIL(I) .HEAD(I) ,TIME(I) . 

7LOW(I)  DiST(l)  WI.  .A(fJ,SP(I)  i=l,NAR65 
,I3,2I5,2I6,3X,F4.1,3X,F6.3,2x'f4) 


SPECIFY  THE  WIDTH  AND  DEPTH  OF  UNIT  FORMATION. 

CALL  FORMTN  (NARC, CASE, UNITF , NCOL , ROWSP , PARAM, WIDTH, TIME, 
*  NROW , DIST , FLOW , FORWD , FORDP , NVEH , NWD , NDP ) 


J  (3).  DETERMINE  THE  APPROPRIATE  VALUE  FOR  PARAMETER. 


*  DETERMINE  THE  VALUE  OF  PARAMETER  WHEN  THE  MINEFIELD  IS  ON  ARC. 

IF  (RESM  .EO.  'Y')  THEN 

IF  (PARAM  .EO.  1)  THEN 

TIME (MINE)  =  TIME (MINE)*  4 
END  IF 
END  IF 

*  DETERMINE  THE  APPROPRIATE  VALUE  FOR  PARAMETER 

N  =  NNODE 
M  =  NARC 

DO  10-  KA  =  1 , NARC 

INODE (KA)  =  TAIL(KA) 

JNODE(KA)  =  HEAD(KA) 

IF  (PARAM  .EQ.  1)  THEN 

•'ST(KA)  =  TIME(KA) 

ELSE  **  (PARAM  .EQ.  2)  THEN 
ARCOST(KA)  =  DIST(KA) 

END  IF 

10  CONTINUE 


J  (4).  DETERMINE  THE  MIN  (PARAMETER)  PATH. 


CALL  SHORTP  (N , M , INODE , JNODE , ARCOST , ISTART , LAST , BIG , NUMP , 

*  I S  PATH ,  XLF.N ,  NP  ) 

*  DETERMINE  THE  ARC  NUMBER  ALONG  THE  PATH. 

DO  60  M3  =  1,  NARC 

LINK (TAIL (M3), HEAD (M3))  =  N0(M3) 

60  CONTINUE 

PATHL  =  0. 

DO  70  MA  =  1,  NUMP-1 

LARC(MA)  =  LINK(ISPATH(MA) ,  ISPATH(MA+1 ) ) 

PATHL  -  PATHL  A  DIST(  LARC(MA)  ) 

70  CONTINUE 

*  WHEN  'PRT'  IS  1-0,  THERE  IS  NO  FEASIBLE  ROUTE  TO  GO. 

IF  (PARAM  .EQ.  1)  THEN 
PRT  =  0.0 

DO  75  MB  =  1 , NUMP-1 

IF  (TIME(LARC(MB ) )  .GE.  999)  THEN 
PRT  =1.0 
END  IF 

75  C.TINUE 

END  It 

*  DETERMINE  THE  TOTAL  TRAVERSAL  TIME  FROM  START  NODE  TO  LAST  NODE. 

SPEED  =  0. 

DO  80  M5  =  1,  NUMP-1 

SPEED  =  SPEED  +  (DIST(M5)  *  SP(M5)  /  PATHL  ) 

SC  CONTINUE 

ADJUST  =  FORDP  /  SPEED 

1PATH  =  XLEN  +  (  ADJUST  *  100.  ) 
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:<*>■ 


PRINT  THE  RESULT  OF  CALCULATION. 


a 


CALL  PRNT 
STOP 
END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

*  A.  SUBROUTINE  FOR  DETERMINING  THE  INITIAL  OPTIONS. 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  OPTION  ( CASE, UNITF.NCOL,ROWSP, PARAM, FORWD,NFORWD, 

*  RESM,  MINE) 

INTEGER  CASE , PARAM , UNITF , ROWS? 

CHARACTER  *1  RESM 


603 

IS 


16 

601 


17 

602 


WRITE (6 ,603) 
FORMAT ( '  ',///) 
PRINT* , 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT* , 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 


WHICH  TYPE  OF  PARAMETER  DO  YOU  DESIRE?1 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA • 

A  A  • 

*  1.  MIN  TIME  PATH  *' 

A  A  i 

*  2.  MIN  DISTANCE  PATH  *' 

A  A  i 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA i 

READ ( 5 i * )  PARAM 

PRINT*,  '**  NOTE  *  FOR  YOUR  REFERENCE,  CURRENT  ANSWER  IS', 
PARAM 

IF  (  PARAM  .EQ.  1  )  THEN 
GO  TO  16 

ELSE  IF  (  PARAM  .EQ.  2  )  THEN 
GO  TO  20 

ELSE  IF  (_PARAM  .NE.  1  .OR.  PARAM  .NE.  2  )  THEN 
PRINT*,  1  ' 

PRINT*,  '***  ERROR  *  5NTER  THE  NUMBER  1  OR  2  ***' 
PRINT*,  '  ' 

GO  TO  15 
END  IF 
WRITE (6, 601) 

FORMAT (  ',////////) 

WHICH  TYPE  OF  UNIT  DO  YOU  DESIRE?  (  1  OR  2  ) ' 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA • 

A  A  i 

*  1.  VEHICLE  UNIT  *' 

A  A  I 

*  2.  DISMOUNTED  TROOPS  *' 

A  A  i 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA ■ 


FOR  YOUR  REFERENCE,  CURRENT  ANSWER  IS1 


PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

PRINT*, 

READ (5,*)  CASE 
PRINT*,  '**  NOTE 
CASE 

IF  (  CASE  .EQ.  '1  .OR.  CASE  .EQ.  2  )  THEN 
GO  TO  17 

ELSE  IF  (  CASE  .NE.  1  .OR.  CASE  .NE.  2  )  THEN 
PRINT*,  '  ' 

PRINT*,  '***  ERROR  s  ENTER  THE  NUMBER  1  OR  2  ***' 
PRINT*,  '  ' 

GO  TO  16 
END  IF 
WRITE(6 . 602) 

FORMAT? !  ',////////) 

PRINT*,  '  WHICH  TYPE  OF  FORMATION  DO  YOU  DESIRE? 

PRINT*,  '  ' 

PRINT*  '  AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA I 

PRINTS  '  *  *  ' 
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*  2.  SINGLE  COLUMN  FORMATION  *' 

*  * 
************************************1 


PRINT*  1  *  1.  MULTIPLE  COLUMN  FORMATION  *' 

PRINT*,  »  *  *  • 

PRINT*  •  *  2.  SINGLE  COLUMN  FORMATION  *' 

PRINT*,  '  *  *  ' 

PRINT*  *  ************************************ i 

READ? 5'*)  UNITF 

PRINT*,  ’**  NOTE  «  FOR  YOUR  REFERENCE,  CURRENT  ANSWER  IS', 
UNITF 

IF  (  UNITF  .EQ.  1  .OR.  UNITF  .EQ.  2  )  THEN 
GO  TO  18 

ELSE  IF  (  UNITF  .NE.  1  .OR.  UNITF  .NE.  2  )  THEN 
PRINT*,  1  ' 

PRINT*,  '***  ERROR  s  ENTER  THE  NUMBER  1  OR  2  ***' 

PRINT*,  1  ' 

GO  TO  17 
END  IF 


IF  (CASE  .EQ.  1)  THEN 

IF  (UNITF  .EQ.  1) 


WRITE (6. 604) 


THEN 


FORMAT ( '  ',/////////) 

PRINT*,  '  ENTER  THE  NUMBER  OF  COLUMN. (5,  10,  20)' 

PRINT*,  '  ' 

PRINT*  '  *******  1 

PRINT*'  '  ' 

PRINT*  1  *******  1 

PRINT*'  1  1 

PRINT*  1  *******  1 

READ(5>)  NCOL 
PRINT*,  1  1 

WRITE (6 , 605) 

FORMAT ( 1  ',/////////) 

PRINT*,  1  ENTER  THE  DISTANCE  BETWEEN  EACH  ROW. (75,50,25) 

PRINT*,  '  ' 


PRINT*  1  /*****  **■ 

PRINT*,  '||  ' 

PRINT*  1  >******  *1 

PRINT*,  1  ' 

PRINT*  1  *******  1 

PRINT*,  '  ' 

READ (5,*)  ROWSP 

ELSE  IF  (UNITF  .EQ.  2)  THEN 
WRITE (6, 606) 

FORMAT ( 1  ',/////////) 

PRINT*,  '  ENTER  THE  DISTANCE  BETWEEN  EACH  ROW.  (75, 5"i, 25) 

PRINT*,  '  ' 

PRINT*,  1  /*******. 

PRINT*,  '||  ' 

PRINT*,  1  /***-«***' 

PRINT* ,  '  ' 

PRINT*  1  *******  1 

PRINT*,  1  ' 

READ (5,*)  ROWSP 
END  IF 

ELSE  IF  (CASE  .EQ.  2)  THEN 

IF  (UNITF  .EQ.  1)  THEN 
WRITE(6 . ol4) 

FORMAT (  ',/////////) 

PRINT*,  1  ENTER  THE  NUMBER  OF  COLUMN.  10,  20,  50)' 

PRINT*,  '  ' 

PRINT*  '  ******** 

PRINT*,  '  ' 

PRINT*  '  *******  1 

PRINT*,  '  ' 

PRINT*,  '  ******** 

READ? 5 , *)  NCOL 
PRINT*,  1  ' 

WRITE (6, 615) 

FORMAT(  ',/////////) 

PRINT*  '  ENTER  THE  DISTANCE  BETWEEN  EACH  ROW. (5, 10, 15) ' 

PRINT*,  '  ' 


/////////) 

ENTER 


THE  DISTANCE  BETWEEN  EACH  ROW. (5, 10, 15) 1 


PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
READ(5,*) 
END 


(j 


ROWSP 

IF 


END  IF 

20  WRITE (6. 607) 

607  FORMAT ( '  ',/////) 

PRINT*,  '  IS  THERE  ANY  MINEFIELD  ALONG  THE  ARC?  (Y  OR  N) 

READ(5 , 1 (Al) ' )  RESM 

IF  (RESM  .EQ.  ’Y1)  THEN 

print*,  •  ’  Inter  the  arc  number  of  minefield  arc 

READ (5,*)  MINE 

ELSE  IF  (RESM  .EQ.  'N* )  THEN 
MINE  =  0 
END  IF 
RETURN 
END 

★★art*********************************************;********************** 


i 


*  B.  SUBROUTINE  FOR  DETERMINING  THE  INITIAL  CONDITIONS. 

************ A********************************************************** 

SUBROUTINE  INIT ( ARCOST , NUMP , ISPATH , XLEN , ADJUST , PATHL , SPEED , 

*  LINK, TPATH,NNODE, BIG) 

DIMENSION  ISPATH (50) .ARCOST(IQO) 

DIMENSION  LINK (100, 100) 

BIG  =  1.0E10 
NUMP  *  0 
XLEN  *  0. 

ADJUST  =0. 

PATHL  »  0. 

SPEED  =  0. 

TPATH  =  0. 

DO  10  I  =  1 ,NN0DE 

ISPATH(I)  =  0 
10  CONTINUE 

DO  20  J  =  1,100 

ARCOST (J)  =  0. 

20  CONTINUE 

DO  50  Ml  *  1,  100 

DO  50  M2  =  1,  100 

LINK (Ml, M2)  =0 
50  CONTINUE 

RETURN 
END 

*********************************************************************** 

*  C.  SUBROUTINE  FOR  DETERMINING  THE  FORMATION  WIDTH  AND  DEPTH. 
*********************************************************************** 


SUBROUTINE  FORMTN  (NARC, CASE ,UNITF,NCOL, ROWSP, PARAM, WIDTH, TIME, 
*  NROW , D I ST , FLOW , FORWD , FORDP , NVEH ,  NWD , NDP ) 

INTEGER  TIME . FLOW 

INTEGER  CASE, PARAM, UNITF. ROWSP 

DIMENSION  TIME ( 100 ) , FLOW  ( 100 ) , DIST ( 100 ) , WIDTH (100 ) 

BIG1  =  999. 


J  (A).  DETERMINE  THE  FORMATION  WIDTH  OF  A  UNIT  OF  MOVEMENT. 
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i 

s 


1 


IF(  CASS  .EQ.  1  )  THEN 
NVEH  ■  50 

IF  (  UNITF  .EQ.  1  )  THEN 
NROW  ■  NVEH  /  NCOL 
COLSP  -  25.0 

FORWD  ■  7 (NCOL- 1)  *  COLSP)  /  1000 
ELSE  IF(  UNITF  .EQ.  2  )  THEN 
NROW  =  50 

NCOL  =  NVEH  /  NROW 
COLSP  «  4.0 

FORWD  =  (NCOL  *  COLSP)  /  1000 
'  ENp  IF 

FORDP  =  ( (NROW- 1 )  *  ROWSP )  /  1000 
ELSE  IF(  CASE  .EQ.  2  )  THEN 

WHEN  UNIT  IS  MULTIPLE  COLUMN  FORMATION,  THE  FORMATION  WIDTH 
WAS  ENTERED. 

IF  (  UNITF  .EQ.  1  )  THEN 
NMEN  *  1000 
NROW  «  NMEN  /  NCOL 
COLSP  a  S • 0 

FORWD  =  (NCOL  *  COLSP)  /  1000. 

FORDP  =  (NROW  *  ROWSP)  /  1000. 

ELSE  IF  (  UNITF  .EQ.  2  )  THEN 
FORWD  =0.006 
FORDP  =  2.5 
END- IF 

NWD  =  FORWD  *  1000 
NDP  =  FORDP  *  1000 

END  IF 


*  (B) .  ASSIGN  BIG  NUMBER  WHEN  ARC  WIDTH  IS  LESS  THAN  FORMATION  WIDTH. 


DO  50  1=1,  NARC 

IF  (PARAM  .EQ.  1)  THEN 

IF  (WIDTH(I)  .LT.  FORWD)  THEN 
.  TIME (I)  =  BIG1 
END  IF 
END  IF 

50  CONTINUE 
RETURN 
END 

*********************************************************************** 
*  D.  SUBROUTINE  FOR  DETERMINING  THE  SHORTEST  PARAMETER  PATH. 

ititititititititkititit'itititititititititkitititkit'kitkitititkitkitititititititititkkkititicititicitkititititititititititititititit 


* 

k 


it 
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SUBROUTINE  SHORTP  <N, M, INODE, JNODE , ARCOST , ISTART, LAST, BIG, NUMP. 
*  ISPATH , XLEN , NP  ) 


MINIMUM  PATH 

FIND  A  SHORTEST  PATH  BETWEEN  TWO  GIVEN  NODES 

INTEGER  INODE (M), JNODE (M) ,  ISPATH(N). 

REAL  WK4, ARCOST (M) 

LOGICAL  IWORK1 ( 34 ) , IFIN 

DIMENSION  IWORK2 ( 34 ) ,IWORK3(34) ,WK4(34) 

34  «  NUMBER  OF  NODES  IN  NETWORK. 


DO  10  I  *  1,N 
WK4(I )  =  BIG 
I WORK1 (I)  =  .TRUE. 
IWORK2 ( I )  =  0 
CONTINUE 

WK4( ISTART)  =  0. 

I  ■  ISTART 


I WORK II START)  -  .FALSE 
NP  ■  0 
XLEN  *  0. 


* 

« 

*  (A)  FOR  EACH  FORWARD  ARC  ORIGINATING  AT  NODE  I  CALCULATE  THE  LENGTH 

*  OF  THE  PATH  TO  NODE  I. 

* 

20  IC  -  0 

DO  30  K  ■  1,M 

IF  (INODE (K)  .EQ.  I)  THEN 
IC  ■  IC  +  I 
IW0RK3(IC)  *  K 
ISPATH(IC)  =  JNODE(K) 

END  IF 

IF  £JNODE^K)+.EQ.  I)  THEN 

IWORK3(IC)  =  K 
ISPATH(IC)  =  INODE (K) 

END  IF 

30  CONTINUE 

IF  (IC  .GT.  0  )  THEN 
DO  40  L  =  l.IC 
K  ■  IW0RK3(L) 

J  »  ISPATH(L) 

IF  (IWORKl(J))  THEN 
D  =  WK4(Ij  +  ARCOST(K) 

IF  (6  .LT.  WK4 ( J ) )  THEN 

WK4(J)  »  D 
IW0RK2( J)  *  K 
END  IF 

END  IF 

40  CONTINUE 

END  IF 

* 

*  (B).  FIND  THE  MINIMUH  POTENTIAL. 

* 


D  =  BIG 
I ENT  =  0 
I FIN  »  .FALSE. 

DO  SO  1=1.  N 

IF  ( IWORK1 ( I ) )  THEN 
IFIN  «  .TRUE. 

IF  (WK4 ( I )  .LT.  D)  THEN 
D  «  WK4 ( I ) 

IENT  a  i 
END  IF 
END  IF 

50  CONTINUE 

★ 

*  (C).  INCLUDE  THE  NODE  IN  THE  CURRENT  PATH. 

* 


IF  (D  .LT.  BIG)  THEN 

IWORK1 (iENT)a . FALSE . 

IF  (IENT  .NE.  LAST)  THEN 
I  «  IENT 
GO  TO  20 
END  IF 

ITT 

IF  (IFIN)  THEN 
NP  «  l 
RETURN 
END  IF 


ZJ)  THEN 


END  IF 

IJ  ■  LEST 
NUMP  -l 

ISPATH(l)  ■  LAST 

K  *  *»  <#M> 1J) " 

ELSE  ,  4 

ZJ-  INODE (K) 

END  IF 

NUMP  -  NUMP  +  1 
ISPATH(NUMP)  *  IJ 
IF  (IJ  .NE.  ISTART)  GO  TO  60 
L  *  NUMP  /  2 
J  »  NUMP 

DO  70  I  -  1 .L 

isPArafifn =I1spath(j> 


ISPATH(I)  »  I! 
ISPATH(J)  *  K 
J  *  J  -  1 
-  WK4(LAST) 


SPATH(J) 


XLEN  -  WK4(LAST) 

RETURN 

END 

*****************************  ****************************************** 
*  E.  SUBROUTINE  FOR  PRINTING  THE  RESULTS  OF  CALCULATION. 
********************* ************************************************** 
SUBROUTINE  PRNT 


INTEGER 
INTEGER 
CHARACTER 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
COMMON  / 


COMMON  / 

COMMON  / 
COMMON  / 
OPEN  (UNIT 
OPEN  (UNIT 


TAIL , HEAD , TIME , FLOW 

CASE , NTYPE , PARAM , SP , UNITF , ROWSP 

resh 

TAIL (150) ,HEAD(1S0) ,TIME(150) 
FLOW(IOO) .DIST (100) ,WIDTH(100) 


QSliSS  >i85raBjiSrww 

IN^DE(j$o|  !  JNO^E  ( 1 50  ) ,  ARCOST  ( 1 50  ) 

ISPATH(50) 

INTS  /  CASE,NTYPE,N,M,NARC, ISTART  LAST, PARAM, 
NUMP , UNITF , NVEH , NCOL , NROW , ROWSP , MINE , 
NWD, NDP , NFORWD 

INTA  /  TAll, HEAD, TIME, FLOW, LINK, LARC, NO, SP, 
INODE , JNODE, I SPATH 

REALS/  BIG , TPATH , XLEN , PATHL , FORWD , FORDP , PRT 
REALA/  WIDTH, DIST, RESM 
»  10,  FILE  *  'TIME' ) 

=  20,  FILE  *  'DIST') 


OPEN  (UNIT  -  20,  FILE  *  'DIST1 


*  (A).  PRINT 

* 


RESULTS  FOR  MINIMUM  TIME  PATH. 


THEN 


IF  (PARAM  .EO.  1)  THEN 
WRITE (  6,711) 

FORMAT( 5X^35 ('-' ),/,5X, 'MIN  TIME 

1F  (<WRITE?Q6,7l2)  NROW, NCOL, ROWSP 
WRITE (10, 7 12)  NROW, NCOL, ROWSP 
FORMAT ( 1  ',4X  'VEHICLE  UNIT' 


TIME  PATH' ,/,5X,15('-')) 


' '  'tSx,' VEHICLE  UNIT** ' .  / ,  5X ,  'FORMATION  . '  ,13 
ROW$ ' ' f 3 , '  COLUMNS ' , / , 17X, ' ROW  SPACE ' , 


'  Runs 

13,'  METERS 
ELSE  IF  (CASE  .EQ,  2)  THEN 
NWD  *  FORWD  *  1000 
NDP  -  FORDP  *  1000 
IF  (UNITF  .EQ.  1)  THEN 
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713 


714 


* 

* 


6,713)  N 

L0,713£  NROW , NCOL , NWD , NDP 


NROW , NCOL , NWD , MDP 
HOW, NCOL, NWD, NDP 
DISMOUNTED  TROOPS1 


JX 


Li?,/, 


WRITE  I 
WRITE! 

FORMAT A  _ 

' FORMATION  •  ROWS  • *15.  '  ,  COlnHWS  ■  t 
16X,'  WIDTH', 15.',  DEPTH', IS,'  METERS') 
ELSE  IF  (UNXTF  .EQ.  2)  THEN 
WRITE?  6,714)  NWD,NDP 
WRITE (10, 714)  NWD, NDP 
FORMAT ( '  ' , 4X , ' DISMOUNTED  TROOPS ' , /  ' SX , 
'FORMATION  s ' , ‘  WIDTH', IS,1,  DEPTH', 16,' 
END  IF 
END  IF 


721 


722 


IF  (RESM  .E 
WRITE 
WRITE 
FORMAT 


Y')  THEN 
"  MINE 


MINE 

MINE  FIELD  (  ARC  NUMBER 

/ , 5X,35( ' - ' )) 

ELSE  IF  (RESM  .EJ.  'N')  THEN 


METERS  1 ) 


M3,') 


END  IF 


WRITE (  S 
WRITE? 10 ,722) 

FORMAT (/,5X, ’NO  MINE  FIELD  ARC' ,/ ,5X,35( '  -  ' ) ) 


*  WHEN  PRT  EQUALS  0.0,  THERE  IS  A  FEASIBLE  ROUTE  TO  GO. 

IF  (PRT  .EQ.  0.0)  THEN 
IY  *  100 

ISUMT  a  XLEN  /  100 
NXLEN  =  XLEN 

MSUMT  =  (  MOD(  NXLEN,  IY  )  *  60  )  /  100 
WRITE?  6,731)  ISUMT,  MSUMT 
WRITE (10, 731)  ISUMT,  MSUMT 

731  FORMAT ( 1  ',/ ,SX, 'SUM- OF  TRAVERSAL  TIME  :  ', 

*  13, 2X, 'H0UR',I5,2X, 'MINUTE' ) 

ITHR  a  TPATH  /  100 
NTPATH  a  TPATH 
MINUT 


\ ,  tuiiu,  iiiuiuBndiui  .  ina  .  , 

13, 2X, 'HOUR' ,I5,2X, 'MINUTE' ) 
PATHL 
PATHL 


732 

WRITE (  6,732) 
WRITE (10, 732) 
FORMAT ( '  ',/, 

733 

* 

* 

WRITE ( 
WRITE (1 
FORMAT! 

6,733) 

.0,733) 

734 

WRITE ( 
WRITE (l 
FORMAT  ( 

6,734) 

,0,734 

735 

WRITE (  6,735) 
WRITE (10, 735) 
FORMAT (1 OX, 13 

736 

WRITE?  6,736) 
WRITE (10,736) 
FORMAT? 1  ' , /  . 

750 

WRITE?  6,750 
-  WRITE (10, 750) 
F0RMAT?1QX, 13 

760 

* 

WRITE (6, 760) 
FORMAT? '  ' , / , 

F7.2,2X, 'KM') 


I  al,  MUMP) 


nwi 

11: 


I  »1,  NUMP 


-i! 


F5.1) 

*  WHEN  PRT  EQUALS  1.0.  THERE  IS  NO  FEASIBLE  ROUTE  TO  GO. 

ELSE  IF  (PRT  .EQ.  1.0)  THEN 
WRITE?  6,741) 

WRITE (10, 741) 

741  FORMAT </,5X, 'THERE  IS  NO  FEASIBLE  ROUTE  TO  GO.') 

END  IF 


*  (B).  PRINT  THE  RESULTS  FOR  MINIMUM  DISTANCE  PATH. 
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811 

821 

* 

822 

831 

832 

833 

834 
860 

★ 

*  PRINT 

A 


WRITS  6,811 
WRITS (20 

TOShfi.-*' 

WRITS!  6,821 
WRITE(20,821 
FORMAT (/,SX, 


',5X, 'MIN 
THEN 


DISTANCE  PATH1  ,/,5X,19( '-' )) 


END 


MINS 
MINS 

.  1  MINS  FIELD  (  ARC  NUMBER  *  * ,13 , ' ) '  , 
/,5X,35( '-' )) 

ELSE  IF  (RESM  .EQ.  'N')  THEN 
WRITE?  6,322) 

WRITE v 20  822) 

FORMAT (/ *SX, 'NO  MINE  FIELD  ARC' ,/,5X,35( )) 

END  IF 

WRITE (  6,831) 

WRITE (20 ,831) 

FORMAT ( '  1 , / , SX, 'NODE  NUMBER  ALONG  MIN  DISTANCE  PATH  ' ) 
WRITE?  6,8i2  ?I,ISPATH(I),  I  *1,  NUMP) 

WRITE ?20, 832) (I , ISPATH(I) ,  I  *1,  NUMP) 
FORMAT(l6x,l3,3X,I3) 

WRITE?  6,833) 

WRITE (20  833) 

FORMAT (' /,SX, 1  ARC  NUMBER  ALONG  MIN  DISTANCE  PATH  •) 
WRITE?  6,8i4  ?I,LARCfI),  I  *1,  NUMP-1) 

WRITE ? 20 , 834 H I , LARC ( I ) ,  I  *1,  NUMP-1) 

FORMAT (10X, 13, 3X, 13) 

WRITE?  6,360)  XLEN 
WRITE (20  860)  XLEN 

FORMAT ( ' ' ' ,/,5X, ‘TOTAL  DISTANCE  ALONG  MIN  DIST  PATH  : ' , 
FS.l,''  KM') 

IF 


921 

920 


THE  INPUT  DATA  FOR  CHECKING. 

WRITE (6  921) 

WRIIEU:92oi  (I) 

IF  (PARAM  .EQ.  1)  THEN 
MU  *  10 

ELSE  IF  (PARAM  .EQ.  2)  THEN 
MU  ■  20 
END  IF 
WRITE (MU, 921) 

WRITE (MU, 920) 

DISf(I),WIDTH(i),SP(i: 

FORMAT ( 1  NO  TAIL  HEAD  TIME  FLOW  DIST' 

'  WIDTH  SPEED ' ) 

FORMAT ( 2X , I 3 , 21 5 , 2 I 6 , 3X , F4 . 1 , 3X , F6 . 3 , 2X , 14 ) 

RETURN 
END 


(NO(I) .TAIL(I) ,HEAD(I) ,TIME(l) ,FLOW(I) , 
DIST? I ) , WIDTH  ( I ) , SP ( I ) , 1*1 , MARC) 
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APPENDIX  E 

COMPUTER  EXEC  PROGRAM 


This  appendix  contains  two  computer  exec  programs  used  to  support  the 
algorithms  in  Chapter  111  and  IV. 

SINGLE  ARC  EXEC  PROGRAM 

This  exec  program  is  used  for  the  computer  program  for  single  arc  attributes. 

&TRACE  OFF 

&TYPE  Please  provide  the  FILENAME  for  your  VS  FORTRAN  program. 

&READ  VAR  &FN 

&TYPE  Do  you  need  to  compile  your  program  ?  (Y) 

&READ  VAR  &R  COMPILE 

&IF  &R  COMPILE  NE  Y  &GOTO  -RUN 

-H  FORfVS  &FN 

&IF  &RC  EQ  0  &SKIP  9 

&TYPE  Your  prograra  did  not  compile;  check  for  errors. 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP1 

&IF  &RESP1  NE  Y  &EXIT  1 

&COMMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP2 

&IF  &RESP2  EQ  Y  &GOTO  -H 

&EXIT  1 

-RUN  &TYPE  Do  you  wish  your  INPUT  to  be  from  the  terminal?  (Y) 

&READ  VAR  &IN 

&IF  ftIN  NE  Y  &GOTO  -RUN2 

-INPUT.FILE  FILEDEF  01  DISK  NB71  DATA  A1 

-INPUT.FILE  FILEDEF  02  DISK  NBN0D1  DATA  A1 

-INPUT_FILE  FILEDEF  03  DISK  NBARC1  DATA  A1 

-INPUT  FILE  FILEDEF  04  DISK  SPEED  DATA  A1 

FILEDEF  05  TERMINAL 

-RUN2  &TYPE  Do  you  wish  your  OUTPUT  to  go  to  the  terminal?  (Y) 

&READ  VAR  &0UT 

&IF  &0UT  NE  Y  &GOTO  -OUTPUT_FILE 
&GOTO  -LOAD 

-0UTPUT_FILE  FILEDEF.  06  DISK  &FN  OUTPUT  A  (LRE'U,  133 
-LOAD  LOAD  &FN  (START 
&IF  &RC  EQ  0  &SKIP  9 

&TYPE  Your  program  did  not  run  correctly;  check  for  errors. 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP3 

&IF  &RESP3  NE  Y  &EXIT  2 

&C0MMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP4 

&IF  &RESP4  EQ  Y  &GOT0  -H 

&EXIT  2 

&IF  &0UT  EQ  Y  &GOTO  -REDO 

&TYPE  Your  output  is  in  the  file  &FN  OUTPUT  A 
&TYPE  Do  you  wish  to  BROWSE  your  output?  (Y) 

&READ  VAR  &RESP 

&IF  &RESP  EQ  Y  &COHMAND  BROWSE  &FN  OUTPUT  A 
&TYPE  Print  your  output  file?  (Y) 

&READ  VAR  &RESP7 

&IF  &RESP7  EQ  Y  &C0MMAND  PRINT  &FN  OUTPUT  A 
-REDO 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y/N; 
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filtKJU)  VAR  ARESP5 

AIF  AHESF5  EQ  V  XEDXT  AFN  FORTRAN  A 

ATYPE  Do  you  wish  to  run  th«  program  again?  (T) 

AREAD  VAR  &RESP6 

AIF  ARESP6  E0  7  AND  AR1SP5  EQ  7  AGOTO  -H 

AIF  ARESP6  EQ  7  AOOTO  -RUN 

AEXIT 

NETWORK  EXEC  PROGRAM 

This  exec  program  is  used  for  computer  program  for  Cartesian  space  network. 

&TRACB  OFF 

&TYPE  Please  provide  the  FILENAME  for  your  VS  FORTRAN  program. 

&READ  VAR  &FN 

AT7PE  Do  you  need  to  compile  your  program  ?  (7) 

AREAD  VAR  &R COMPILE 

AIF  AR  COMPILE  NE  7  AOOTO  -RUN 

-H  FORTVS  AFN 

AIF  ARC  EQ  0  ASKIP  9 

ATYPE  Your  program  did  not  compile;  check  for  errors. 

ATYPE  Do  you  wish  to  XEDIT  the  program  file?  (7) 

AREAD  VAR  ARESP1 

AIF  ARESP1  NE  Y  AEXIT  1 

ACOMMAND  XEDIT  AFN  FORTRAN  A 

ATYPE  Do  you  wish  to  run  the  program  again?  (7) 

AREAD  VAR  ARESP2 

AIF  ARESP2  EQ  Y  AGOTO  -H 

&EXIT  1 

-RUN  ATYPE  Do  you  wish  your  INPUT  to  be  from  the  terminal?  (Y) 

AREAD  VAR  AIN 

AIF  AIN  NE  Y  AGOTO  -RUN 2 

-TNPUT_FILE  FILEDEF  01  DISK  NETA  DATA  A1 
-INPUT  FILE  FILEDEF  02  DISK  NETB  DATA  A1 
FILEDEF  05  TERMINAL 

-RUN2  ATYPE  Do  you  wish  your  OUTPUT  to  go  to  the  terminal?  (Y) 

AREAD  VAR  AOUT 

AIF  AOUT  NE  Y  AGOTO  -OUTPUT_FILE 
AGOTO  -LOAD 

-OUTPUT  FILE  FILEDEF  06  DISK  AFN  OUTPUT  A  (LRECL  133 
-LOAD  L5AD  AFN  (START 
AIF  ARC  EQ  0  ASKIP  9 

ATYPE  Your  program  did  not  run  correctly.*  check  for  errors. 

ATYPE  Do  you  wish  to  XEDIT  the  program  rile?  (Y) 

AREAD  VAR  ARESP3 

AIF  ARESP3  NE  Y  AEXIT  2 

ACOMMAND  XEDIT  AFN  FORTRAN  A 

ATYPE  Do  you  wish  to  run  the  program  again?  (Y) 

AREAD  VAR  ARESP4 

AIF  ARESP4  EQ  Y  AGOTO  -H 

AEXIT  2 

AIF  AOUT  EQ  Y  AGOTO  -REDO 

ATYPE  Your  output  is  in  the  file  AFN  OUTPUT  A 
ATYPE  Do  you  wish  to  BROWSE  your  output?  (Y) 

AREAD  VAR  ARESP 

AIF  ARESP  EQ  Y  ACOMMAND  BROWSE  AFN  OUTPUT  A 
ATYPE  Print  your  output  file?  (Y) 

&READ  VAR  &RESP7 

AIF  ARESP7  EQ  Y  ACOMMAND  PRINT  AFN  OUTPUT  A 
-REDO 

ATYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 

AREAD  VAR  ARESP5 

AIF  ARESP 5  EQ  Y  XEDIT  AFN  FORTRAN  A 

ATYPE  Do  you  wish  to  run  the  program  again?  (Y) 

AREAD  VAR  ARESP6 

AIF  ARESP6  EQ  7  AND  ARESP5  EQ  Y  AGOTO  -H 

AIF  ARESP6  EQ  7  AGOTO  -RUN 

AEXIT 
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APPENDIX  F 

UNIT  TYPE  AND  FORMATION 


Table  18  shows  the  integer  value  and  the  unit  type  it  represents. 


TABLE  18 
TYPE  OF  UNIT 


Integer  Value 

Type  of  Unit 

1 

Tracked  Vehicle 

2 

Wheeled  Vehicle 

3 

Dismounted  Troops 

|  Table  19  shows  the  integer  value  and  the  unit  formation  it  represents. 

1 


TABLE  19 

UNIT  FORMATION 


Integer  Value 

Unit  formation 

1 

Deployment 
(  multiple  column  ) 

2 

Undep.loyment 
(  single  column  ) 
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APPENDIX  G 
RATING  CONE  INDEX 


Figure  G.l  shows  an  example  of  bow  the  rating  cone  index  can  be  used  [Ref.  4;  p. 

20]. 
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